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ANNOUNCEMENT GR VOLUME OF ARS P 


Major space age role for 
fluorine rocket fuel oxidizers 


Potential applications for fluorine in rocketry extend 
over the entire propellant spectrum—from cryogenics 
to storable liquid to solid propellants 


As America’s primary producer of fluorine and halogen fluorides 
and leader in fluorochemical research, General Chemical helps sup- 
port America’s rocket and missile development programs in three 
key areas. 

For investigation of cryogenic systems in the super-high energy 
range, General supplies large quantities of liquid fluorine. For inves- 
tigation of advanced, high-energy storable liquid propellant systems, 


¥ te General supplies chlorine trifluoride (CTF) and other halogen 
25,000 POUNDS OF LIQUID FLUORINE, from General's fluorides. To serve both these areas, General has large existing pro- 
Metropolis, Ill, Works — ready for over-the-highway duction capacity, and can meet volume requirements. 


shipment in tank trucks of our own design to rocket 
and missile research centers throughout the nation. 


In solid propellant systems, General is conducting intensive inves- 
tigation under Government contract. Incorporation of the fluorine 
atom in the oxidizing molecule is part of this research. 

If you would like technical data, we will be glad to send you our 
free bulletins, “Fluorine” and “Chlorine Trifluoride.” Simply request 
them on company letterhead. 


CHLORINE TRIFLUORIDE has been produced commercially 
by General Chemical for the past 10 years. Reactors 
like this one, located at General’s Baton Rouge Works, 
have a capacity of more than 100 tons of CTF per year. 


PART OF THE GIANT FLUORINE CELL ROOM at General Chemical’s Metropolis, Ill, 
Works. This is the largest privately owned fluorine production facility in the world. 


: First in 
GENERAL CHEMICAL SCIENTISTS make use of the most ads fluorine GENERAL CHEMICAL DIVISION 
vanced scientific equipment to investigate fluorine com- chemistry 40 Rector Street, New York 6, N. Y. 


pounds as high energy oxidizers for solid propellants. 


For 


OcroBE 


= | FORECAST | 
| 
: 


For data on the environment of space, the 


JOURNAL OF GEOPHYSICAL RESEARCH 


is a most important source of basic information. 
The Journal of Geophysical Research publishes more 
articles containing basic data on the environment 
provided by the Earth, its atmosphere, and the space 
beyond than any other journal. In addition to basic 
scientific data on the environment, it contains articles 
on the relevant physics. For example, during 1959 
papers were published on the following subjects: 
Van Allen radiation, cosmic radiation, atmospheric 
density, temperature and composition at satellite 
altitudes, magnetic fields in space, solar radiation, 
and many other subjects of value in systems design. 


Papers in other areas traditional to the broad field of 
geophysics are also included. In 1959 the Journal 
contained over 2400 pages; it is anticipated that sub- 
stantially more material will be published in 1960. 


Groups of papers presented at symposia are frequently 
and promptly published in the Journal, making the 
summaries and opinions of leading scientists quickly 
available in print. In 1959, the Journal contained the 
proceedings of symposia on Scientific Effects of 
Artificially Introduced Radiations at High Altitudes 
(Argus results), Exploration of Space, and Fluid 
Mechanics in the Ionosphere. 


Monthly—$20.00 for calendar year 1960 on subscription 


AMERICAN GEOPHYSICAL UNION 


Department M2 
1515 Massachusetts Avenue, N. W., Washington 5, D. C. 


Other publications of the AGU... 


TRANSACTIONS: quarterly, including IGY Bulletin, $4.00 for the calendar year 
GEOPHYSICAL MONOGRAPH SERIES: No. 5 now available, Physics of Precipitation, 435 pp., $12.50 
IZVESTIYA, Academy of Sciences, USSR, English edition: monthly, available for 1957, 1958, 1959, and 1960, 


$25.00 for each calendar year 


GEODEZIYA I KARTOGRAFIYA (Geodesy and Cartography): monthly, 1959 series, $20.00 
TRUDY, Geophysics Institute, Academy of Sciences, USSR: Nos. 37-40 (in cooperation with Consultants 


Bureau), prices range from $6.00 to $8.00 


OcroBER 1960 


The Task — Provide celestial accuracy for massive sensing devices. 
The Solution — Hydrostatic bearings for smooth, rigid, 
“frictionless” support. 


In Polaris submarines and other military applications, these 
hydrostatic bearings, developed by Nortronics, are particularly suited to 
virtually frictionless support of a large structure, which can be 
positioned accurately by a servo system with no unevenness, runout or 
binding — yet extremely stiff under applied loads. 

A NOTE FOR TALENTED ENGINEERS: 

We invite your inquiry as to important positions at Nortronics’ 
Precision Products Department. Excellent salaries for qualified men 
in both our Systems and Component Groups. 


Contact E. P. Reardon 


Ti PRECISION PRODUCTS DEPARTMENT 
NOR RONIGS NORWOOD, MASSACHUSETTS 
Field Offices 
Highway #46 2486 Huntington Drive 


A D IviSIOn of Teterboro, New Jersey San Marino, California 
Telephone: ATlantic 7-0461 


NORTHROP CORPORATION Heights 871-0 TWX-Alhambra 9619-U 
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AEROSPACE 
DICTIONARY 


Your complete 


guide to the 
new language 
of space 


exploration 


by Frank Gaynor 


with an introduction 
by Wernher von Braun 


As man probes further into space, a 
vital new language comes swiftly 
into being. This up-to-the-minute 
reference work answers the needs 
of the student and general reader, 
as well as those in government and 
industry who require a knowledge 
of the essential terminology in 
space exploration. 


Thousands of clear, concise en- 
tries, alphabetically arranged, 
present the most authoritative infor- 
mation on the language of rocketry 
and astronautics, guidance systems, 
satellites, telemetering devices, 
manned space flight, re-entry phe- 
nomena, celestial mechanics and 
aeronautic science. 


The dictionary’s editor, Frank 
Gaynor, had the assistance of the 
United States Department of De- 
fense as well as civilian agencies 
and research centers in preparing 
the book’s entries. America’s most 
distinguished authority on rocketry 
and space travel, Dr. Wernher von 
Braun, has also assisted in the 
work’s preparation and contributed 
a provocative introduction “The 
Why of Space Travel.” 


AD-7 
PHILOSOPHICAL LIBRARY, Publishers 
15 East 40th Street, New York 16, N. Y. 


Please send me .... copies of AEROSPACE 
DICTIONARY @ $6.00 per copy. To expedite 
shipment | enclose 


Check Money Order 
(0 Check here if you wish C.O.D. 
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_ ASSOCIATION WITH LEADING SCIENTISTS 
_ GRADUATE THESIS PROGRAM 


For employment 

write to 

Personnel Director 
Division 60-91 


scientific laboratory 


OF THE UNIVERSITY OF CALIFORNIA 
LOS ALAMOS, NEW MEXICO 


Scope of ARS JOURNAL 


This Journal is devoted to the advance- 
ment of astronautics through the dissemina- 
tion of original papers disclosing new scientific 
knowledge and basic applications of such 
knowledge. The sciences of astronautics are 
understood here to embrace selected aspects 
of jet and rocket propulsion, spaceflizht 
mechanics, high speed aerodynamics, flizht 
guidance, space communications, atmospheric 
and outer space physics, materials and striic- 
tures, human engineering, overall syst-m 
analysis, and possibly certain other scient ific 
areas. The selection of papers to be prinied 
will be governed by the pertinence of the topic 
to the field of astronautics, by the current or 
probable future significance of the research, 
and by the importance of distributing the in- 
formation to the members of the Society «nd 
to the profession at large. 


Information for Authors 
Manuscripts must be as brief as the pro er 


presentation of the ideas will allow. !°x- 
clusion of dispensable material and con: se- 
ness of expression will influence the Edit:rs’ 
acceptance of a manuscript. In terms of 
standard-size double-spaced typed pages a 
typical maximum a is 22 pages of ‘ext 
(including equations), age of referer: es, 
1 page of abstract ants illustrations. 
Fewer illustrations permit text, and ice 
versa, Greater length will be accept:ble 
only in exceptional cases. 

Short manuscripts, not more than one 
quarter of the maximum length stated for ‘ull 
articles, may qualify for publication as 
Technical Notes or Technical Comme: nts, 
They may be devoted to new developments 
requiring prompt disclosure or to communts 
on previously published papers. Such manu- 
scripts are published within a few months o 
the date of receipt. 

Sponsored manuscripts are publisied 
occasionally as an ARS service to the indus- 
try. A manusc-ipt that does not qualify for 
publication, according to the above-stated 
requirements as to subject, scope or lenyth, 
but which nevertheless deserves widespread 
distribution among jet propulsion enginecrs, 
may be printed as an extra part of the Journal 
or as a special supplement, if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as a 
footnote on the first page of the article. 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be donkte spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 
200 word abstract. State the authors’ posi- 
tions and affiliations in a footnote on the 
first page. Equations wed symbols may be 
handwritten or typewritten; clarity for the 
printer is essential. Greek letters and unusual 
symbols should be identified in the margin. If 
handwritten, distinguish between capital and 
lower case letters, and indicate subscripts and 
superscripts. References are to be grouped at 
the end of the manuscript and are to be given 
as follows. For journal articles: Authors 
first, then title, journal, volume, year, page 
numbers; for books: Authors first, then title, 
publisher, city, edition and page or chapter 
numbers. Line drawings must be clear and 
sharp to make clear engravings. Use black 
ink on white paper or tracing cloth. Lettering 
should be large enough to be legible after 
reduction. Photographs should be glossy 
prints, not matte or semi-matte. Each illus- 
tration must have a legend; legends should be 
listed in order on a separate sheet. 

Manuscripts must be accompanied by 
written assurance as to security clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility reste 
with the author. 

Preprints of papers presented. at 


meetings are aut lly ed for 
publication. 
Sub it iptsind li (origi- 


nal plus first carbon, with two sets of 
illustrations) to the Managing Editor, ARS 
500 Fifth Avenue, New York 


ARS JOURNAL is published monthly by 
the American Rocket Society, Inc. and . 
American Interplanetary Society at Seth & 
Northampton Sts., Easton, Pa., U. A. 
Editorial offices: 500 Fifth Ave., New . ork 
36, N. Y. Price: $12.50 per year, $2.00 
per single copy. Second-class postage paid 
at Easton, Pa., with additional entry at New 
York, N. Y. This publication is author- 
ized to be mailed at the special rates 
of postage prescribed by Section 132.122. 
Notice of change of address should be sent 
to the Secretary, ARS, at least 30 days prior 
to publication. Opinions expressed herein 
are the authors’ and do not necessarily refiect 
the views of the Editors or of the Socicty. 
© Copyright 1960 by the American Rocket 
Society, Inc. 
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Builders of more large, thin wall, high 
strength solid propellant rocket engine 
cases and nozzles for development pur- 
poses than any other company in America. 


A small experienced organization geared 
to handle your development and prototype 
requirements for static and flight tests in 
the shortest possible time. 
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New Books from ACADEMIC PRESS 


INTRODUCTION 
TO ASTRODYNAMICS 


By Rosert M. L. Baker, Jr. and Maup W. MakEMSON 
University of California at Los Angeles 


September 1960, 358 pp., illus., $7.50 


This is the first introductory textbook based on a university 
course in astrodynamics. In concise, clear language it de- 
scribes all of the practical problems and methods encountered 
in the field. The nomenclature is carefully defined in a glos- 
sary of over three hundred terms. 

The book will serve as a textbook for introductory under- 
graduate courses and as a reference and handbook of astro- 
dynamical techniques for engineers and scientists involved 
in space technology. It differs from classical texts in celestial 
mechanics by relating theory and techniques directly to con- 
temporary space vehicle problems. 


HYPERSONIC FLOW THEORY 


Wattace D. Hayes 
Princeton University and Space Technology Laboratories 
Ronap F. PRosBsTEIN 
Brown University and Avco Research Laboratory 


1959, 464 pp., illus., $11.50 


“.. the nature of the present volume makes it an im- 
portant contribution to the contemporary scene as well as 


a book that will remain vital for many future years.” 
—ARS Journal 


“Those readers familiar with the other work of the authors 
will not be surprised to note that, in presenting a comprehen- 
sive picture of hypersonic flow theory, Prof. Hayes and Prob- 
stein have also succeeded in presenting some aspects of general 
compressible flow theory in a new and stimulating light. 

“*.. . it will long be regarded as the authoritative text on the 
theoretical foundations of the subject. It can be sirongly re- 
commended to anyone directly concerned with problems of 
high speed gas dynamics.” —Aircraft Engineering 


SPACEFLIGHT TECHNOLOGY 


The First Commonwealth Spaceflight Symposium 
Organized by the British Interplanetary Sociely, 1959 


Edited by KENNETH W. GATLAND 
The British Interplanetary Society 


November 1960, aboul 400 pp., illus., approx. $12.00 


The imporiant spaceflight research in progress throughout 
the British Commonwealth is thoroughly discussed by aciive 
participanis in space science from Woomera to Jodrell Bank. 
The subjects covered include equipment necessary for manned 
satellites, satellite tracking, space navigation, satellite re- 
covery, nuclear thermal fission rockets, scientific instrumenta- 
tion, and heaiing problems of entry into planetary atmos- 
pheres. 


Advances in SPACE SCIENCE 


Edited by Freperick I. Orpway, IIT 
George C. Marshall Space Flight Center 
National Aeronautics and Space Administration 


Editorial Advisory Board 


WERHNER VON Braun, FREDERICK C. Durant, III, 
Eucen SANGER, R. SHEPHERD, GEORGE 
Sutton, and ETrENNE VaAssy 


“the scientists who function in this field must be com- 
petent specialists who have a moderately complete knowledze 
of the subject as a whole. Advances in Space Science should 
be extremely helpful to specialists for securing information 
concerning recent contributions made in their specialty to 
astronautics as well as for obtaining a general picture of the 
subject as a whole.” —Science (Review of Volume |) 


VOLUME 1 1959, 412 pp., illus., $12.00 


VOLUME 2 1960, 450 pp., illus., $13.00 
Experimental Physics Using Space Vehicles, 
by P. SonettT 
Tracking Artificial Satellites and Space Vehicles, 
by Karu G. HENIZzE 
Materials in Space, by Freperick L. BacBy 
Plasma Propulsion Devices, by Morton CaMac 
Electrostatic Propulsion Systems for Space Vehicles, 
by Ernst STUHLINGER and Rosert N. SEItTz 
Attitude Control of Satellites and Space Vehicles, 
by Rosert E. Roperson 
AUTHOR INDEX—SUBJECT INDEX. 


VOLUME 3 1961, in preparation 


SPACE TRAJECTORIES 


A Symposium Sponsored by 
The American Astronautical Society, 
The Advanced Research Projects Agency, 
and Radiation Incorporated 


Chairman of the Symposium: 
T. C. Hetvey, Radiation Incorporated 


Seplember 1960, 298 pp., illus., $12.00 


This symposium was planned to achieve a comprehensive 
survey of the field and the widest possible exchange of ideas 
between the leading scientists who participated. The area 
to be covered by each paper was carefully outlined before- 
hand in order to provide a truly comprehensive treatment of 
the entire field. This volume, therefore, will be useful to 
specialists in the field and will serve as a valuable reference 
book for those working in related sciences. 


ACADEMIC P RESS, New York and London 


111 Fifth Avenue, New York 3 
926 


17 Old Queen Street, London, S.W. 1 
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ACADEMIC PRESS, Publishers 


are pleased to announce a new serial publication 


Progress in 


Astronautics and Rocketry 


Sponsored by the American Rocket Society 


Series Editor: Martin SUMMERFIELD, 


Princeton University 


With this Progress series, the American Rocket Society is 
inaugurating a new approach to the basic task of publishing 
the scientific output of an association with many diverse 
specialties. Each Progress volume will be focused on a 
single specialty. 


Volume: Solid Propellant Rocket Research 


A Selection of Technical Papers based mainly on a 
Symposium of the American Rocket Society held at 
Princeton, New Jersey, January 28, 29, 1960 


MECHANICAL PROPERTIES OF GRAINS 


A. B. Boyd, W. M. Burkes, and J. E. Medford, Grain Design and 
Development Problems for Very Large Rocket Motors 

A. M. Freudenthal, and L. A. Henry, On ‘“‘Poisson’s Ratio” in 
Linear Viscoelastic Propellants 

M. L. Williams, Mechanical Properties and the Design of Solid 
Propellant Motors 

N. N. Au, A Method of Strength Analysis of Solid Propellant 
Rocket Grains 

J. Vandenkerckhove, and G. Lampens, Stress and Strain Analysis 
of Cylindrical Case-Bonded Grains 


STEADY-STATE BURNING MECHANISMS 


Martin Summerfield, G. S. Sutherland, M. J. Webb, H. J. Taback, 
and K. P. Hall, Burning Mechanism of Ammonium Per- 
chlorate Propellants 

D. W. Blair, E. K. Bastress, C. E. Hermance, K. P. Hall, and 
M. Summerfield, Some Research Problems in the Steady- 
State Burning of Composite Solid Propellants 

W. Nachbar, A Theoretical Study of the Burning of a Solid 
Propellant Sandwich 

R. F. Chaiken, and W. H. Andersen, The Role of Binder in 
Composite Propellant Combustion 


COMBUSTION OF METALS 


Irvin Glassman, Combustion of Metals: Physical Consider- 
ations 

W. M. Fassell, C. A. Papp, D. L. Hildenbrand, and R. P. Sernka, 
The Experimental Nature of the Combustion of Metallic 
Powders 

Derck A. Gordon, Combustion Characteristics of Metal Particles 

Claude P. Talley, The Combustion of Elemental Boron 

William A. Wood, Metal Combustion in Deflagrating Propellant 


Edited by Martin SUMMERFIELD 
October 1960, 692 pp., illus., $6.50. 


THEORIES OF UNSTABLE COMBUSTION 


F. T. McClure, R. W. Hart, and J. F. Bird, Solid Propellant 
Rocket Motors as Acoustic Oscillators 

Reuel Shinnar, and Menachem Dishon, Heat Transfer Stability 
Analysis of Solid Propellant Rocket Motors 

A. G. Smith, A Theory of Oscillatory Burning of Solid Propellants 
Assuming a Constant Surface Temperature 

Sin-I] Cheng, Combustion Instability in Solid Rockets Using 
Propellants with Reactive Additives 

R. W. Hart, J. F. Bird, and F. T. McClure, The Influence of 
Erosive Burning on Acoustic Instability in Solid Propellant 
Rocket Motors 


EXPERIMENTS ON UNSTABLE BURNING 


W. Grant Brownlee, and Frank E. Marble, An Experimental 
Investigation of Unstable Combustion in Solid Propellant 
Rocket Motors 

E. M. Landsbaum, W. C. Kuby, and F. W. Spaid, Experimental 
Investigations of Unstable Burning in Solid Propellant 
Rocket Motors 

T. A. Angelus, Unstable Burning Phenomenon in Double-Base 
Propellants 

E. W. Price, Review of Experimental Research on Combustion 
Instability of Solid Propellants 

R. H. Wall, Resonant Burning of Solid Propellants: Review of 
Causes, Cures and Effects 


SOLID PROPELLANT IGNITION 


R. F. McAlevy, III, P. L. Cowan, and M. Summerfield, The 
Mechanism of Ignition of Composite Solid Propellants by 
Hot Gases 

A. D. Baer, N. W. Ryan, and D. L. Salt, Propellant Ignition by 
High Convective Heat Fluxes 

Rodney B. Beyer, and Norman Fishman, Solid Propellant Igni- 
tion Studies with High Flux Radiant Energy as a Thermal 
Source 


In preparation: 

Volume II: Liquid Rockets and Propellants 

Edited by Loren E. Bollinger, Martin Goldsmith, and Alexis W. 
Lemmon, Jr. 


Volume III: Space Power Systems 

Edited by Nathan W. Snyder 

Volume IV: Electrostatic Propulsion 

Edited by David Langmuir, Ernst Stuhlinger, and J. M. Sellen 


ACADEMIC P RESS, New York and London 


111 Fifth Avenue, New York 3 
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The Image Furnace as a 
Research Tool 


HE CONCENTRATION of radiant energy by optical 

means has been known since antiquity. In Aristophanes’ 
comedy “The Clouds,” Stepsiades explains to Socrates how he 
intends to defeat his creditors by aiming a burning glass prop- 
erly as he meets them, so as to melt the wax tablets on which 
the debt is recorded (1).! 

The first military use of this principle is attributed to 
Archimedes during the siege of Syracuse in 212 B.C. He is 
reputed to have “set fire to the invading ships of Marcellus’ 
Navy by means of a burning glass composed of small square 
mirrors moving every way upon hinges, which, when placed 
in the sun’s rays, were directed upon the Roman fleet, so as to 
reduce it to ashes at the distance of a bow shot” (1). Archi- 
medes’ activities in defense of his home city are probably 
also the first case of the use of science to devise new weapons 
for the armed forces, and resulted in prolongation of the 
campaign for 3 years. 

Nearly two millennia were to pass before the concept was 
to be used for research purposes. In 1695, Averoni and Tar- 
gioni of Florence used a burning glass to decompose a diamond, 
a substance previously considered inalterable (2). During the 
next 50 years, Tschirnhauss in Germany, Parker in England, 
and Geoffroy and Homberg in France found the lens (or solar 
image) furnace an invaluable tool in their studies of the proper- 
ties of matter and fre (2-4). Its use as a high temperature 
tool of the new science of chemistry was climaxed by 
Lavoisier’s construction of a large double lens furnace in 1774 
to demonstrate that diamond was just another form of carbon 
(5). 

One of the laboratory demonstrations used by Tyndall dur- 
ing his lecture on the properties of radiant energy was prob- 
ably among the first examples of the arc image furnace (6). 
Fig. 1 illustrates the ignition of a chlorine-hydrogen mixture 
by an electric are through a double parabolic mirror system. 

The preceding application of radiant energy concentra- 


Received Sept. 1, 1960. 
! Numbers in parentheses indicate References at end of paper. 


NEVIN K. HIESTER and 
R. E. DE LA RUE 


Stanford Research Institute 
Menlo Park, Calif. 


tion exemplifies one of the important advantages of im:ge 
furnaces, i.e., their ability to supply uncontaminated heat at a 
point remote from the source. This, coupled with the fact 
that the materials being studied can be enclosed in a traiis- 
parent envelope and heated rapidly either in a vacuum or in a 
precisely controlled atmosphere without crucible contamina- 
tion, has led to the recently increased interest in such devices. 

Modern solar furnace research probably started in Germany 
in 1909, and sporadic activity, at about 10-year intervals, lias 
continued since then in Germany, the United States and 
France (2). A more intensive effort has been under way at 
several laboratories in recent years. 

Some completed research work has been reported by Trombe 
and Foex and their co-workers at Mont Louis in France (7) and 
by Mii’s group in Japan (8), but very little actual results have 
been published in this country. 

This has probably been due to several factors, the most im- 
portant of which is that expected high temperatures were not 
attained. This is possibly because many of these furnaces 
are built around war surplus searchlight mirrors whose optical 
qualities leave something to be desired, and because of the 
added complexity of a sun tracking mechanism whose con- 
tinued precise performance is a must. In addition, the 
vagaries of weather and the requirement that the sun not be 
too near the horizon also limit the working period. 

For these reasons emphasis has recently shifted to arc image 
furnaces where one has a relatively stationary source of radiant 
energy, and where the device can be used at any time and yet 
provide a flux and irradiation area comparable to that from a 
solar image furnace based on searchlight mirrors. For ex- 
ample, Stanford Research Institute’s solar furnace with an 
unusally good paraboloid of a 60-in. aperture and a heliostat 
sun follower can, under the best weather conditions, provide 
an irradiance of 160 cal/em?2/see over a 6-mm diameter focus 
(9). At the same time, the SRI are image furnace using a 
double ellipsoidal set of mirrors of 16-in. aperture can supply a 
flux of 160 cal/em?/sec and a focal diameter of 9 mm, with an 
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input of near 9 kw (10). The 160 cal/cm?/sec radiant fiux 
corresponds to a black-body temperature of 3250 K. As a 
result the trend toward are image furnaces as research tools has 
recently accelerated. 

The purpose of this paper is to provide a guide to the litera- 
ture covering some of the important theoretical and practical 
aspects of image furnaces and their operation. In addition, 
some of the research areas where these devices have been used 
will be briefly described in the hope that it may trigger ideas 
for other ways in which they can be used to solve research 
problems. 


Theory 


In an image furnace, regardless of whether the radiating 
source is the sun, the crater of a carbon arc or a luminescent 
flame, the optical system does no more than collect as much of 
the emitted energy as possible and focus it in a concentrated 
form in a single region. The furnace itself is at the plane of 
maximum radiant flux per unit area. The prediction of the 
size of the focal spot and of the flux is the goal of image fur- 
nace theory. The optics used govern the theoretical degree of 
concentration, but the actual concentration depends on many 
practical factors. 

The majority of solar furnaces in the United States are 
single continuous paraboloidal mirrors, either following the 
sun directly or through a heliostat. The theory for these 
types of furnaces has been well worked out (2, 11-18), includ- 
ing the effects of such sources of thermal losses to the radiant 
energy as: Atmospheric absorption, glass absorption, imper- 
fect reflectivity of reflecting surface, geometrical imperfec- 
tions of reflecting surface, shadowing effects and imperfec- 
tions in reflector alignments. 

These theories show, disregarding inefficiencies, that the 
highest flux in the focal plane is a function of paraboloid aper- 
ture to focal length ratio (or opening expressed in terms of the 
rim angle 6) and covers an image whose diameter is a function 
of the focal length. Thus, for a sun of uniform flux intensity, 
the maximum theoretical flux is 46,100 times sin? 6 times inci- 
dent solar flux, and the image diameter d is 0.0093 times focal 
length. The variation of flux, outside the image, in the focal 
plane has been calculated by Le Phat Vinh (17), Mii et al. 
(16) and Simon (19); and the first two have also calculated 
flux variations in planes parallel to the focal plane. 

In actuality the sun is brighter in the center—the so-called 
limb-darkening effect—and this must be taken into account. 
De La Rue et al. (20), Mii et al. (16) and Jose (21) have cal- 
culated the flux through the focal spot for the nonuniform 
case, and the first two have also considered the distribution for 
planes adjacent to the focal plane. 

As pointed out earlier, there are factors which contribute toa 
lower performance of a solar furnace than theory would 
predict. This design point has been considered seriously by 
Hiester et al. (2, 14, 22) and Bliss (13, 15). 

In connection with the large solar furnace proposed by the 
Air Force for construction at Cloudcroft, N. Mex., a number 
of theoretical studies were undertaken for the Air Force 
Missile Development Center at Holloman AFB. Several of 
these were concerned with the effect on performance of using 
mosaics of spherical (or toroidal) mirror segments to approxi- 
mate a paraboloid (18,23,24). This was necessitated by the 
impracticality of fabricating a continuous, geometrically per- 
fect, curved surface of the aperture desired. Baum and 
Strong (25) considered another approach to the system, 
reduction of the aperture by finding an optical system that 
places more of the intercepted radiant energy in the central 
sun image on the focal plane. They show that double mirror 
systems, having continuously curved surfaces (but where the 
geometry of each changes from point to point, in a related 
manner), do provide some small improvement.? The other 


_* This, of course, would be at the expense of increased fabrica- 
tion difficulties. 
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Fig. 1 Tyndall demonstrating an arc image furnace. 
(Redrawn by D. Lew) 


key paper in this series is by Jose (26) covering the theory for 
matching the heliostat to its paraboloid. 

Although, as pointed out earlier, are image furnaces are 
providing considerably more experimental results, to date the 
only theoretical discussion of their performance is that of 
Null and Lozier (27).8 


Solar Image Furnace 


The types of solar furnaces can be classified in a number of 
ways: 


Optics 
Single paraboloid 
Double paraboloid with focusing lens 
Double lens 
Paraboloid with ellipsoid 
Combined variable mirrors 
Cylindro-paraboloid 
Concentrating Surface 
1 Continuous geometric curved surface 
2 Flat segments tangent to the surface 
3 Curved segments approximating and tangent to the 
geometric surface 
4 Oriented flat segments (analogous to a Fresnel lens) 
Energy Collection 
a Concentrator movable, follows sun directly 
b Concentrator fixed, heliostat follows sun 


Thus the French solar furnace shown in Fig. 2 might be 
classified as type A3b. 

Most of the high temperature furnaces built recently are 
continuous paraboloidal mirrors, either following the sun 
directly or through a heliostat, types Alaand Alb. The only 
important exceptions are the Cla lens-type furnace at 


3 Additional work in this area is under way at Stanford Research 
Institute. 
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Fig. 2 General view of the 35-ft furnace located at Mont Louis 
in the Pyrenees 


California Institute of Technology, which utilizes an array of 
19 lenses that ultimately concentrate the radiations at a 
unique focus (28), the type A4a segmented flat plate furnace 
at Massachusetts Institute of Technology (29), and the type 
A3b segmented spherical concentrator furnace at the Quarter- 
master Research and Development Laboratories (30). 

Usually, in large units auxiliary mirrors are used to facilitate 
access to a stationary focus (30,31). Auxiliary mirrors in this 
ase are usually simpler to operate and actually less expensive 
than the delicate and sturdy mechanisms that must be pro- 
vided to rotate accurately the entire paraboloid and focus. 

Considerable effort has been extended in various design 
studies for solar furnaces in the past few years. Most of these 
design studies have been in connection with the large solar fur- 
nace proposed by the Air Force for construction at Cloudcroft, 
N. Mex. These studies include such topics as optical design 
principles (15), basic optical considerations in the choice of a 
design for a solar furnace (25), the loss of energy by absorption 
and reflection in the heliostat and parabolic condenser (32), a 
solar furnace using a horizontal heliostat array (33), design of 
a solar furnace and the calculations of concentration using 
spherical-mirror elements (23), thermal effects on solar fur- 
nace mirrors (34), design of a heliostat mirror for a solar 
furnace (26), design of the condenser of a solar furnace using 
nonparabolic elements (18), and the economic factors in fur- 
nace design (22). Other studies concerning mechanical fac- 
tors and operating aspects are contained in (35 and 36). 

The choice of reflector materials is governed by reflectivity 
and corrosion resistance, as well as by the optics andthe spec- 
trum desired. A discussion of reflectivity of metals and the 
corrosion characteristics has been given elsewhere (35,37,38). 
Those in use today are silver, aluminum, rhodium and stellite 
surfaces. Back-silvered glass has been determined to be the 
most favorable. The spectral distribution for the Stanford 
Research Institute solar furnace has been measured by 
Marcus and Wohlers (9). 

As important as the optical design of a solar furnace is the 
guiding system, or sun tracking mechanism. When a solar 
furnace is used for high temperature research, it is essential 
that the position of the hot zone on the target be very stable. 
This can be achieved by keeping the concentrator or heliostat 
focused on the sun. Manual focusing results in jerky motion 
and is therefore unsatisfactory. Mechanical rotation has to 
be guided so as to follow precisely the apparent motion of the 
sun. Various guiding systems have been described in the 
literature (8,39,40,41). 

Control of temperatures and flux is invaluable if proper use 
of the solar furnaces is to be made. In operation of a solar 
furnace, there is little control over the flux received by the 
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target, except by use of various mechanical methods. Among 
these methods would be: Using a control shutter placed be- 
tween the focal plane and the concentrator; shading a portion 
of the concentrator and/or heliostat; or, if possible, changing 
the focal length by adjustment and deformation of the in- 
dividual mirror segments of the concentrator. Laszlo (42) 
used a unique Venetian-blind-type shutter to control the tem- 
perature and flux. Conn (43) used a shutter in the form of a 
cylinder which moves along the axis of a paraboloidal mirror 
focus. In addition to the mechanical method of controlling 
flux, flux control can be attained by off-focus operations (20). 

With the increased use of image furnaces for research, tlie 
measurement of high temperatures and fluxes is becoming 
more important. Conn (44) points out the necessity of shut- 
ting off, by two rotating sectors, the incident and reflecti-d 
energy from the sample, so that the true temperature of the 
specimen can be measured with a disappearing-filament pv- 
rometer. Kostkowski (45), Blau (46) and Brenden et al. (17) 
have also reviewed recent developments in this field. In an 
article discussing proposed instrumentation for the large solr 
furnace planned for construction at Cloudcroft, Shank (4%) 
has discussed temperature measurement in solar furnaces. 

The development of fluxmeters for high intensity thermal 
radiation has been reported by Farber (49) and Gardon (50), 
and Hughes (51) has described the pertinent factors in the 
design of such instruments. Blau (46) also reports on two 
flux measuring devices as well as equipment for measuring 
emissivities at high temperatures. He, Gardon (50), Lohet al. 
(52) and Laszlo (42) give flux distributions across the focal 
spots of various solar furnaces. Overall flux measurements at 
the focus as a function of rim angle are given by Farber (49). 

In recent years a number of solar furnaces have been built. 
A survey of these installations in the United States, and of 
their research plans, is available in references (53 to 55). 
Reference (22) contains a table of estimated overall furnace 
efficiencies for actual furnaces. These workers, and others 
(13), have concluded that, in general, an overall furnace factor 
of 70 per cent is probably the highest efficiency that can be 
reached. 

A survey of the European solar furnaces is contained in 
(44). The world’s largest solar furnace, located at Mont 
Louis in France, is described by Trombe (31). Several fur- 
naces have been built in Japan by the Government Industrial 
Research Institute, and are described in (56). 

A directory of world activities and bibliography of significant 
literature on solar furnaces is given in (57). Development in 
the field of solar furnaces is also reported in the Journal of 
Solar Sciences and Engineering, and new furnaces are normally 
announced by The Sun at Work, a news publication of the 
Association for Applied Solar Energy. 


Are Image Furnaces 


Developments in arc image furnaces have been closely 
associated with advances in the lamp projectors and carbons 
used in the motion picture industry. As the demand for 
brighter light sources increased, the manufacturers developed 
a variety of carbons, lens, reflectors and accessories for the in- 
dustry. Today, from this wide choice of component parts, 
various types of arc image furnaces can be assembled. 

Are image furnaces are generally classified by the optics. 
A diagram of each type is shown in Fig. 3. The relative merits 
of each are discussed by Null and Lozier (27); therefore, the 
characteristics of each type will only be summarized. 

Type I—Single elliptical mirror: This is the simplest form 
of an image furnace giving a large area of low irradiance or a 
small area of high irradiance, depending on whether the are 
crater is located at the first or second focus. 

Type II—Doubdle elliptical mirror: This provides cen- 
tinuous exposures, and is one of the most practical systems, 
producing the highest irradiances with the furnace 
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aperture facing either up or down. It has shuttering ad- 
vantage over type 

Type 11I—Double parabolic mirrors: This provides essen- 
tially the same performance as type II. These have been 
widely used because of the availability of surplus Army and 
Navy searchlight mirrors. 

Type IV—Condenser-relay lenses: This has the largest 
exposure area, and, as in type I, there is no restriction on the 
size of the exposure equipment. The level of radiation is the 
lowest of all types described. 

A modification of type IIT has been proposed by Ploetz (58). 
The device consists of two identical reflecting paraboloids 
mounted face to face, so that the focal point of each falls at 
the vertex of the other. This arrangement allows both the 
source and image to be external to the optical system. 

A survey of some arc image furnaces and their uses is given 
by Davis (59) and Butler (60). Operating and performance 
data for carbon are, double elliptical mirror-image furnace and 
a comparison of other furnaces are given in (27 and 61). 

At present, the flux, temperature and size of the image spot 
is limited by characteristics of the commercial lamps and car- 
bons. For commercially available are sources and carbons, 
the maximum image diameter is about 10 mm and flux about 
225 cal/em?/see (62). It has been reported (27) that the 
Gretener Super Vent arc lamp does produce an image diameter 
10 to 50 per cent larger in proportion to the carbon diameter 
than other lamps. Fluxes as high as 350 cal/cm?/sec, equiv- 
alent to a black-body temperature of about 4000 K, utilizing 
this lamp and experimental positive carbons have been ob- 
tained. 

Experiments have been conducted both here and in Europe 
in operating arcs at very high currents (63,64). In the future, 
higher fluxes over much larger areas should be realized by the 
use of larger cored carbon. 

Experiments have also been conducted utilizing pressurized 
tungsten electrodes (65). Plans with this furnace at the 
University of California include operating in pressurized argon. 

Calculations based on plasma temperatures of high cur- 
rent argon ares indicate fluxes up to 2000 cal/em?/see may 
be obtained (60). 
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Work is also being conducted at Southwest Research In- 
stitute on a double parabolic image furnace, in which the 
source is an air blown,high current are (66) which should result 
in high irradiances. 

The life of a single carbon, commercially available, is from 
20 to 30 min when operated at the prescribed voltage and cur- 
rent. If continuous operation is required, the use of two 
lamps and a quick changeover device (similar to that used 
for motion picture theaters) can be used (27). One of the 
authors (67) has successfully used this method for the growth 
of single crystals. The use of a continuous-burning lamp 
could also accomplish the same purpose. Such carbon arc 
lamps have been described (68-70). At present, no commer- 
cial units are available. As the operating voltages and cur- 
rents are raised to gain further increases in radiation outputs, 
the burning rates of the carbon goup. Burning rates of 140 to 
160 in. per hr have been reported (27). For the normal length 
carbon this would allow operating time of about 5 min, a 
period sufficiently long for some experiments. 

Control of the flux and temperature in an arc image furnace 
can be obtained by a number of methods. For commercial 
lamps a number of types of carbons are available which can be 
operated at different voltages and currents, thereby giving 
different arc intensities. Also, the voltage and current can be 
varied over a wide range for individual carbons to give a 
range of arc intensities. Various mechanical methods, as 
described for solar furnaces, can also be used. For example, 
with the double ellipsoidal system, a rotating shutter can be 
placed at the common focus (10). Flux and temperature 
control can be achieved by defocusing either the source or the 
image receiver. Both modes of defocusing produce about 10 
per cent decrease in central brightness per millimeter of axial 
movement off-focus with the rotating-type lamp and carbons 
(27). 

Reasonable temperature control at a given temperature can 
be expected in an are image furnace. Fluctuations in the 
irradiance of an are image furnace owing to nonhomogenous 
regions in the carbon are often as small as +5 per cent. This 
is only possible with the best feed and control mechanism (60). 
Sobon (71) describes an automatic system for this purpose. 
This system has been installed on the furnace reported in 
(67). 

A comparison has been madg of the distribution of irradiance 
and illuminance across the image with three different types of 
optical systems utilizing elliptical mirrors and high intensity 
carbons (61). Fig. 4 shows this irradiance information in watts 
permm. (The flux conversion factor is 1 w per mm equals 
23.9 cal/cm?/sec.) Also reported is the spectral distribution 
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of irradiance at center of the focal spot for a double elliptical 
mirror furnace. Brown (72) measured the absolute spectral 
intensity distribution at the center of the focused image of a 
carbon are in the region of 3000-9000 A. 


Research With Solar Image F wemaces 


The results reported from solar furnace experiments have 
predominantly involved oxide systems, and the majority 
of this work has been performed at the Solar Energy Labora- 
tory at Mont Louis. Trombe (73), at the Tucson Confer- 
ence on the Use of Solar Energy—The Scientific Basis, out- 
lined some of their experiments on the fusion of refractory 
oxides and their treatment of minerals. He describes the 
large-scale fusion of alumina and sand, and experiments on the 
fritting and stabilization of zirconia. Trombe also used the 
large solar furnace to prepare citrate-soluble phosphate by the 
defluorination of apatite through fusion. Refractory minerals 
were also broken down by the intense heat, and zirconia was 
prepared from zircon, and forsterite from serpentines. More 
recent results are contained in a survey of high temperature 
research in France, reported by Trombe and Foex (7). 

Similar but less extensive work on oxides has been under 
way in Japan and is described by Mii (8). Laszlo (74) used a 
solar furnace to melt mullite and zirconia. Duwez and Loh 
used the California Institute of Technology solar furnace to 
melt specimens in their studies of the zirconia-thoria phase 
relationships (75). 

Conn (76) describes and reviews the work on the production 
of molten zirconia in a solar furnace. A number of devices for 
permitting the measurement of high temperature properties of 
materials, such as thermal conductivity, and expansion, heat 
content and emissivity are described by Glaser (77), who also 
makes brief reference to the behavior of metal samples 
heated in air in the solar furnace. 

An entirely new approach to the use of the solar image fur- 
nace in research has been made by Marcus and Wohlers (78). 
Image furnaces are really radiant energy collectors, and, in the 
case of solar concentrators, a large fraction of the collected 
energy lies in the visible and ultraviolet region where photo- 
chemical reactions occur. Such reactions are not, in general, 
very efficient, but if the photon flux can be greatly increased 
by such a concentration technique, the total conversion might 
also be increased. : 

An experimental study was therefore undertaken to deter- 
mine whether photochemical reactions would proceed by the 
same mechanisms in the focus of a solar image furnace as they 
do at lower light intensities. The first experiment in this 
new field was the study of the photoreduction of the ceric 
ion at the higher irradiances provided by the solar furnace at 
Stanford Research Institute. The collector of this concentra- 
tor faces down directly over the heliostat, so that liquids can 
be irradiated in open or closed vessels placed at the focus. 
Based on the results of many runs at different irradiances, it 
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was concluded that stoichiometry and kinetics of the photo- 
reduction of ceric ion in perchloric acid solution is the same 
at these high light intensities as observed with low light 
levels at room temperature. It was also found that, al- 
though there is local boiling at the location of the focus, the 
remainder of the liquid system, even that portion through 
which the energy is transmitted, stays cool, and therefore, it is 
possible to undertake liquid phase reactions in such a do- 
vice. 


Research With Arc Image Furnaces 


Considerably more detailed information is available in tive 
open literature on a variety of uses for arc image furnaces. A 
brief summary for each of the more pertinent is provided in 
the following paragraphs. 


Solid Propellant Ignition 


The ease of ignition of a solid propellant is dependent upon 
many factors; important ones are the amount of energy sup- 
plied to the material, and the duration of its application. In 
the classical methods of dropping the propellant in a molten 
metal or using a weighed amount of black powder to ignite 
it, the results are essentially relative, since it is impossible to 
determine the amount of heat actually transferred to the ma- 
terial. Image furnaces are known to provide measuralile 
quantities of radiant energy at very high flux rates, instan- 
taneously over a small area, with the energy pulse accurately 
controlled, and to give several advantages for studies of the 
ignition behavior of such materials. 

Stanford Research Institute has designed and utilized an 
are image furnace installation in which a specimen of solid 
propellant can be irradiated with a predetermined amount 
of energy (10). A schematic drawing of the double ellipsoidal 
mirror image furnace used for this purpose is shown in Fig. 
5. An Ashcraft high intensity motion picture projection lamp 
house supplies up to 120 cal/cm?/sec on the propellant speci- 
men inside the glass-faced exposure chamber. (This is 
equivalent to a black-body temperature of approximately 
3000 K or 5000 F.) 

High speed motion pictures taken of propellant ignition 
occurring at atmospheric pressure, show that voluminous 
quantities of particulate material are evolved from the surface 
during and after the pyrolysis brought about by exposure to 
the high flux radiant energy. The initial reaction occurs 
in the vapor phase, and the exothermic energy resulting is 
then transferred back to the surface. If the exposure is of 
sufficient length to supply a required minimum amount of 
this energy, a steady-state condition is brought about and 
sustained burning occurs. 

Ignition of a particular propellant requires a certain amount 
of energy known as the threshold energy. When less than 
this threshold energy is supplied, ignition does not occur. 
Variations of the time of exposure, with constant irradiance, 
is used to bracket the threshold energy. The variables 
affecting propellant ignitibility, such as pressure, tempera- 
ture, thermal flux and exposure time, surface condition, 
environmental gas composition, and past history can readily 
be studied. Typical results when several of these are varied 
are shown in Figs. 6 and 7. Fig. 6 shows that, at reduced 
pressures, for example, a minimum duration of energy transfer 
(about 400 millisec) is necessary for ignition to succeed, 
regardless of flux. This is related to the kinetics of the gas 
phase reaction. Fig. 7 shows that the threshold ignition 
energy is markedly affected by the presence of an external 
source of oxidizer. The authors suggest that both thus 
effect and that of pressure provide strong support for the gas 
phase theory of propellant ignition. 


Radiation Ignition of Fine Fuels 


The combustion of organic solids, in the absence of sc!f- 
contained oxidizer, involves two distinct processes: The 
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decomposition of a solid into volatile products, followed by 
the gas phase oxidation of these products. When a material 
such as cellulose is irradiated in an image furnace, the high 
temperatures produced near the surface provide volatiles at 
high rates. Ignition of this air-gas mixture will take place 
when, and only when, the gaseous fuel has a high degree of 
reactivity, and the relative concentration of the fuel and 
oxygen lie between two well-defined limits. 

The work at the Naval Radiological Defense Laboratory 
on the ignition of fine fuels by irradiation has been done with 
condenser relay, lens system, arc image furnace (79). The 
fine fuels are simulated with a series of experimental papers 
manufactured from pure alpha cellulose in two densities, 
many different thicknesses, and with varying amounts of 
carbon black. 

Three distinct types of ignition were soon recognized, each 
characterized by well-defined phenomena. 

Glowing ignition: A bright red glow appears under low 
irradiances and long exposure times throughout a constant 
equilibrium temperature rise. 

Spontaneous flaming ignition: There is a sudden appearance 
of a flame front which lasts only over the duration of the ex- 
posure produced with high irradiances and short times. It 
is independent of thickness. 

Sustained flaming ignition: Continued burning, necessarily 
preceded by spontaneous flaming, is as in the preceding case. 
This type of flaming depends upon the thickness of the mate- 
rial, in which not only must the temperature rise be rapid, 
but high equilibrium temperatures must be maintained within 
the material. 

Correlation procedures have been developed, as shown in 
Fig. 8, to successfully relate the ignition of fuels whose thick- 
ness varied by a factor of 10, densities by a factor of 2, and 
over a wide range of irradiances and exposure times. In this 
figure, the ordinate is in dimensionless units of energy, the 
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abscissa in dimensionless time units, and the individual curves 
are in thickness units. Such a correlation pattern as this on 
the fundamental processes of flash pyrolysis is essential to 
provide a basis for new methods of inhibiting ignition and 
extinguishing the resulting fire.‘ 


4 Stanford Research Institute is actively developing a modifica- 
tion of the arc image furnace for fundamental investigations in 
the area of fire research. 
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Thermal Diffusivity of Metals at High Temperatures 


If an insulated infinite slab is irradiated on one side only in 
an image furnace, it is possible to theoretically calculate the 
gross disposition of heat at any point within the slab from a 
knowledge of the thermal properties of the material (80). 
The adaptation of this concept to,an experimental technique 
requires that the area of irradiance must be approximately | 
in. in diameter, and that the flux be sufficiently high so that 
the specimens can be raised to 1000 C. The condenser relay, 
carbon are image furnace at Naval Radiological Laboratory 
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Fig. 9 Thermal diffusivity of steel measured with image 
furnace 
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Fig. 10 Specific heat of molybdenum compared with data (solid 
line) in American Institute of Physics Handbook, 1957, section 
4, p. 42 
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Fig. 11 Total hemispherical emissivity for polished platinum 
(normal emissivity values from Armour Research Institute) 
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provides such an area with a maximum irradiance of approxi- 
mately 6 cal/em?/sec. The temperature rise at two points in 
a cylindrical specimen of the metal or alloy, mounted so that 
it simulates a portion of an infinite insulated slab, is measure: 
during the time the front is irradiated. 

From the changes in temperature at the front and back of 
the specimen and heat transfer theory, the thermal diffusivit . 
data are calculated. An example of the results obtained for 
steel is shown in Fig. 9. It should be noted that the metho | 
is sensitive to phase changes, such as the sudden change i): 
slope at the Curie temperature. 


Specific Heat of Metals at High Temperatures 


A novel method has been developed at the Naval Radiolog- 
ical Defense Laboratory (81) for the determination of the 
heat capacity of metals and alloys based on the use of tl x 
are image furnace. A small, thin test specimen is suspende:| 
by its own fine thermocouple wires in an evacuated chambc’, 
the walls of which are optically black on the inside and water- 
cooled, and irradiated through a double quartz window at oue 
end of the test chamber. If the temperature of the specim«.n 
is higher than the walls, if the pressure is lower than that ::t 
which convective losses occur, and if the conduction of he:it 
through the supporting wires is negligible, all heat losses are 
by radiation. If, now, the specimen has been treated so that 
its surface is optically black, it will radiate with a Planck 
spectral distribution, with a hemispherical emissivity of 1.0. 
Since the heat loss is proportional to three quantities—the 
mass, the rate of temperature change with time and the spe- 
cific heat—this specific heat can be calculated. Results for 
molybdenum are shown in Fig. 10. 


Hemispherical Emissivity of Metals 


Hemispherical emissivity is defined as the ratio of the total 
thermal energy radiated from unit area of a solid surface, at a 
given temperature, to that radiated by a black body at the 
same temperature. It is clear that hemispherical emissivity 
will be different than the normal, but there is no simple re- 
lationship between the two. 

The experimental arrangement for the direct measurement 
of total hemispherical emissivity, as developed by Butler and 
his co-workers at the Naval Radiological Defense Laboratory, 
is to suspend a specimen of the metal in an evacuated chamber 
whose walls are optically black and are maintained at a con- 
stant temperature (82). Enclosed shutters inside the cham- 
ber can be opened to admit radiation through a quartz 
window for heating the specimen, and subsequently closed 
during the cooling, when the measurements are made. 

Typical results for the hemispherical emissivity of polished 
platinum are shown in Fig. 11. It is to be noted that the 
values are consistently higher than for a similarly polished 
specimen when measured in the customary normal manner. 
Measurements made in the same way with polished copper 
indicated a similar departure from normal values. Butler 
suggests that the differences in the two sets of data are more 
indicative of the basic difference in the definition of the meas- 
ured quantities than in any definite differences in specimen 
preparation. 


Oxidation—Erosion Behavior of Graphite 


The National Carbon Research Laboratory has conducted 
an extensive program for WADD to determine weight loss 
and decrease in length of impregnated graphite at high tem- 
peratures and moving streams of air (83). The method is to 
place the specimen of graphite, in rod or hemispherical cone 
form, at the focus of double elliptical mirror arc image furnace, 
and to irradiate the end of the material which has been care- 
fully positioned. A stream of air is directed on the end of th 
sample from a tube protruding through the center of a collect- 
ing mirror. 

With an are current of approximately 190 amp the irradi- 
ance over an area of 1 in. in diameter is 61.5 cal/cm?/sec, a= 
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measured with a calorimeter. Standard exposure times are 
30 sec. Weight and length measurements before and after 
exposure in the furnace are made for each specimen, and in 
some cases profiles of the eroded specimens recorded. 


Reflectivity and Emissivity of Graphite at 3800 K 


A method for the measurement of the reflectivity of graphite 
near the sublimation point has been developed for WADD at 
the National Carbon Research Laboratory (84). 

The radiation from a high current, carbon are furnace is 
focused by a double elliptical mirror system on the anode 
of a standard pyrometric arc in the image plane of a second 
mirror. These arcs employ spectroscopic electrodes and 
operate at a brightness temperature close to 3800 K, which is 
messured with a photocell having a spectral response cor- 
responding to that of the human eye. Two shutters are 
synchronized to separate the emitted from the reflected 
energy, as shown in Fig. 12. With both shutters open, both 
the emitted and reflected energies are recorded by the de- 
tector. With the first shutter closed, and the second open, 
only the emitted energy is recorded. The difference be- 
tween the two readings gives the reflected illuminance of the 
burning arc. Comparative measurements of reflected energy 
are made by observing the freshly coated surface of mag- 
nesium oxide. It was concluded from these studies that the 
crater surface, at the operating temperature of the arc, 
reflected 2 per cent of the incident light, and hence the 
emissivities were 98 per cent (83). 


Crystal Growing in an Are Image Furnace 


‘Two of the greatest problems in the growth of single crystals 
of refractory materials are preventing contamination and 
changes in stoichiometry at elevated temperatures. At the 
temperatures required for crystal growth these materials are 


R, A s 


A - Image Furnace Arc, D, and D, - Diaphragms 
S, and S, - Synchronous Shutters, L - Lens, C - Photocell 
E - Exposed Sample, M - Motor, R, and R, - Furnace Reflectors 


Fig. 12 Schematic diagram for reflectance measurements at 
3800 K 


usually extremely reactive, so that they cannot be contained 
in crucibles without introducing considerable contamination. 
When this problem has been overcome, for example in the 
Verneuil or flame fusion process, the atmosphere produced by 
the heat source is such that oxidation or reduction of the ma- 
terial occurs, and stoichiometry cannot be adequately con- 
trolled. 

A crystal growing facility has been developed at Stanford 
Research Institute (67) which employs the Verneuil tech- 
nique for crystal growth, but which utilizes an are image 
furnace as the heat source. Such furnaces are uniquely 
suited for this application, since the energies required to 
melt most refractory compounds can be obtained at sharply 
defined zones. Since the heat is concentrated directly on the 
crystal, no other furnace elements reach high temperatures, 
minimizing out-gassing. By employing a glass furnace 
chamber, highly purified and precisely controlled oxidizing, 
reducing or neutral atmospheres, as well as vacuum, can be 
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Fig. 13 Diagrammatic sketch of arc image crystal growing furnace 
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provided to permit close control of crystal stoichiometry at 
elevated temperatures. Contamination is virtually elimi- 
nated, since there is no crucible, and once high purity pow- 
ders are obtained, purities can be maintained, and in some 
cases, improved during crystal growth. The glass chamber 
also permits continuous observation of the crystal during 
growth. 

The are image, crystal growing furnace is illustrated 
schematically in Fig. 13. Two commercial arc image sources, 
each provided with an elliptical mirror, are positioned facing 
each other, so that their common focus can coincide in one 
plane. A water-cooled mirror in this position directs either 
of the beams vertically to another similar elliptical mirror 
which concentrates the energy downward to the crystal mount. 
The flux at the focus is approximately 60 to 80 cal/cm?/sec 
over an area } in. in diameter in the furnace chamber. 
With this system a continuous operation can be maintained 
by alternating the sources when recarboning is necessary. 

The crystals are grown on the tip of a }-in. diameter re- 
fractory rod mounted inside a Pyrex cylinder and partially 
shielded by an alumina tube. The powder of the refractory 
compound is fed through a fine orifice, just above the focal 
plane, from a tube extending through the center of the mirror 
to a feed hopper driven by a variable. amplitude vibrator. 
The desired environmental gas flowing down the tube directs 
the powder onto the molten surface. The pedestal is rotated 
and is drawn downward at the same rate as the crystal grows. 


High Temperature Microscopy 


An entirely new concept in high temperature microscopy 
has been reported by Newkirk and Brenden (85). An are 
image furnace was constructed for use with a microscope 
having a working distance of 7.25 in. and a numerical aperture 
of 0.1. A drawing of the equipment is shown in Fig. 14. 

Maximum irradiance of the specimen while under observa- 
tion is 90 cal/em?/sec. The interaction of uranium dioxide 
with silicon or silicon carbide was photographed at 2200 C 
under a magnification of 22 X with a motion picture camera. 
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“OPTICAL RELAY SYSTEM 


Fig. 14 Arc image furnace for high temperature microscopy 


Plus X or Kodachrome-type film were used with exposures 
made at the rate of 24 frames per sec, each lasting 45 sec. 
Analysis of the motion picture film can provide rate inform:- 
tion on the stability of solids in various atmospheres and in 
contact with other reactive solids. 


Ablation 


The success of reinforced plastics as heat shields for short- 
time hypotherma] environments has resulted in a large number 
of studies on the mechanism of ablation. Mixer and Mary- 
nowski (86) concerned themselves with the chemical reac - 
tions involved. They viewed the total hypothermal effe:t 
as the resultant of three components—a velocity effect, 1 
thermal eifect and a plasma effect, and used a plasma jet for 
providing all three of these. The separate effects of the se: - 
ond and third components were simulated by exposure in t] e 
are image furnace and an RF discharge, respectively. They 
used the same furnace as developed for the propellant ignition 
studies previously described, with one modification to perm t 
a sampling of the gaseous decomposition products evolved :n 
the exposure chamber. Mass spectrometric analyses of 
these products were helpful in deducing some of the chemistry 
involved in ablation. 
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Astronomical Constants and Their 
Importance in Lunar Trajectory 


Determination 


TROSS' 


Aeronutronic Div., 
Ford Motor Co. 
Newport Beach, Calif. 


Whenever the problem of calculating orbits for space vehicles arises, the accuracy of such calcula- 


tions receives close scrutiny. No matter what technique is adopted for the orbit calculation, uncer- 
tainties enter into the calculations which cannot be eliminated by numerical techniques and mathe- 
matical representations. These uncertainties result from insufficient accuracy in the constants 
utilized for the representation of such effects as Earth’s bulge and equatorial radius, the lunar and 
solar masses and their position coordinates. In this paper an attempt is made to appraise the 
effects astronomical constants have on lunar vehicle impacts. This study shows that lunar impacts 
can only be predicted to fall into an area of about 78 miles in radius. More accurate predictions can 
only be made when the terrestrial and astronomical constants are improved. 


LTHOUGH concentrated efforts are being directed to im- 
prove the geocentric and selenocentric constants at the 
present time, only limited contributions can be made until 
such time as the existing and future satellite data can be 
analyzed completely. Until the time when these constants 
are improved, the present values with their respective uncer- 
tainties must remain the basis for orbit determination activi- 
ties. It is very important, however, to appraise the signifi- 
cance of each constant and the effect of all known uncertainties 
on some quantity, such as the impact point of a ballistic tra- 
jectory. The results of such an investigation on lunar trajec- 
tories are presented in this report, and the constants to which 
lunar trajectories are most sensitive are identified. 


Perturbation Effects and Their Importance in 
Precision Lunar Trajectories 


What perturbations must be included in the equations of 
motion for accurate lunar trajectory calculations? Everyone 
who has ever faced the problem of writing equations of motion 
has had to answer this question. In this section some of the 
considerations which should enter into such an analysis are 
discussed. 

A convenient form in which the equations of motion may be 
considered is to add accelerations resulting from perturba- 
tions to the two-body equation of motion for a lunar vehicle, 
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whose mass is negligible, moving under the influence of Earth’s 
gravitational field only. Thus the equation of motion takes 
the form (3 and 7)? 


fo. = + P [1] 
where 
r;.; = vectorial distance from body ¢ to body j 
m;: = mass of body 7 
P = those perturbations not yet identified 


Numerical subscripts used in this paper to identify bodies are: 
1, vehicle; 2, Earth; 3, moon; 4, sun; 5,6, ..., solar bodies. 

A lunar vehicle experiences accelerations caused by the 
masses of the moon, sun and other solar bodies. The per- 
turbative accelerations resulting from these effects take the 
form 


Ti,2 Tia 
Maximum values for the differential accelerations are given in 
Table 1. From these tabulated values it is noted that per- 
turbations resulting from planets other than Earth are quite 
small and may therefore be ignored even in a high precision 
lunar trajectory determination procedure. 

The Earth is an oblate spheroid, a fact not taken into ac- 
count in the two-body equations of motion. Therefore, 
another perturbation owing to Earth’s bulge must be con- 


2 Numbers in parentheses indicate References at end of paper. 
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Pu = + KBr) [3] 


where J and K are related to the DeSitter constant and are 
the second and fourth harmonic terms in Earth’s gravitational 


potential. Their respective coefficients may be given as 
A(x) = A(y) = A(z) — 2 = 1 — 5U? 
B(x) = By) = B(z) — 12 + 28U? 


3 — 42U? + 


where U is a function of z. 

The rectangular geocentric inertial coordinates are defined 
so that the z-y plane identifies the Equator plane with zx 
directed toward the vernal equinox, z directed toward the 
north pole, and y completing the right-hand system. Ac- 
celerations resulting from the terrestrial harmonic terms as a 
function of r have been tabulated in Table 2. From this 
tabie it is apparent that the bulge harmonics of the second 
order are quite important for the entire lunar trajectory cal- 
culation, but the fourth-order harmonics become negligible 
after the vehicle is more than 8 g-radii away from the geo- 
center, so that these terms are not needed beyond that dis- 
tance. 

Thus far, solar system perturbations have been analyzed 
and identified. It is of interest to determine the effect these 
perturbations have on the lunar impact point. To analyze 
this effect a standard trajectory® is calculated in which the 
nominal values of all constants are employed. Subsequent 
trajectories called “altered trajectories’ are considered in 
which all parameters are the same as those used in the 
standard trajectory, only the perturbative term under con- 
sideration is set to zero. The impact point of the standard 
trajectory is then compared with the impact point of the 
altered trajectory and displacement on the lunar surface is 
determined. See Table 3. 

It will also be of interest to determine maximum relative 
accelerations imposed on the vehicle by various perturbations 
during its transit to the moon. For this purpose Table 4 has 
been prepared. 

Thus far the perturbations caused by the terrestrial bulge, 
the sun, solar system planets and the moon have been con- 
sidered. The moon’s effect upon lunar vehicles is, of course, 
not only that of a point mass, but its figure may well be a con- 
tributor of additional perturbations. A number of authors, 
such as Jeffreys (10), think that the moon’s figure may be a 
triaxial ellipsoid whose longest axis is directed toward Earth. 
Based on this assumption, estimates for the lengths of the 
longest, the intermediate and the shortest axes, represented 
by a, b and c respectively, can be made (2) 


a = 1738.57 + 0.07 km 
b = 1738.21 + 0.07 km 
c = 1737.49 + 0.07 km 


The moments of inertia about each axis are symbolized by A 
Band C, respectively, and experimentally determined values 
for their ratios are (10) 


(C — A)/C = 0.6269X10-* + 2.7 X 10% 
(B — A)/C = 0.2081X10-* + 9 X 1077 


The acceleration owing to the moon’s triaxiality is then 


° A standard trajectory is one for which the velocity vector at 
impact intersects the lunar surface at normal incidence. 

‘In the event the altered trajectory does not impact on the 
moon’s surface the point of nearest approach is used in place of 
the impact point. 
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sidered. Bulge perturbations are obtained through the rela- 


Table 1 Solar system perturbations on lunar vehicles 
Perturbative Acceleration on vehicle 
(Earth g’s) geocenter vehicle 
distance- 
Object 32 g-radii 64 g-radii 
Sun 16.5 XxX 107 33.0 xX 107 
Moon 10-4 0.16 
Venus 0.15 x 10-° 0.30 x 10-9 
Mars 0.25 10- 0.53 x 
Jupiter 0.21 10-” 0.43 « 10-% 


Table 2 Accelerations resulting from terrestrial 


harmonics 

1 0.6494 107? 0.7232 10- 
4 0.4058 x 1078 0.2824 x 10-6 
8 0.1015 10-3 0.1766 1077 
64 0.1585 « 10-5 0.4310 « 107" 


(U2? = sin? 6 = 1) 


Table 3 Perturbational effects on lunar impact 


Deviation from standard 
trajectory impact, 


Constant nautical miles 
J 368.0 
K 9.83 
ms 374.0 


J, K, ms, ms 5500.0 


(i.e., two-body orbit) 


Table 4 Maximum accelerations resulting from various 


perturbations 
Acceleration relative to geccenter, 
Perturbing —Earth g’s-— 
object At Earth 32 g-radii At moon 
(64 g-radii) 
Earth 1.0 O82 x 
Second terres- 
trial har- 
monic 0.6 10-2? 0.6 xX 107 04 X 10-8 
Fourth terres- 
trial har- 
monic 07 xX 10+ O7 <107* 
Moon x10 O86 
Sun 0.16 X 0.33 x 
Venus 0.15 x 10-* 0.30 x 10-° 
Mars 
Centrifugal 
force owing 
to precession 0.3 xX 10-" 
written 
—— K(131)C X 
J. 
U,(r3,1) AU, (r3,1) 
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certainty +0.0154 sec is added. 


Table 5 Impact deviations owing to uncertainty in astronomical constants 


from standard impact, 


Constant Value used (1) nautical miles Remarks 

J (1623.41 + 4) x 10-* 8.5 

m (332488 + 43)m2 0 

ms (122888 + 39)10~7me 0.016 

™s Ephemeris value* + 0.0154 sec 4.6 

a, 6.378270 (+ 11 XK 10-*) megameters 47.41 from velocity uncertainty 
32.67 from altitude uncertainty 
78.32 from altitude and velocity 


* The apparent value from the “American Nautica Almanac and Lunar Ephemeris” (8) is used in the calculation, and the un- 


Maximum deviation 


uncertainty 


where r;,; is the vector from the selenocenter to the vehicle and 
K(r3,4), Ui(t3,1) and U2(r3,) are functions of the direction co- 
sines of f3;. From these values of the ratio of the moments 
of inertia an upper bound for the bulge perturbation can be 
calculated. At the moon’s surface this bulge perturbation is 
1.3 X 10-* Earth g’s; halfway between Earth and the moon 
it is about 5 X 10-" Earth g’s, and at Earth’s surface it is 
about 1.2 X 10~-" Earth g’s. It can therefore be concluded 
that the lunar bulge should be considered as a perturbation in 
computing orbits for lunar satellites, but for Earth-moon tra- 
jectories this perturbation may be ignored. 

Since the equations of motion may be written to apply to a 
nonrotating system of axes, the question of the axis rotating 
slowly owing to precession and nutation is next considered. 
The value of the precession is about 50 sec per year. The cen- 
trifugal acceleration at the distance of the moon (64 g-radii 
approximately) is then about 1.4 X 10~" g-radii/min? (Table 
4). 

Planetary and solar position coordinates are referred to the 
mean equinox of 1950.0. Thus these coordinates are referred 
to a rotating coordinate system; its rate is about 50 sec of arc 
per year. The precession displacement for the sun will be 2.4 
x 10-* times the perturbation terms of the sun, which is 
small. The nutational terms are even of a lesser order of mag- 
nitude and can be neglected. The largest of the short-period 
nutational terms gives rise to a maximum acceleration of 2 X 
10— Earth g’s. Thus, it can be seen that sufficient justifica- 
tion exists for treating a system, which in reality is rotating, 
as inertial, for the practical purposes of calculation. 


Constants 


For the purpose of this study, every effort has been made to 
utilize the latest and most reliable constants possible. The 
values of the astronomical constants employed in the lunar 
vehicle orbit determinations discussed in this report are: 

Sun’s mass (1) 


m4 = 332488(1 + 0.00013) m. = 332488 + 43 m2 
Moon’s mass (1) 

m; = [1/81.375(1 + 0.0003)] m. = (122888 + 39)10~7 mp 
Second harmonic in Earth’s gravitational potential (9) 
J = (1623.41 + 4)10-* 

Fourth harmonic in Earth’s gravitational potential (9) 
K = (6.37 + 0.23)10-* 
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Geocentric gravitational constant (1) 
k. = 0.07436574 g-radii’/?/min 
k. min = 13.447052 min 
Earth’s equatorial radius® (1,6) 


a. = (6.378270) (+10 X 10~*) megameters 


‘ Earth’s mass is chosen as unity. 


Uncertainty in the Constants 


A procedure similar to that outlined previously in determin- 
ing the effect perturbations have on the lunar impact point 
is employed to determine the effect uncertainties in the astro- 
nomical constants have on lunar trajectories. 

Variations in lunar trajectory impact sites resulting from 
the uncertainties in the constants are determined by first adopt- 
ing a standard reference trajectory determined on the basis of 
the nominal value of the constant. Deviations from the 
reference trajectory impact are then determined by altering 
the constant under consideration to its limiting value within 
the recorded uncertainty; trajectories so determined are 
called “altered trajectories.”” Thus, the trajectory variations 
owing to the uncertainties in m3, m4, J’ and 7,4 are determined.® 
The uncertainty in Earth’s equatorial radius a, is reflected 
in the gravitational constant and in the conversion of burnout 
velocity and altitude from laboratory units into astronomical 
units. Earth’s radius is given as a, = 6.37827 (+11 X 10~*) 
megameters = 20926040.8229 + 230.1865 ft. With this 
value it is found that 1 g-radius/k,—! = 25936.342 + 0.2853 
fps. Suppose that a burnout velocity of § = 34704.2391 fps 
is given; then, upon converting this velocity it is found that 


1.338069202 g-radii/k,-! > > 1.338039796 g-radii/k.— 


Thus the uncertainty in Earth’s equatorial radius may result 
in an error in conversion of the burnout velocity units, thereby 
causing the vehicle to follow an altered trajectory whose impact 
differs from the design impact point. The maximum impact 
deviation resulting from a velocity uncertainty is about 47 
nautical miles. See Table 5. 

In converting the burnout altitude from laboratory units to 


5 At the time these calculations were performed, the latest value 
of a, which is a, = 6.378145 (1 + 11 X 10-6) megameters had not 
been determined. Consequently the calculations discussed in 
this paper are based on the older value. 

6 x4 is the moon’s horizontal parallax. 
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Az Ay Ax Ar 
(1074) (1073) (1073) (1073) 
+660 +150 —-330 +330 
+600 +135 ~-300 +300 
+540 +120 —-270 +270 
+480 +105 —-240 +240 
+420 +90 +210 
+360 +759 -180 «+180 
+300 +60 —-150 +150 ‘ 
+240 +45  -120 +120 
+180 +30 90 + 90- 
+120 +15 + 60- 
+ 60 0 -30 + 30 
° 5 


t x 103 (min) 


Fig. 1 Maximum trajectory deviations owing to uncertainties in Earth’s gravitational radius 


astronomical units, the uncertainty in Earth’s equatorial 
radius is again reflected. Suppose the burnout altitude is 
given as h = 2092604.08229 ft, then upon converting this 
altitude into astronomical units it is found that 


1.1000121 g-radii > h > 1.0999879 g-radii 


The nominal value is h = 1.1 g-radii. The maximum impact 
deviation owing to an altitude uncertainty is about 33 nautical 
miles. See Table 5. Combined, the maximum deviation re- 
sulting from both altitude and velocity uncertainty is about 
78 miles, in contrast to other uncertainties, such as the 
second harmonic term in Earth’s potential, the sun’s mass, the 
moon’s mass and the moon’s horizontal parallax, all causing im- 
pact variations which never exceed 10 nautical miles. It is, 
therefore, evident that the uncertainty in Earth’s gravitational 
radius may cause the most significant impact variations. 

A summary of the maximum impact point variations owing 
to the uncertainty in the astronomical constants is given in 
Table 5. 


Orbit Deviations 


!t has been shown that the uncertainty in the astronomical 
constants causes an uncertainty in the orbit the vehicle fol- 
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lows from Earth to the moon. To illustrate how a maximum 
deviated orbit may vary from the nominal orbit, differences 
between two such 3}-day trajectories were computed, and the 
results are shown graphically in Fig. 1. 
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Heat Diffusion in Gases, Including 
Effects of Chemical Reaction 


C. FREDERICK HANSEN' 
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Moffett Field, Calif. 


The diffusion of heat through gases is treated where the coefficients of thermal conductivity and 
diffusivity are functions of temperature. The diffusivity is taken proportional to the integral of 
thermal conductivity, where the gas is ideal, and is considered constant over the temperature inter- 
val in which a chemical reaction occurs. The heat diffusion equation is then solved numerically for 
a semi-infinite gas medium with constant initial and boundary conditions. These solutions are 
in a dimensionless form applicable to gases in general, and they are used, along with measured 
shock velocity and heat flux through a shock reflecting surface, to evaluate the integral of thermal 
conductivity for air up to 5000 K. This izitegral has the properties of a heat flux potential and re- 
places temperature as the dependent variable for problems of heat diffusion in media with variable 
coefficients. Examples are given in which the heat flux at the stagnation region of blunt hypersonic 


bodies is expressed in terms of this potential. 


HIS paper will consider some aspects of heat conduction 

through gases. Specifically, it will be confined to the case 
of pure diffusion; that is, heat transfer by mass convection 
and by radiation will not be treated. Finally, it will consider 
briefly the role of the diffusion process in heat transfer to the 
stagnation region of high speed vehicles. 


Heat Diffusion in Ideal Gases 


The classical solutions to the heat diffusion equation have 
generally made use of the assumption that the coefficient of 
thermal conductivity is a constant [see (1 and 2) for example ].? 
However, for the problem of interest here, the conduction of 
heat through gases, thermal conductivity must be treated as 
a strong function of temperature (3). In such a case, it is 
useful to replace the usual dependent variable, temperature, 
with a quantity ¢ 


kdT (1] 


which is the coefficient of thermal conductivity integrated over 
temperature. Then the differential equation of heat conduc- 
tion takes the usual form of the diffusion equation 


d¢/dt — = 0 [2] 


where a is the thermal diffusivity. The dependent variable ¢ 
will be called the heat flux potential, since the heat flux at any 
point in the medium is just the gradient of ¢ 


q = grad ¢ [3] 
Presented at the ARS Semi-Annual Meeting, June 8-11, 1959, 
San Diego, Calif. 
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Note that just as the reference level used for temperature is 
arbitrary, so the lower limit of the integral (Eq. [1]) is arbi- 
trary. However, it will prove convenient to choose this lower 
limit somewhere near absolute zero for the case of conduction 
in gases. 

Consider now the one-dimensional heat flow through a semi- 
infinite gas medium, where the initial conditions and the 
boundary conditions are constants. The initial potential 
throughout the gas will be designated by ¢,, and the boundary 
value by go, the diffusivities by a,, and ap, respectively (see 
Fig. 1). Normalized coordinates will be used such that go and 
a are both taken to be unity. The results can easily be gen- 
eralized to account for arbitrary units of these boundary con- 
ditions, of course. 

The idealization of constant initial and boundary conditions 
will be approximately realized in a physical situation where 
the reservoir is a slab of material with a very large heat 
capacity and a large thermal conductivity, such as a metal. 
Then the wall can soak up heat fast enough to maintain nearly 
constant conditions at the interface. The slab might be sud- 
denly immersed in a constant temperature gas, for example, or, 
alternatively, the gas might be heated suddenly by a plane 
shock wave reflecting from the solid interface. 

Boltzmann (4) has shown that the solution to this problem 
may be expressed with perfect generality as a function of 
a/Vt. Therefore, the time and distance variables are trans- 
formed to the single dimensionless parameter y 


y = 2/V4act [4] 


where d@ is the diffusivity of the gas at the boundary. The co- 
efficient ad) is unity, of course, in the present normalized co- 
ordinate system, but this is not an essential feature of the 
transformation. In terms of this parameter y, the partial 
differential equation of heat conduction becomes the dimen- 
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sionless, total differential equation 


= 6 [5] 


The factor 4 in the parameter y (Eq. [4]) is an arbitrary 
stretching factor. It is chosen merely so that the solutions to 
Equation [5] reduce to the usual error function form when 
the diffusivity is constant, that is when a equals a. The 
merit in using the normalized form for the potential ¢ (that 
is, v/¢o) is that one set of integrations for Equation [5] will 
suffice for all possible boundary conditions. All that is re- 
quired before one proceeds with the integration is that the 
diffusivity a be evaluated. Two cases will be considered, one 
in which the gas is ideal and inert, and the second in which the 
gas is in local equilibrium but is chemically active with a large 
heat of reaction. 

According to the kinetic theory of inert gases (5), the co- 
efficient of thermal conductivity is approximately proportional 
to a power n of the temperature, where n is close to 3. Thus 
the integral of thermal conductivity is proportional to the 
n + 1 power of temperature 


S kdT To [6] 


Now the thermal diffusivity equals the conductivity divided 
by the heat capacity per unit volume Cp. Since the density 
p is inversely proportional to temperature, for an ideal gas, 
the diffusivity is proportional to the n + 1 power of tempera- 
ture also 


k/Cp « 7] 


it follows that diffusivity is proportional to the potential ¢, 
and in the present normalized coordinates this means that 


a=¢ [8] 


Sometimes a linear relation between a and ¢ with a finite 
intercept will best fit measured values, but this can always 
be transformed to the direct proportionality of Equation [8] 
by appropriate adjustment of the lower limit of the integral, 
Equation [1]. Generally, this limit will be close to absolute 
zero. 

It may be remarked that even if one accounts for the 
variations of heat capacity which occur in a real gas, the pre- 
ceding relations are valid to the order of approximation that 
Prandtl number is a constant. This can be seen from the fact 
that diffusivity is just the kinematic viscosity divided by the 
Prandtl number 


_ u/p 
Pr [9] 


and the kinematic viscosity is again approximately propor- 
tional to the n + 1 power of temperature, at least in gases 
composed of neutral particles. 

When the result of Equation [8] is used, the heat diffusion 
equation takes the simple but nonlinear form 


dy 


ay? + 2 0 [10] 


Analytie solutions to Equation [10] are known, but none 
which satisfy the boundary conditions, namely 


= 1 [Ila] 
lim 9 > ¢ [11b] 
yoo 


However, it is relatively easy to integrate the equation 
numerically, starting from a given value of the boundary 
derivative go’ = (dg/dy)o, and terminating as the solution 
asymptotically approaches a limit ¢.. The solutions are 
something like error functions stretched slightly out of shape. 
Three of these solutions are shown in Fig. 2 as functions of 
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y/ VJ Yo. The solutions ¢ are shown divided by the error 
functions, 1 + (¢.. — 1) erf (y/ V. en). The ¢ functions rise 
more steeply than the error functions near the origin, go 
through a maximum deviation, then approach the error func- 
tion as a limit. 

The value of ¢,, is uniquely related to the derivative at the 
origin, go’. This relation is shown in Fig. 3. The derivative 
¢o’ is equivalent to a dimensionless heat flux at the boundary, 
whereas the normalizing function 1 + / 1/4 go’ is the value 
which ¢,, would have if diffusivity were a constant. We shall 
return to this relation between ¢,, and ¢’ after considering the 
case of the chemically reacting gas. 
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Heat Diffusion in a Chemically Reacting Gas 


A chemical reaction in the gas behaves as a reservoir which 
soaks up heat as temperature is increased, and liberates heat 
when temperature drops. Consequently, the specific heat is 
very large if the heat of reaction is large compared to RT. 
Hirschfelder (6) has shown that the coefficient of thermal con- 
ductivity, for gases in local equilibrium, is also very large; 
in fact it is approximately proportional to the specific heat. 
Thus, the Prandtl number is relatively constant and the dif- 
fusivity is a function of temperature which is not greatly 
affected by the chemical reaction. The integral of thermal 
conductivity, on the other hand, is greatly increased as a re- 
sult of the reaction. This situation is illustrated in Fig. 4. 

The diffusivity increases linearly with g up to the point 
where the chemical reaction occurs; there it flattens out until 
the reaction is about complete, and then it increases again in a 
more or less linear manner. The solid lines are theoretical 
estimates, calculated from results in (7) for the case of air in 
which oxygen dissociation occurs. The pressure dependence 
of the curve for ¢ less than 2¢, has been removed by normaliz- 
ing both ordinate and abscissa with the factor 


ge = [12] 


= boundary temperature 

T. = temperature where thermal conductivity is a maxi- 
mum, i.e., where diffusivity is most nearly inde- 
pendent of the conductivity integral 


The temperature 7. may be calculated from equations 
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Fig. 5 Chemical equilibrium heat diffusion 


derived by Butler and Brokaw (8) for the thermal conduc- 
tivity of reacting gas mixtures. 

It should be noted that the calculation of the relation be- 
tween diffusivity and the conductivity integral does not do- 
pend on precise numerical values of the transport coefficients. 
For example, the estimate shown for air in Fig. 4 is based on 
simple kinetic theory and very approximate collision cross 
sections (7). However, the corrections introduced by more 
exact calculations affect both the diffusivity and the conduc- 
tivity simultaneously in such a way that the functional rela- 
tion between the two is maintained. 

The relation for reacting gases in general will be similar ‘o 
that for air. The dashed curve on Fig. 4 is a limit which is 
approached as the heat of reaction becomes very large. T) is 
limiting relation has been used in the integrations whi:h 
follow; that is, a is assumed to equal ¢ up to the point v., 
and thereafter is taken to be constant. Up to the point », 
then, the solutions are the stretched out error functions which 
were discussed earlier, and which are shown in Fig. 2. At +, 
these solutions are joined by the solution for constant diffusi v- 
ity which has a matching slope at the junction (d¢/dy). 


re. (de y? 


(ex J) [13] 
Ve Ve 


In Fig. 5, the limit of the potential ¢., is shown for a reacting 
gas as a function of the boundary heat flux go’. The solution 
now depends on the value of ¢., which in turn depends on the 
boundary temperature, the pressure and the chemical reaction 
being considered. For example, in the case of air and a 
boundary temperature of zero deg C, the values 20, 40 and 60 
for ¢- correspond to pressures about 10~‘, 10° and 10? atm, 
respectively. The solutions for values of ¢. corresponding to 
other pressures and different boundary temperatures can be 
obtained fairly accurately by interpolation between the nu- 
merically integrated solutions, such as the curves plotted in 
Fig. 5. The inert gas solution (Fig. 3) corresponds to ¢, equal 
to infinity (Fig. 5). It should be noted that these solutions 
can only be used up to the point where the chemical reaction 
goes to completion. Beyond this, the increase in diffusivity 
must again be taken into account. The domain of validity is 
a characteristic of each specific reaction. Within this limita- 
tion, the curves of Fig. 5 apply to chemically reacting gases 
in general. 


Heat Flux Potentials for Air 


The foregoing solutions (Fig. 5) have been put to use, in 
conjunction with experiment, to evaluate the heat flux po- 
tential for air as a function of temperature (3). In this ex- 
periment, a plane shock wave is reflected from the end wall of a 
shock tube. The heat flux to the wall is measured, and this 
fixes the abscissa for the graphical solution in Fig. 5. The 
temperature at the wall and the pressure of the gas after the 
shock reflection determine the value of ¢, (see Eq. [12]). 
Hence the value of the ordinate in Fig. 5 is fixed, and the value 
of ¢.. is determined. The temperature of the gas associated 
with this heat flux potential is taken to be the equilibrium 
temperature after the shock reflection, and this is a known 
function of the measured shock velocity. This procedure, of 
course, implies that the experiments are conducted under 
conditions where the chemical relaxation times are short com- 
pared to the test interval. 

Results of the experiments and calculations are shown in 
Fig. 6. The heat flux potential is shown divided by the inert 
gas value, where the coefficient n is taken to be } 


2 T cal 
inert = = = 3. /2 4} 
Pinert = 3 koTo (7) 3.2 X 10-*7 [1 
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and 7 isin deg K. Thus the ratio ¢/¢inert Should be unity if 
the coefficient of thermal conductivity is proportional to the 
one half power of temperature. This appears to be approxi- 
mately true up to the temperature where oxygen dissociation 
begins. At the pressures involved in the present shock tube 
experiments, this occurred at about 2500 K. At higher 
temperatures, the heat flux potential becomes more than 
twice as large as the inert gas value as a result of the 
dissociation of oxygen. The solid curve on Fig. 6 shows 
a theoretical estimate of the heat flux potential, based 
on the coefficients of thermal conductivity presented in 
(7) The experimental points were obtained by the 
method described in (3), except that somewhat better 
instrumentation has been used to improve the reliability 
of the data. The small departure from ideal gas values 
which occurs at low temperatures is a result of the effect 
on thermal conductivity of vibrational energy excita- 
tion and of the formation of nitric oxide. However, the oxy- 
gel dissociation process has the large heat of reaction re- 
quired to produce a sizable increase in the heat flux potential. 
The theory and experiment clearly show the same behavior 
owing to this effect. 

lt will be of interest next to examine an application of the 
heat flux potential to some problems of heat transfer in hyper- 
sonic flow. 


Applications to Hypersonic Flow 


Normally, heat transfer in fluid flow exhibits a complex in- 
terdependence between conduction and convection processes. 
However, at the stagnation region of hypersonic vehicles, it 
has been shown (9) that, to a first approximation, the conduc- 
tion and convection effects are separable, and solutions can be 
expressed in a form which has a simple physical interpretation 
in terms of the heat flux potentials. The heat flux at the wall 
is given by the difference between the potential at stagnation 
conditions g; and its value at the wall ¢. all divided by a 
characteristic length 6 


q = (¢t — ¢w)/6 [15] 


The quantity 6 can be interpreted as a thermal boundary layer 
thickness defined such that the gas temperature at the edge 
of the layer is a maximum, approximately the stagnation 
temperature (Fig. 7). According to this model then, the heat 
transfer at the stagnation region occurs essentially by a diffu- 
sion process, although the mass convection processes are re- 
sponsible for establishing the magnitude of 6. For an axially 
symmetric stagnation region, for example, the thickness 6 
is approximately given by [see (9) ] 


6 r pi Be 


where the subscript ¢ refers to stagnation conditions and the 

subscripts 1 and 2 refer, respectively, to the free stream and 

to conditions just following the shock wave (see Fig. 7). The 

Reynolds number Re is based on free stream conditions and 

= a length equal to the effective radius of curvature of the 
ody. 

None of the factors in the expression for the thermal 
boundary layer thickness (Eq. [16]) is strongly influenced 
by chemical reaction. In the integral ¢, however, the inte- 
grand k is greatly increased, and the upper limit of the integral 
T. is greatly decreased by the reaction 


kdT [17] 


It can be shown that these two effects almost compensate one 
another, and to a first approximation ¢, is just a function of 
the vehicle velocity, the Prandtl number and a weighted 
average coefficient of viscosity 

u?f/2Pr [18] 
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where 
= Tt Tt 
a= f, /f, [19] 


The viscosity coefficient, however, is not strongly influenced 
by chemical reactions, unless a high percentage ionization is 
produced [see (7)]. Consequently, stagnation region heating 
is relatively insensitive to dissociation effects, and is not a 
good quantity to measure for the purpose of evaluating trans- 
port properties or real gas effects in dissociating gases. It 
also becomes clear why calculations based on ideal gas proper- 
ties have been observed to have the correct magnitude in 
spite of the unrealistically high stagnation temperature pre- 
dicted by this method. Another method, found to give 
reasonable results where chemical reactions occur, assumes 
heat flux proportional to the enthalpy gradients. This is, of 
course, consistent with the fact that thermal conductivity is 
approximately proportional to the specific heat in the domain 
of the reaction, so that the integral ¢ is approximately propor- 
tional to enthalpy. 

The heat transfer to cylindrical stagnation regions can like- 
wise be expressed in terms of the heat flux potentials. Equa- 
tion [15] applies to this case without change, except that 6 is 
slightly larger than for the axially symmetric case [see (9) ]. 


1 V/ Re (2 1)" (2)" (20) 
6 r 4 Pi Me 
The analysis can also be extended to yawed cylindrical stag- 
nation regions, such as the leading edge of a swept hypersonic 
wing. It is found that the chemical reaction effects are can- 
celed in somewhat the same manner as in the stagnation 
heating of axially symmetric bodies. In the limit as the 
stagnation temperature becomes large compared to the wall 
temperature, the ratio of heat flux at angle of yaw q(A) to the 
value at zero yaw q(0) is approximately given by 


= cos A (cos A+ (3) Pr sin? [21] 
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For a Prandtl number of 3 and yaw angles up to 70 deg, this 
is very close to the cos*/? ) relation proposed by Goodwin, 
Creager and Winkler (10) for the total heat transfer to yawed 
cylinders. This effect has been checked experimentally in 
hypersonic flow (9,11). 


Concluding Remarks 


Although the preceding solutions are based on some de- 
cidedly first-order approximations to the Navier-Stokes equa- 
tions, it is felt that the essential functional relationships which 
influence stagnation region heat transfer are probably retained, 
inasmuch as the results compare favorably with more rigorous 
calculations which are now available, such as those per- 
formed by Fay and Riddell (12). Perhaps the principal value 
of these approximations, at least in domains where the rigorous 
solutions have been worked out, lies in providing a somewhat 
simplified physical insight into the problems. It may be 
noted, however, that the concept of the heat flux potential is 
free of restrictions concerning the behavior of the Lewis num- 
ber, and it can be applied to problems where multiple diffusion 
coefficients are involved in the heat transfer process. This 
situation exists, for example, where ionization or two or more 
simultaneous chemical reactions occur. However, at ioniza- 
tion temperatures the effects of emission and adsorption of 
radiation in the gas may be difficult to analyze separately from 
heat diffusion. In this case, temperature again becomes the 
normal dependent variable to use in the problem, and the 
functional dependence of thermal conductivity on tempera- 
ture becomes important. This function could be approxi- 
mately obtained from values of the integral ¢, but the dif- 
ferentiation of the data involves a considerable loss of ac- 
curacy. It is advisable to treat a favorable comparison be- 
tween theoretical and experimental values of ¢ as satisfying 
merely a necessary but not sufficient check of the theory. 

In summary, some solutions have been presented for the 


nonsteady diffusion of heat through a semi-infinite gas having 
constant initial and boundary conditions. It is found that the 
heat flux potential function ¢ is a natural parameter to use xs 
the dependent variable for problems involving heat diffusion 
through media in which thermal conductivity is a function of 
temperature. This potential has been evaluated by experi- 
ment for air up to 5000 K, including the influence of oxygen 
dissociation, and the results agree reasonably well with theo- 
retical estimates. Finally, some relations between the heat 
flux potentials and heat transfer at the stagnation region of 
high speed vehicles have been discussed. 
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Gradient Theory of Optimal 
Flight Paths 


HENRY J. KELLEY! 


Grumman Aircraft Engineering Corp. 
Bethpage, N. Y. 


An analytical development of flight performance optimization according to the method of gradi- 
ents or “‘method of steepest descent” is presented. Construction of a minimizing sequence of flight 
paths by a stepwise process of descent along the local gradient direction is described as a computa- 
tional scheme. Numerical application of the technique is illustrated in a simple example of orbital 
transfer via solar sail propulsion. Successive approximations to minimum time planar flight paths 
from Earth’s orbit to the orbit of Mars are presented for cases corresponding to free and fixed bound- 
ary conditions on terminal velocity components. 


HE PAST decade has seen considerable progress in 

techniques for the determination of flight paths which are 
optimal in the sense of various performance criteria. Treat- 
ments have employed almost exclusively the classical ‘“in- 
direct’’ method of the calculus of variations which is based 
on the reduction of variational problems to differential equa- 
tions. A number of works on this subject are listed in (1 
through 9).2. An excellent bibliography is presented in the 
survey paper of (10). 

Although many interesting results have been forthcoming 
from analytical solutions of the Euler-Lagrange differential 
equations governing optimal flight, the idealizing assump- 
tions usually invoked limit their applicability in practical 
situations. Under more realistic assumptions, a numerical 
attack on these equations is required, and in this approach a 
serious difficulty may arise in the satisfaction of two-point 
boundary conditions. [See, for example, (11, 12 and 13).] 
This difficulty becomes a limiting factor where the order of 
the differential equations governing the basic system is four 
or higher. 

Attention is directed in the present work to one of the direct 
methods of the calculus of variations, namely the method of 
gradients or “‘method of steepest descent,”’ which offers cir- 
cumvention of the two-point boundary value difficulty. The 
method also possesses the attractive feature of simultaneous 
optimization with respect to configuration parameters. 

The notion of descent along the gradient direction was 
originally introduced by Hadamard in connection with 
mathematical existence proofs (14). Only in recent years has 
it found practical application to multivariable minimum 
problems of ordinary calculus (15) and to solution of systems 
of algebraic equations (16, 17) and integral equations (18). 

The main idea of the present treatment stems from material 
presented by Prof. R. Courant in a 1941 address to the Ameri- 
can Mathematical Society (19). An application of the gradient 
method to fixed end-point variational problems has been 
given by Stein (20). Our extension of the gradient idea to the 
case which includes differential equations as subsidiary condi- 
tions is heuristic in character. 


Problem Formulation 


For present purposes it will be assumed that the system of 
differential equations to be satisfied along the flight path is 
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given in first-order form 


[1] 


These equations relate velocities and positions, forces and 
accelerations, mass and flow of propellants and coolants, and 
the like. z,, may be termed problem variables, and y the 
control variable. Differentiation with respect to the inde- 
pendent variable, time ¢t, is denoted by a superscribed dot. 

An important class of problems is that in which the per- 
formance quantity to be minimized is expressed as a func- 
tion of the final values of the variables z,, and ¢ 


P = .. t,) [2] 


At a specified initial time & as many as n boundary condi- 
tions on the z,, may be stipulated. Since an entire function 
y(t) is at our disposal, we may reasonably consider problems 
in which numerous conditions are imposed upon the xz, at 
various subsequent t values. In the following we will restrict 
attention to conditions imposed at the terminal point of the 
flight path. Among the n + 1 quantities consisting of the n 
final values of the 2» plus the final time t;, no more than n 
relations may be specified in order that the value of P not be 
predetermined. 


im = 


Neighboring Solutions—Variations 


We now assume that a solution of Equations [1] is available 
which satisfies the boundary conditions but which does not 
minimize P. Denoting the solution by 7m = £n(8), y = 5(0), 
we examine behavior in the neighborhood of this solution by 
setting tm = Fn + 6am, y = 5 + dy and linearizing 

The partials of g, are evaluated along x, = ¥mn, y = ¥ and 
are therefore known functions of the independent variable ¢. 
The functions 6z,, and dy are the variations of x, and y in the 
neighborhood of fn, 7. 

A formal solution of Equations [3] may be written in the 

form 


= t 
= 52 po&mp(t) + Mm(t, — r)dy(r)dr 
p=1 0 
m=1,...,n [4] 


where the first member represents solution of the homogeneous 
system of equations and the second a superposition of control 
variable effects. The functions u,, are Green’s functions or in- 
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fluence functions; u»(7, ¢ — 7) may be thought of as the 
solution for 6z, corresponding to éy a unit impulse (Dirac 
delta function) introduced at time 7 (21). 

Since interest centers on final values of z,, we evaluate the 
expressions [4] at ¢ = t;; however, to cover the possibility of 
a variable end point, t; = 7; + 6ét;, a first-order correction 
term must be included 


p=1 
p= 


Computation of the Functions z,, 


Since computation of the functions u,,(7, ¢ — 7) over a 
complete range of both arguments is unnecessary, only their 
evaluation at t = 7, required for subsequent calculations, it is 
reasonable to seek a means for performing the special compu- 
tation which avoids the labor of the more general one. The 
following development relates the functions u»(7, 7; — 7) to 
solutions of a system of equations adjoint to the system [3] 
through an application of Green’s theorem. The scheme em- 
ployed is due to Bliss, as reported by Goodman and Lance 
(22). 

We rewrite Equations [3] employing a subscript notation 
suitable to our immediate purpose 


=> + ay i=1,. 


[6] 
and write the system of equations adjoint to this system 
t=1,...,0 [7] 


which is obtained by transposing the matrix of coefficients and 
changing the sign. 
The solutions of the two systems are related by 


d n 


t=1 
After integration of both left and ne members between 
definite limits & and 7;, we find 


ay at [9] 


This is the one-dimensional form of Green’s theorem (22). 

We now consider numerical solutions of the adjoint system 
with all boundary values specified at t = i,. To the special 
solutions corresponding to 


Ady) = 
=] i= 


t=1 


we assign the symbols \,(”)(é)._ In this fashion n expressions 
for the values of the éz,,(z;) are obtained from [9] 

m = a [11] 


By comparison with Equations [5] it may be seen that the de- 
sired relation between the u»(7, 2; — 7) and the dA," is the 
following one 


= 


i=1 


Um (1, — = m=z=t1,...,” [12] 


and, also, that the ,.,(2Z/) of Equations [5] is equal to A," (t). 

In the preceding development the choice of symbols for 
the variables of the adjoint system is deliberate, for Equations 
[7] are precisely those governing the Lagrange multiplier 
functions of the “indirect” theory. We note the important 
distinction, however, that the coefficients of [7] employed in 
the “indirect” theory are evaluated along a minimal solution 
of Equations [1], whereas in gradient computations they 
correspond to nonminimal paths. 


Descent Parameter 


Following Courant (19), we now introduce a parameter o 1s 
a second independent variable, and seek a functional d»- 
pendence of the control variable y(t, «) on this parameter such 
that the derivative of the performance quantity dP/do is 
negative. In fact, within the restrictions imposed by tle 
boundary conditions and the system equations, we shill 
attempt to make the slope of descent dP/do “as steep 1s 
possible.”’ 

To enable operation within the restrictions just mentione1, 
we break down the control variable y as follows 


= ot, 2) + [13] 


Here the f,(t) are a set of known linearly independent func- 
tions of t. Our intention is that the coefficients a, of the 
second member of [13] be sacrificed to the fulfillment of 
boundary conditions, the number r being chosen appro- 
priately for this purpose. This will leave the function ¢ free 
for the minimization of P. 

We now take the derivatives of various quantities with re- 
spect to o and evaluate them at o = 4, corresponding to the 
nonminimal solution z,,(t, = In, y(t, ) = introduced in 
the preceding section. For ¢ — ¢ = Ago small, the varia- 
tions appearing in Equations [5] may be identified as 


Eg + (dim,/do) do 
m 


im = + my = + (d2m,/do) Ao 
m= 1,...,” [13] 


ty = ty + ot, =iyt+ (dt;/do) Ao {16] 
3+ by = + (Oy/0c)Ac 
$+ af,(t) + +> Ac 
q=1 do 


ll 


Equations [5] then take the form 


Om, 


Oa p= 1 
f _ at 


m=1,..,2 [18] 


Boundary Conditions 


We consider boundary conditions of the separated type, 
i.e., equations relating either initial values or final values. 
Boundary values may be variable on a surface, typified by 


Ley, ty) = O [19] 


in which case the following linear combination of derivatives 
must vanish 


d3 O35 O35 vy O35 dt; 
do Ot, do Ox, do + Ot; do 
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In this expression the partial derivatives of 5 are evaluated at 
= Xu, Ly = vy and ty = ty. 

In the case of fixed boundary conditions of the form zu, = 
f,, = constant, ty = 7; = constant, the relations to be satisfied 
take the simple form 


or 


do do [21] 


As many as n final conditions, fixed or of the form [19] 
may be specified, as mentioned earlier. Where fewer than this 
number are specified, we speak of ‘‘free” or “‘open’”’ boundary 
values. 

At the initial point similar freedom of choice may exist 
among the n initial values of the z,,, the only difference here 
being that ¢ will usually be fixed. 


Gradient of P 


Upon inspection it may be noted that there are n equations 
of the form [18] and at most 2n equations for boundary 
values, making a possible total of 3n equations. These relate 
2n + 1 derivatives of the xm,, Xm, and t,, the r derivatives of 
the a, and integrals containing 0¢/dc. We accordingly 
choose the number r of the a, asr = n — 1 — 8, 
where s is the number of open boundary conditions. Thus 
the system of equations will be determinate for arbitrary 
0¢, Oc. 

These equations may be arranged as a linear simultaneous 
system 


AZ=B [22] 


where Z has 8n — s elements consisting of the derivatives with 
respect to o Of 2m» Xm, t; and a,. The matrix A is square 
and contains among its elements the quantities Enp, Jm,, the 
3 partials, and integrals of products yunf,, all of which are 
known. The column B includes integrals containing 0¢/0¢ in 
their integrands. The solution 


Z=A"B [23] 


may be obtained through matrix inversion or equivalent proc- 
esses. The matrix A expresses the relationship between small 
shifts in boundary values and small adjustments in the con- 
trol function (through the a,) to deal with them. Hence in 
normal circumstances A will be nonsingular. 
The derivative with respect to o of the quantity P to be 
minimized 
dP‘ OP du, 


OP din, , OP dt, 
do dt; do 


[24] 


may now be expressed in terms of known quantities and inte- 
grals containing 0¢/0e according to the solution [23]. It 
will take the form of a linear combination of the integrals 


to Oo 
which may be expressed as a single integral 


|= Catin(7, ty — [26] 


do 


By analogy with a characteristic property of a vector gradi- 
ent, we are now prepared to identify the gradient of P with 
respect to the function ¢. [See p. 222 of the Courant-Hilbert 
English edition, (23), also (19).] If @ were a vector possessing 
a finite number of components ¢:, 7 = 1, .. ., 7, the derivative 
of P with respect to o 
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could be expressed as 


= = grad P. > [28] 


do do 


If ¢(¢, c) is a continuous function, and P a functional of ¢, 
as in the case of present interest 


dP _ 
Plo sear [29] 


and the function occurring in product with 0¢/d¢, which we 
denote [P]¢, may be regarded as the gradient of P by ex- 
tended definition. 

Thus for the problem at hand it is evident from Equation 
[26] that 


= Costas; ly 7) [30] 


is the gradient of P. 

Were the solution x, = fn, y = ¥ such as to minimize P, 
contrary to our assumption, the gradient [P], would vanish. 
If our development were to parallel the classical or ‘“‘indirect’’ 
approach, the construction of a solution would be sought from 
the vanishing of [P], for which P is stationary and possibly a 
minimum. Thus [P]¢ isin some sense an Euler expression and 
the equation [P], = 0 and the relations [23] have an equiva- 
lence to the Euler-Lagrange equations and transversality 
conditions of the “‘indirect’’ theory. 


Descent Process 


Returning momentarily to the elementary geometric con- 
cept of a vector gradient, we regard P(q,, . . ., @;) as a surface, 
and starting from P(g, . . ., ¢;) we move a point along this 
surface so that P and ¢; become functions of a time parameter 
a. Then the velocity of ascent or descent along a line on the 
surface is as given by Equations [27 and 28]. We now 
choose the line along which the descent is as steep as possible, 
characterized by 

do; oP 
do 0¢: 


where k positive corresponds to ascent and k negative to 
descent. 

It is clear that in this continuous process, wherein the point 
moves according to the system of ordinary differential equa- 
tions [31], the process will for ¢—> © approach a position for 
which grad P = 0, if P is bounded below. 

This elementary idea may be generalized to the present 
variational problem according to the extended interpretation 
of the gradient of P, Equation [29]. We set 


= —[Ple [32] 

and starting from the nonminimal solution ¢ = z,@ = ¢, tm = 
Em, Ag = Gq, re-evaluate these quantities on a continuous basis 
as the parameter o increases. 

Thus the continuous version of descent along the gradient 
requires numerical treatment of a partial differential equation 
for ¢(t, «) with determination of ¢(t, ~) the ultimate objec- 
tive. 


Stepwise Version 


As an alternative to the continuous procedure given by 
Equation [32], we may elect to proceed stepwise, correcting a 
set of approximations to the solution [P], = 0 by corrections 
proportional to the negative of the gradient 

git) [P]eAo [33] 
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Fig. 1 Orbital transfer schematic 


After choice of step size Ao and evaluation of ¢“+ and the 
the functions would be obtained from numeri- 
cal integration of the basic system Equations [1] and P“+ de- 
termined from the values at t;@+, 

Since the determination of the gradient [P], is expensive in 
terms of volume of numerical computations, it would seem 
desirable to exploit each calculation of local gradient direction 
of the utmost, taking Ao as large as possible. A procedure em- 
ployed in some applications (17) is to follow the local gradient 
direction until the function P reaches a minimum. Such a 
procedure could be implemented by numerical integration of 
Equations [1] for a number of values of Ao and selection of 
the value of Ao for minimum P. 

A possible pitfall of such a scheme of operation is that the 
boundary conditions, having been satisfied only in linearized 
version (Eq. [20]), may long since have been violated before 
minimum P is reached (24). Thus the choice of Ac involves 
a compromise which must be decided in the particular appli- 
tion at hand. If the minimum P rule is employed for choice 
of Ac, the boundary conditions must be restored by a correc- 
tion cycle designed to recover from the departures. 

After such restoration has been accomplished, for example 
via the coefficients a,, the solution so obtained takes on the 
role of ®, 7 for the computation of a new local gradient direc- 
tion. 


Convergence 


The convergence of the descent process has been investi- 
gated by Stein for a case not complicated by subsidiary condi- 
tions (20); his conclusion was that the process will converge 
if the functional whose minimum is sought is bounded below. 
In the process described in the preceding sections, the ques- 
tion of convergence is intimately related to the success of the 
technique for maintaining or correcting the boundary values 
at the terminal point. 

The possibility of correcting small departures in terminal 
values of the x, through small changes in the coefficients a, 
can be shown to hinge on the nonvanishing of the determinant 
of the matrix A of Equation [22]. In certain cases a tendency 
of this determinant to become small may be observed as a 
minimum of P is approached. Where the difficulty is acute, 
the behavior may be likened to breakdown of a first-order 
differential correction scheme for guidance along the trajec- 
tory. 

This type of behavior is a close relative of the conjugate 
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end-point phenomenon of classical variational theory. If an 
extremal through point 1 has on it a contact point 2 with the 
envelope of a family of extremals through 1, then point 2 is 
said to be conjugate to point 1. The term conjugate end points 
refers to specification of boundary values at two such points 
as 1 and 2. The reader is advised to consult the text of Bliss 
(25) for the analytical basis of conjugate point theory. 

The relevant test of the classical theory, the necessary con- 
dition of Jacobi, employs a determinant as a criterion. The 
requisite analysis amounts to investigation of optimal dif- 
ferential corrections of terminal values. Thus a matrix which 
is the optimal correction analog of the matrix A is set up, aid 
the criterion of the Jacobi test is the nonvanishing of the «- 
terminant of this matrix at points along the extremal up ‘o 
the terminal point of interest. Vanishing at the terminal poi it 
indicates that the end points are conjugates, and the necessa:y 
condition is fulfilled in borderline fashion. Computational 
procedures for applying the Jacobi test are of some intere:t; 
however, the acquisition of a “test specimen’’ solution requires 
that any convergence difficulty encountered in the desceut 
process first be overcome. 

Thus a possibility of convergence difficulty arises in conj.- 
gate end-point cases if the functions f, selected happen to 1e- 
semble the optimal correction functions of the indirect theory. 
There are certainly other possibilities for unfortunate choice 
of correction functions. Computational experience will |e 
required to establish guidelines on this matter. 


Configuration Parameters 


In many practical engineering applications, optimal per- 
formance is sought not only in terms of flight path selection 
but also in terms of parameters influencing vehicle configura- 
tion. We have, for clarity, avoided complicating the preced- 
ing analytical work by such considerations. It is an easy 
matter, however, to add terms to Equations [5] of the form 
(O2m,/e;)6e; and to carry them through along with con- 
straints relating e;, thus forming a basis for simultaneous 
optimization of configuration and flight path. The slopes 
O2'm,/Oe; are probably best determined by numerical integra- 
tion of Equations [1] for small changes in ¢;. 


Solar Sailing Example 


For the purpose of exploring the computational aspect of 
the gradient optimization technique, we have chosen a planar 
case of transfer between planetary orbits by means of the in- 
teresting solar sailing scheme. The potential capabilities of 
solar sail propulsion have been investigated in the papers of 
Garwin (26), Cotter (27), Tsu (28) and London (29). This 
problem has the simplicity appropriate to an exploration of 
method, yet sufficient complexity to render analytical solu- 
tion quite difficult unless drastic simplifications are intro- 
duced. 

The equations of motion and kinematic relations are given 
in a notation nearly the same as that of Tsu (28). With 
reference to the schematic of Fig. 1, these are as follows: 

Radial acceleration 


2 2 2 


Circumferential acceleration 


— — a(Ro/R)? sin cos? 6 [35] 
Radial velocity 
R = 93 = Uu [36] 
Circumferential angular velocity 
=0/R [37] 
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MARS ORBIT 
EARTH ORBIT 
R 
ner. 
SUN 
= R Ag R a R lcos 6| [34] 


] 


Since the heliocentric angle ¥ does not appear in the first 
three equations, and will not appear in the statements of 
boundary conditions to be considered, Equation [37] may be 
ignored for purposes of gradient optimization. This amounts 
to the assumption that terminal matching of the heliocentric 
angles of vehicle and “‘target’’ planet is accomplished by selec- 
tion of launch time. 

Seeking minimum-time transfer, we identify the functional 
Pas 

P = ty; [38] 


The functions u, v, R are the variables z,, of the theoretical 
development, and the sail angle @ appears in the role of the 
control variable y. 

As initial conditions we specify velocity components wu, v 
and radius R corresponding to motion in Earth’s orbit ap- 
proximated as a circle 


= 0 [39] 

u(0) = w = uz = 0 [40] 
= = ve [41] 
RO) = Ro = Re [42] 


We consider terminal conditions corresponding to arrival at 
the orbit of the planet Mars (also taken as a circle) with pre- 
scribed velocity components 


u(t,) = Uys [43] 
v(ts) = [44] 
R(t;) = Ry = Ru [45] 


For fixed boundary values of u, v and R, the equations cor- 
responding to Equations [5] of the preceding theoretical de- 
velopment are 


buy = + hi dts [46] 
iy 

= + = 0 [47] 

aR, = + = 0 [48] 


In this case the number of functions f, and coefficients a, 
required is 


r=n—-1—s=38-1-0=2 [49] 
We select the functions f, as 
fi =1 (50) 
fx) =? (51) 
and the control function 6(¢) is broken down as 
at + at? [52] 


The system of equations corresponding to Equation [22] 
simplifies to 


(f." + (f, + fi; = 


iy Og 


dr [54] 


iy da, iy daz , _ dt; 

0 


o 
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This 3 X 3 case may conveniently be inverted analytically. 
There seems little point in listing the inverse elements here, 
however. 


de f, 
das | _ i 

dts 
ido | JI Lvo 


The slope of descent dt;/do is given by 


dt i 
do 0 Oc 

and the gradient of P by 


= + Csoue + [58 
Accordingly we set 
= —[P]o = Cams + + [59] 


and proceed with stepwise descent as in Equation [33]. 

A particularly suitable case for a first illustration of com- 
putational technique is the one in which terminal velocity 
components are unspecified—‘“free’’ boundary conditions. 
Here Equations [46 and 47] may be deleted and 


r=n-—-l—s=3-1-2=0 {60] 
so that no functions f, are needed. Hence 
and 
OP us 
61 
Oa 93; I6 
us? dr [62] 
in this case. 


Computations have employed numerical values of the 
various constants from Tsu’s paper, with 


a = 0.1 cm/sec? = 3.28 X 107 fps? 


This value corresponds to about 10~4 g thrust acceleration 
developed by the sail when oriented broadside to the sun (6 = 
0) at Earth’s orbit radius, or about 17 per cent of the sun’s 
gravitational attraction. 

Results of descent computations for the case of “open” 
terminal velocity components are shown in Figs. 2 through 4. 
The control program of the original flight path (Figs. 2 and 3), 
chosen arbitrarily, was far from optimal in that the radial 
velocity component at crossing of Mars’ orbit was small. 
The greatest reduction in flight time—more than half of the 
original—is seen to be obtained in the course of the first de- 
scent (Fig. 4). In three descents minimum flight time has 
been attained for practical purposes, although small changes 
in the detailed structure of the control program are still in 
evidence. 

Results for the case of terminal velocity components 
matched to the target planet 


us = UM vs = UM 


are presented in Figs. 5 through 8. The tendency of the 
terminal values to depart from the prescribed values is shown 
in Fig. 5. These were restored via an iterative correction 
process employing increments in the coefficients a; and a» of 
Equation [52]. Typically, two or three iteration cycles were 
required to correct each point. Descent curves are shown in 
Fig. 6. The approach to the minimum-time solution is de- 
picted in Figs. 7 and 8. 
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Fig. 2 Successive approximations to optimal transfer path, 
terminal velocity components open 
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Fig. 3 Successive approximations to optimal sail angle program, 
terminal velocity components open 
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Fig. 4 Descent curves—solar sail transfer, terminal velocity 
components open 
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Fig. 5 Departure of terminal values, ‘‘matched’’ terminal 
velocity components 
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Fig. 6 Descent curves—solar sail transfer, ‘‘matched’’ 
terminal velocity components 
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Fig. 7 Successive approximations to optimal sail angle program, 
‘‘matched’’ terminal velocity components 
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Fig. 8 Successive approximations to optimal transfer path, 
‘‘matched’’ terminal velocity components 


The first attempt at computations for this “matched ve- 
locity’ case employed functions f, constant and linear with 
time. This met with near zero determinant difficulty of the 
type discussed earlier. The combination of linear and square- 
law corrections indicated in the foregoing was successfully 
used to avoid this difficulty. 


Concluding Remarks 


Attention has been confined in the preceding development 
to the main ideas of the gradient technique. The extension to 
cases involving several control variables offers no particular 
difficulty. The limited computational experience reported 
here suggests that the gradient technique may be a useful one 
in applications, and particularly in those presenting difficulty 
when the classical “indirect’’ approach—numerical integration 
of the Euler-Lagrange equations—is used. Apart from the 
feature of surmounting the two-point boundary value dif- 
ficulties of such cases, the gradient method may be particu- 
larly appropriate in distinguishing between minimal solutions 
of the Euler-Lagrange equations and those which are merely 
stationary. 
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Nomenclature 

Lm = problem variables 

y = control variable 

Ym = functions of zm and y appearing in basic system 
equations 

t = time 

P = function of final values of zm to be minimized 

6Ly,, bY = variations of rm and y 

Emp = functions appearing in the solution of Equations 
[3] 

Mm = Green’s functions in the solution of Equations [3] 

T = variable of integration 
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= variables of the adjoint system, Equation [7] 
= descent parameter 

= auxiliary control variable, Equation [13] 

= known linearly independent functions of ¢ 
= coefficients of the functions fy 


Ao = increment in o 

3 = function relating boundary values of zm and ¢ 

A = matrix of coefficients, Equation [22 

Z = column matrix having dzm,/do, drm ,/de, dt;/do 
and dag/do as elements, Equation [22] 

B = column matrix containing integrals of ym(0¢/dc) 
products, Equation [22] 

oe = coefficients of the linear combination of um, Equa- 
tion [26] 

[P]¢ = gradient of P, Equation [30] 

k = proportionality constant of descent 

ei = configuration parameters 

u = radial velocity component (see Fig. 1) 

v = circumferential velocity component (see Fig. 1) 

R = radial distance to sun (see Fig. 1) 

y = heliocentric angle 

Ay = acceleration due to sun’s gravitational attraction 
at Earth’s radial distance 

a = thrust acceleration of radially oriented solar sail 
at Earth’s radial distance 

6 = sail angle measured from radial orientation 

Subscripts 

i . = denotes final value 

0 = denotes initial value 

m, Prt general indexes 

= Earth 

M = Mars 

(-) = denotes time derivative 

(~) = denotes nonminimal solution 

(7) = denotes value at ith step 

(m) = denotes special functional notation, explained by 


Equation [10] 
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Optimum Thrust Programming of 
Electrically Powered Rocket Vehicles 


in a Gravitational Field 


Cc. R. FAULDERS! 


North American Aviation, Inc. 


Downey, Calif. 


The general problem of optimum thrust programming of an electrically powered rocket under 
the condition of constant jet power is considered. The thrust vector is assumed to be parallel to the 
instantaneous velocity vector at all times. The various optimization problems possible under these 
restrictions are shown to be equivalent to maximization of the change in velocity for specified pro- 
pellant mass and arbitrary range, maximization of range for specified propellant mass and arbi- 
trary change in velocity, aad maximization of the change in velocity for specified range and specified 
propellant mass. The calculus of variations is employed to obtain analytical expressions for the 
thrust acceleration program for the foregoing problems with a constant gradient of the tangential 


component of gravitational force. 


Limiting values of this gradient for which the tangential com- 


ponent of gravity can be assumed constant in the derivation of optimum thrust programs are deter- 


mined. 


eC aL powered rocket engines, such as the ion 
rocket and the plasma rocket, are characterized by a 
power source that is separate from the propellant. For a 
fixed power setting, therefore, the rate of propellant expendi- 
ture and the exhaust velocity can be varied over wide ranges, 
and a variety of thrust programs can be achieved. Electrical 
rockets are presently limited to very low thrust levels of the 
order of 10° or 10‘ Earth g’s per unit mass of the complete 
space vehicle. For this reason, an entire mission would 
generally be carried out under power, with operating times 
measured in days or weeks. 

The various possible requirements for optimum thrust 
programs are summarized in Table 1, assuming that the total 
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time of powered flight is a specified parameter, and that the 
thrust is always in the tangential direction. The change in 
vehicle velocity is indicated by AV, the change in position, or 
range, by As, and the propellant mass by mp. 


Table 1 Optimum thrust program requirements | 
(total time specified) 

Case AV As Mp | 

1 maximum arbitrary specified 

2 specified arbitrary minimum | 
3 arbitrary maximum specified 

4 arbitrary specified minimum | 

5 maximum specified specified | 

6 specified maximum specified 

“f specified specified minimum | 
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It is to be noted, incidentally, that instead of specifying 
time, the variable to be maximized or minimized in any of 
these cases could be specified and time could be minimized. 
Case 1, for example, thus would be altered to the problem of 
specifying velocity and propellant mass, leaving range arbi- 
trary and minimizing the time required to achieve the estab- 
lished objectives. 

Since only first variations of the appropriate functions are 
to be taken to obtain solutions to these problems, sufficiency 
conditions for true maxima or minima will be omitted from 
the mathematical formulation. In this sense, Case 2 is 
equivalent to Case 1, Case 4 to Case 3, and Cases 6 and 
7 to Case 5. Also, minimum time problems become equiva- 
lent to the corresponding reciprocal problems with time 
fixed. Therefore, three distinct sets of requirements can 
be designated as in Table 2. 

Case I corresponds to the classical problem of maximizing 
the burnout kinetic energy with a fixed total energy and pro- 
pellant mass. A practical application might involve the at- 
tainment of escape velocity with minimum propellant. utili- 
zation. Fox, Michielsen and Irving (1-3)? have analyzed 
Case I for zero gravity and found the required thrust program 
to be one of constant thrust acceleration. Case II can be 
likened to sending a scientific probe a maximum distance 
fron Earth or from the sun, with arbitrary final velocity. 
Thrust programming for maximum range is also discussed 
briefly by Irving in (3), again for zero gravity. A version of 
Case III wherein the increments in velocity and range are 
specified and the propellant is to be minimized is applicable 
to such problems as a change in position within a satellite 
orbit of specified radius or transfer between planets with com- 
bined acceleration and deceleration. 

Analytical expressions are derived here for the thrust 
programs corresponding to the three cases discussed, with the 
influence of gravity taken into account. The jet power is 
held constant at the maximum limit for reasons noted in 
(3). Although the analysis is restricted to trajectories where 
the thrust vector is always parallel to the velocity vector, 
the results are applicable to a large number of practical tra- 
jectories, as indicated in (4, 5, 6 and 10). In particular, 
Earth escape and some orbital transfer missions can be carried 
out with tangential thrust; normal components of thrust 
become significant only when steering requirements are con- 
trolling factors, as in certain Earth-Mars missions (7-9). 
For convenience in the analysis, the gravity force will be 
split into normal and tangential components (Fig. 1) with 
the normal component controlling the change of direction of 
the flight path and the tangential component adding, in alge- 
braic sense, to the thrust. The variation of this tangential 
component of gravity with distance along the flight path 
enters directly into the optimization problem. Although the 
total gravitational force varies inversely with the square of 
the radius from the center of mass, no a priori statement can 
be made as to the variation of the tangential component of 
gravity. An inverse square variation of this component 
would appear only if the trajectories were logarithmic spirals, 
and in general they are not. The assumption is made in this 
paper that the tangential component of gravity varies linearly 
with distance. An examination of digital computer data for 
powered trajectories between Earth and the moon with con- 
stant tangential thrust acceleration (4) indicates that actual 
variations of the tangential gravity component agree quite 
closely with this assumption. 


Variational Differential Equation and End 
Conditions 


Denoting the thrust per unit mass by a and the component 
of gravity force along the flight path by the derivative of the 
gravitational potential U with respect to distance, the equa- 


* Numbers in parentheses indicate References at end of paper. 
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Table 2 Distinct optimum thrust program requirements 
(total time specified) 


Case AV As Mp 
I maximum arbitrary specified 
II arbitrary maximum specified 
III maximum specified specified 


tion of motion can be written 
V=a+dU/ds {1] 
$=V [2] 
where dU/ds is a function of s only, and s is positive in the 
direction of the thrust vector. The propellant mass rate is 
related to the thrust acceleration by 


a = cm/(My — m) [3] 

and the constant jet power is given by 
P = (c*/2)tm [4] 
Combining Equations [3 and 4] to eliminate c and integrating 
According to Equation [5], the propellant mass is uniquely 

determined by the function 

y = 2P/(M, — m) [6] 


where 
[7] 
Adopting the Mayer formulation, the general problem de- 
fined by Cases I, II and III can be stated as one of determin- 
ing the thrust program a(t) for which the function 


F = G(Vr, 87, 7, T) +f." 
— V) + — Ja [8] 


is stationary. 1, v2 and v3 are unspecified functions of time, 
and G is that function of the variables at time 7 which is to 
be maximized or minimized. The three terms in the integrand 
that are set off by parentheses represent the constraints on 
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Fig. 1 Schematic sketch indicating force components 
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the problem, Equations [1, 2 and 7]. Writing the Euler- 
Lagrange equations 


for the integrand of Equation [8] results in 
=0 [10] 
n(d?U/ds*) + v2 = 0 [11] 
+ = 0 [12] 
v3 = 0 [13] 

Combining Equations [10 through 13] 

= (d?U/ds*)a [14] 


Boundary conditions for which the function F is stationary 
can be deduced by setting the total variation of F equal to 
zero. Thus, holding time 7’ fixed and recalling Equation [9] 


(SE +n) ave t+ aor + (22 x 


For Case I, where velocity is to be maximized while range 
is unspecified, G becomes V7 and és7 is nonzero. Thus, from 
Equation [15], (v2)7 is zero or, from Equations [10 and 12] 


In addition, 6V7 is nonzero. The requirement that the co- 
efficient of 6Vr be zero, however, contributes no useful in- 
formation. 

For Case II, maximum range and arbitrary velocity, G be- 
comes s7 and éVrisnonzero. Again from Equation [15], (»1)7 
is zero and, from Equation [12] 


ar =0 


In Case III, G again becomes V7, but only the coefficient 
of 5Vr can be equated to zero owing to the specified range. 
As with Case I, however, the stipulation that (»; + 1)7 equals 
zero does not impose any end condition on a(t). This result 
would be expected, for Equations [1,2,7 and 14] comprise 
a fifth-order set of equations. Sufficient boundary condi- 
tions for the solution follow from V, s and y at time zero and 
the specified s and y at time 7’. 


Linear Variation of Tangential Component 
of Gravity 


Case I 
Defining 
= @U/ds* [18] 
Equation [14] becomes 
a = aa [19] 


Recognizing the appropriate end condition for Case I, 
namely Equation [16], the solution of Equation [19] is 
cosh aT (1 — t/T) 


a 
cosh aT’ [20] 


Substituting Equation [20] in Equation [5] and solving for 
the initial acceleration 


V4aT [(w — 1)/w] 


sinh aT [21] 
EM, cosh aT cosh? aT 


According to Equation [18], a constant tangential com- 
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ponent of gravity corresponds to a of zero. Taking the limits 
of Equations [20 and 21] as a approaches zero 


a/a = 1 [22] 


— 1) 
23 
V P*/EMr \ ( w ) 
The equation of motion can be obtained from Equations 
{1 and 2] in the form 


§ — a’s = a + (dU/ds)o — a8 [2-4] 


With the thrust acceleration given by Equation [20], tle 
solution to Equation [24] can be obtained by means of stan:- 
ard techniques applicable to linear differential equatiors. 
The velocity can then be obtained by taking the derivative of 
s with respect totime. Thus 


av. aT Vo 
(ao/2a) (sinh =) (a,/2a) 
(ao/2) 
As = cosh? aT — 1 Vo : 
(ao/ 2a?) ( cosh aT ) + (ao/2a) 
(dU/ds)o 


(d/2) (cosh aT — 1) [26] 


where @ is given by Equation [21]. 
Case II 


With the end condition given by Equation [17], the 
solution of Equation [19] applicable to Case IT is 


@ _ sinh af(1 — ¢/T) [27] 
A sinh aT 
From Equation [5] 
= aT aT [28] 
EMr sinh aT sinh? aT 


With an a of zero, Equations [27 and 28] reduce to 


= T [29] 

Vol ) (30) 


The solution of Equation [24] and subsequent differentia- 
tion of the solution with respect to time leads to the ex- 
pressions 


= sinh aT’ + @/2a) (cosh a7 — 1) + 
(dU/ds)) . 
sinh aT [31] 
As aT Vo 
) = (c aT — aT + 
(dU/ds)o 
(ao/2) (cosh aT — 1) [32] 


where dp is given by Equation [28]. 

A comparison of the first two terms of Equation [25] with 
the first term of Equation [31] indicates that, for zero initial 
velocity and zero initial tangential gravity force, the velocity 
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increment for Case I is greater than that for Case II by the 
ratio 


_ V (cosh? aT’)(sinh? a7) — (aT)? 
AVu sinh? aT 


[33] 


where the values of a» for Cases I and II have been obtained 
from Equations [21 and 28], respectively. At a7’ of zero, 
Equation [33] indicates a ratio of velocities of 2/+/3, whereas 
at a7’ approaching infinity, the velocity increments for 
Cases I and II are the same. At large a7’, therefore, maxi- 
mum velocity increment and maximum range are realized 
simultaneously. 

‘The ratio of the thrust contributions to the distance incre- 
ments for Cases I and IT is given by 


As; sinh? a7 
den ~/(eosh? aT) (inh? aT) — (aT 
or the inverse of the ratio of velocity increments. 
Case III 
‘he general solution of Equation [19] is simply 
a = C, cosh aT + C2 sinh aT [35] 
which reduces to 
a= C, + Cat [36] 


as w approaches zero. In this case, the constants of integra- 
tion C; and C2 are determined from the specified range and 
propellant mass. Equivalently, specified values of range 
and velocity increment could be used to determine the con- 
stants. A problem that would fall under the category of 
Case III is one wherein the velocity is to be maximized but 
the range is specified as some function of the velocity. 


Selection of Typical Values of (aT) 


Optimum thrust programs for the condition of a constant 
gradient of the tangential gravity component (Eqs. [20, 27 
and 35]) are dependent on the dimensionless parameter a7’. 
This parameter, essentially the square root of a nondimen- 
sional gravity gradient, is a function of the particular tra- 
jectory characteristics, such as the initial position and velocity 
vector, the initial thrust to weight ratio of the vehicle, the 
thrust acceleration program, and the time of flight or dis- 
tance of travel. It is of interest to determine representative 
values of aT’ for some actual powered flight space trajectories. 

For spiral trajectories from an Earth satellite orbit to the 
moon, or to Earth escape velocity, values of a7 can be 
roughly estimated. For a constant thrust per unit mass 
that is small compared with the local acceleration of gravity, 
the initial tangential component of gravity is twice the initial 
thrust per unit mass, and the final tangential component of 
gravity is approximately zero (4). Therefore 


Vo?T? 
roAs 


or 


[37] 


Where Vso is the initial satellite velocity. Range is related 


to Vso by 
= 2As 
and the increment in radius is 


Ar & As(sin ¢)ave 
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Table 3 Values of a7’ for Earth escape 


(u/To?) (Ar/ro escape (aT 

10-4 100 0.1 | 
10-3 30 0.5 | 
10-2 10 1.0 

10-1 5 3.0 


* Approximate value. 


Table 4 Values of a7’ for small radius ratios 


(r2/r1) 
1.25 0.4 
1.5 0.5 
2.0 0.75 


« Approximate value. 


Again from (4), sin ¢ varies between approximately zero at 
takeoff from a satellite orbit and approximately unity at escape. 
With an average value of sin ¢ of the order of one half, 
Equation [37] becomes 


(u/To”) To 


Current estimates of the thrust acceleration of electrically 
powered rockets lead to values of the order of 10~* to 10~* 
Earth g’s. The ratio appearing in Equation [38], however, 
is a function not only of the ratio of thrust to sea-level 
weight, but also of the radius of the initial satellite orbit. 
For values of a/(u/7o?), estimates of Ar/ro for escape were 
obtained from (4). Values of a7 computed from Equation 
[38] are indicated in Table 3. 

Estimates of a7’ can also be obtained for a tangential 
thrust trajectory extending from the Earth orbit to other 
nearby orbits, such as that of Mars, and with only the solar 
gravitational field considered. The initial thrust to mass 
ratio is assumed to be 10-4 Earth g’s, corresponding to a 
thrust to mass ratio of the order of one tenth the local solar 
gravity. In this case, the initial tangential component of 
gravity is again of the order of twice the initial thrust per 
unit mass. However, the final tangential component of 
gravity may be as large as one half this value, assuming that 
the thrust is always directed away from the starting point. 
Thus, following the development of Equation [37] 


w  Vso°T? 
TAs 


= 16 


[38] 


[39] 
The expressions for range and increment in radius become, 
in this case 

Voor = As 
From (4), for small thrust to weight ratios 


(sin 2ao/(u/To?) 


Ar = As(sin ¢)o 


Finally, therefore 
1 Ar 
aT? = 2 [40] 
On the basis of the foregoing analysis, the correspondence 
between values of a7’ and radius ratios of the order of unity 


was taken to be that shown in Table 4. 
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In Fig. 2, the optimum nondimensional thrust programs for 
Case I, computed from Equation [20], is plotted for values 
of aT from 0 to 3.0. Thrust programs for Case II, Equation 
[27], are plotted in Fig. 3 for the same values of a7’. Initial 
thrust accelerations are given in terms of powerplant char- 
acteristics by Equations [21 and 28]. 

A worthwhile comparison can be made of the velocity 
increment resulting from the use of the optimum thrust ac- 
celeration program with the velocity increment that would 
be obtained with simply a constant thrust acceleration. 
Similarly, the range obtained with the optimum program 
can be compared with that for constant thrust acceleration. 
Power level, operating time and mass ratio will be assumed 
to be the same in each case. The results of such comparisons 
are indicative of the importance of determining an optimum 
thrust program. Integration of Equation [24] with a con- 
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Fig. 2 Optimum thrust acceleration programs for Case I 
(maximum velocity increment) 
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Fig. 3 Optimum thrust acceleration programs for Case II 
(maximum range) 


stant thrust acceleration leads to 


sinh a7 [41] 
As Vo 
2(cosh aT — 1) + sinh aT + 
a (cosh aT — 1) [42] 


where dp is given by Equation [23]. For zero initial velocity 
and initial tangential gravity force, Equations [21, 23, 25, 2s, 
32, 41 and 42] can be combined to give 


AVi _ (sinh eT)(cosh aT) + aT 
AV 2 sinh aT 


[43 


Asii V (sinh aT)(cosh aT) — aT 
ASa=ay 2 cosh aT’ — 1 

These two ratios are plotted in Fig. 4. With an aT of, 2.0, 
Fig. 4 indicates a 9 per cent increase in velocity and a 24 
per cent increase in range owing to replacement of a con- 
stant thrust acceleration with the optimum program. At 
large values of #7’, both ratios approach WaT/2. 


[44] 


Discussion of Results 


According to Fig. 2, the assumption of a constant tangential 
component of gravity is adequate for determination of thie 
optimum thrust program for Case I when aT’ is less than 
0.1. With Case II, however, it can be seen from Fig. 3 that 
the same assumption is valid for a7’ up to about 0.5. Opti- 
mum thrust curves for Case III might be expected to follow 
the same general variation with a7’ as that exhibited in Figs. 
2 and 3. A safe estimate for the maximum value of a7’ for 
which the gravity force gradient can be neglected might be, 
therefore, of the order of 0.1 for Case ITT. 
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Fig. 4 Velocity increments and range obtained with optimum 
thrust program compared with values resulting from use of 
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The parameter a is real, and Figs. 2 and 3 are applicable 
when the gravity force is in opposition to the thrust and de- 
creases in the positive direction of motion, or when the gravity 
force adds to the thrust and increases in the positive direction 
of motion. The latter condition would arise, for example, 
if a vehicle were to be accelerated toward the origin of a 
central force field. With a gravity force opposing the thrust 
but increasing in the direction of flight, or adding to the 
thrust but decreasing in the direction of flight, a would be 
imaginary and the thrust variations with time would be 
sinusoidal with a period of 27r/a. There are no obvious 
situations, however, other than the oscillatory phase of low 
thrust takeoff from a satellite orbit, where an imaginary a 
would be encountered with gravity steering alone. 

‘The decrease in thrust acceleration to zero corresponding 
to maximum range (Fig. 3) is accompanied by an increase 
in the exhaust velocity to infinity, as can be seen from Equa- 
tions [3 and 4]. In practice, a limiting exhaust velocity 
wil! exist, the value depending upon the design of the rocket 
engine. The thrust acceleration program to be used when 
account is taken of the minimum thrust limitation is, there- 
fore, of interest. A complete analysis of this problem is 
beyond the scope of this paper. It can be shown, however, 
that a stationary solution is obtained with the thrust accelera- 
tion given by Equation [14] and the end condition, Equation 
[17], up to the time when the limiting exhaust velocity is 
reached and then holding constant exhaust velocity to burn- 
out at total time 7. The end condition would apply at the 
end of flight, not at the time when the limiting curve is 
reached. That is, the curve of thrust acceleration prior to 
reaching maximum exhaust velocity would, if extended past 
the limiting curve, go through zero at the time maximum 
range is to be achieved. Still greater range could probably 
be obtained by suddenly changing the power to zero or some 
reduced level at the time the limiting exhaust velocity is 
reached. Even with a discontinuity in power, however, the 
foregoing remarks concerning the variable exhaust velocity 
portion of the thrust acceleration curve are still applicable. 

A requirement more common than maximizing distance of 
travel, or range, would be that of maximizing the change in 
radial position from the center of the gravitational field. 
When the flight path extends out beyond the radius where 
escape velocity is attained, a correspondence can be expected 
between maximum radius and maximum length of flight 
path. For flights that begin and end deep within the gravi- 
tational field, however, the reduction in thrust acceleration 
to zero that is necessary for maximum range may result in 
a tightening of the spiral trajectory to the point where the 
radius increment is much less than could be obtained with 
some other thrust program. A brief analysis of the latter 
case is given in the Appendix. It is shown there that when 
the thrust acceleration is small compared with the total 
gravitational acceleration at all times, and when the vehicle 
is initially in a circular satellite orbit, maximum change in 
radius is obtained with a constant thrust acceleration. 

The assumption of a constant gradient of the tangential 
gravity component does not apply, certainly, to all possible 
trajectories that can be generated with tangential thrust. 
For more complex variations of the tangential gravity com- 
ponent with distance (e.g., parabolic, exponential), analytical 
solutions can be obtained for specific cases by expansion of 
Equations [1 and 14] in powers of some appropriate non- 
dimensional parameter. Regardless of the actual variation 
of the tangential gravity component, the second derivative 
of thrust acceleration with time is positive as long as the 
thrust remains in the positive direction and a is real, as 
indicated by Equation [14]. End conditions for Cases I 
and II, furthermore, are independent of the gravitational 
force. On this basis, the thrust programs of Figs. 2 and 3 
would be expected to be reasonable approximations for non- 
linear variations of the tangential gravity component as 
long as the average slope is used in the determination of aT’. 
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To obtain an exact solution, numerical integration of 
Equations [1, 2 and 14] would be necessary. For this pur- 
pose, an iterative scheme might be used wherein a trajectory 
would first be computed with an assumed thrust program, 
the resulting variation of the tangential gravity component 
would be used in the computation of a revised thrust pro- 
gram, and the process repeated until convergence was 
realized. The results presented here should prove useful 
in determining the thrust program to be assumed for the 
first iteration. 


Conclusions 


The optimum thrust acceleration program for the case of 
maximum velocity increment with a specified propellant 
mass and arbitrary range is constant thrust acceleration when 
the tangential component of gravity is constant. With a 
variable gravity component, the variation of the thrust 
acceleration involves a zero time derivative at burnout. 

For maximum range, arbitrary velocity increment and a 
specified propellant mass, the thrust acceleration decreases 
linearly to zero when the tangential gravity component is 
constant. For arbitrary variation of the gravity component, 
the burnout thrust is zero. 

The assumption of zero gravity in the determination of the 
optimum thrust program is valid for most engineering pur- 
poses as long as the gravity component variation can be 
represented by a straight line with the square root of the 
nondimensional gradient (a@7’) of the order of 0.1 or less. 


Appendix: Thrust Acceleration Program 
for Maximum Change in Radius 


If the thrust acceleration is small compared with the total 
gravitational force at all times, and if the vehicle is initially 
in a circular satellite orbit, the equations of motion can be 
linearized and the solution averaged over a single revolution 
(4) to give 


V=-—-a [A-1 
dr a 


From Equation [A-2], maximum change in radius is obtained 


when the integral 
T 
ava [A-3] 


is maximized. With a specified mass of propellant, the 
integral 


~ 
ll 


f [A-4] 


must be held constant, according to Equation [5]. From 
Equation [A-1] 


f V Vat [A-5] 
J= f Vdt [A-6] 
The function 
f=VV+,)AP 
must satisfy the Euler-Lagrange equation. Thus 
_d of _ 
dtoV 
or 
= 0 
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From Equation [A-1] therefore 
a=0 [A-7] 


and the thrust acceleration is a constant. 


Nomenclature 

a = thrust per unit mass of vehicle 

c = exhaust velocity, measured with respect to vehicle 

E = total energy expended, EF = PT 

f = integrand for which Euler-Lagrange equations are 
written 

F = stationary function 

G = function of dependent variables at time 7 

I = integral defined by Equation [A-3] 

J = integral defined by Equation [A-4] 

m = mass of propellant expended at time ¢ 

Mp = total mass of propellant 

My = initial mass of vehicle plus propellant 

Mr = final mass of vehicle, M) = Mr + mp 

P = jet power 

r = radius measured from center of attracting body 

s = distance measured along flight path 

t = time 

T = total time of flight 

U = gravitational potential 

V = velocity of vehicle 

Vso = satellite velocity at initial radius 

a = Vd? U/ds? 

r = Lagrangian multiplier 

m = constant of a particular gravitational field 


Solid Propellant Detonability 


v1, v2, vs = unspecified functions of time ¢ 


¢ = flight path angle measured from circumferential 
direction 

= 2P/(M, — m) 

w = mass ratio, w=M)/M,r 

Subscripts 

0 = initial value at ¢ = 0 

: if = final value ati = T 

I = Case I 

II = Case II 
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Using a shock sensitivity test, it has been possible to formulate broad patterns of behavior for 


propellants under certain conditions. 


Conventional double-base propellants have a sensitivity 


somewhat less than that of conventional insensitive military explosives. This sensitivity increases 
slightly with an increase in temperature and markedly with decreasing density. Pore type is im- 
portant; samples with nonconnected pores increase in sensitivity as temperature decreases. Many 
composite propellants containing conventional rubber-like binders are not easily detonated. Com- 
posites containing high energy binders are readily detonated. An investigation of the attenuation 


of shock by Lucite has been made. 


The results of this study make it possible to determine the 


pressure required to initiate detonation in a propellant under the conditions of the booster sensi- 


tivity test. 


S THE quantity of solid propellants used in individual 
rockets, particularly in guided missile propulsion systems, 
grows larger and larger, the safety of such motors with respect 
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to accidental explosion or detonation becomes of considerable 
concern. Because of the confines of shipboard usage, the 
Navy must give special attention to this aspect of safety. 

Early in the Polaris program projects to study propellant 
safety were established with a number of organizations, among 
them the Naval Ordnance Laboratory. One aim of the NOL 
project was to determine the shock sensitivity of solid prope|- 
lants. 

Experience with explosives has shown a good correlation «f 
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impact test results with field handling experience and also with 
the results of shock sensitivity tests, also known as booster or 
gap tests (1).4 The simple and inexpensive nature of impact 
testing has been an advantage in the preliminary screening of 
explosives. Application of the test to propellants has failed 
to produce results which can be correlated with handling ex- 
perience or with shock sensitivity tests. As a matter of fact, 
since the impact sensitivities of most propellants imply that 
they are more sensitive than most booster explosives, it is ap- 
parent that the use of these results as a basis for screening 
would have denied the Armed Forces most, if not all, of its 
propellants. 

For this reason and, in addition, to evaluate directly the 
shock induced detonability of solid propellants, we have used a 
gap test which is a modification of one generally used for liquid 
monopropellants and developed by Edwards at this Labora- 
tory (2). Briefly, gap testing is a procedure by which a high 
energy shock wave is transmitted from a “donor” to an 
“acceptor.” In its usual form, using a standardized donor, 
the distance between the two is varied and the results are re- 
ported in terms of that distance or “gap” at which there is a 
50 per cent probability of detonation of the acceptor. 

Fig. 1 represents, schematically, the fashion in which the 
tests are conducted at NOL (3). The gap cards are of cellu- 
lose acetate and are each 0.0254-cm (0.01 in.) thick and 5.08 
em (2 in.) in diameter. For gaps greater than 1.27 cm 
(0.50 in.) the cards can be replaced by a combination of 1.27- 
cm (0.50 in.) Lucite disks and 0.0254-em (0.01 in.) cards to 
give the desired thickness. We have found no discernible 
difference in substituting 1.27 cm (0.5 in.) of solid Lucite for 
50 cards. The charge diameter exceeds the critical diameter® 
for most explosives. 

The damage to the plate is used to determine whether a det- 
onation has occurred. With most explosives a detonation 
very neatly punches a hole into the center of the plate; the 
diameter of the hole is close to that of the charge. The re- 
mainder of the plate is essentially undamaged. In the 
absence of a detonation, no damage to the plate occurs. With 
propellants, the critical diameter appears often to approach or 
to exceed that of the charge used. In such systems misleading 
(false positive) results may be produced by an unstable detona- 
tion. For such systems it is good practice to measure shock 
velocities. Methods for doing this are described elsewhere in 
detail (4). Briefly, thin wires are inserted into the charge at 
predetermined points. Connected to these wires are units 
which measure the time required for the shock or detonation 
to travel between pins. The stability of detonation is judged 
by the constancy of the velocity as the shock propagates along 
the charge. We have been able to correlate the plate damage 
very nicely with the shock velocity at the plate. For veloc- 
ities over about 4000 m per sec, the damage is essentially the 
same as with high explosives. With increasing detonation 
velocity, the hole size increases up to about 23 in. with 
Composition B, which has a detonation velocity of about 7800 
m per sec. Within the range 3000-4000 m per sec, the edges 
of the plate around the hole will show an outward bulge. From 
1200 to 1300 m per sec there will be a derby hat-like bulge in- 
stead of a hole, the depth being larger for higher velocities. 
At less than 1200 m per sec there is no discernible effect. 
Since for all safety purposes the detonation velocity is sec- 
ondary, all stable detonations are reported as “‘goes.”’ 

The test procedure is roughly as follows: A test is run at 
zero gap. If detonation occurs, the next test is made at eight 
cars, and the number of cards is doubled with each additional 
test until the sample gives a negative result. Each successive 
shot is then made half way between the closest “go” and 
“no-go” card values. This procedure is continued until a 
positive and negative test have been obtained which differ by 
only one card. Without going into further detail, it takes 

* Numbers in parentheses indicate References at end of paper. 

4 l.e., that diameter below which a detonation will not propa- 
gate. 
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Fig. 1 Shock sensitivity test apparatus 


only two to four more shots to ascertain within one half card 
what the 50 per cent point is. It is our feeling that the re- 
sults reported herein are precise within one to two cards. 


Test Results 


Impact vs. Gap Testing 


Let us now examine some of the results which have been 
obtained. First let us take a look at a comparison between 
the impact sensitivities of explosives and of some typical 
propellants. Bearing in mind that Tetryl and pure RDX are 
considered to be too sensitive for use as high explosives, it is 
seen immediately from Table 1 that none of the propellants 
tested would have been accepted by this‘standard. Double- 
base and, by now, many composite propellants have a history 
of safe use and freedom of accidents. There do not seem to be 
any reports of trouble during handling of these materials after 
they are through the manufacturing stages. Safety records of 
this sort are hard to reconcile with the impact data; it is well 
not to place too much reliance upon them. There is, however, 
one important area in which impact data may be very useful. 
Research at this laboratory has shown a good correlation be- 
tween impact data and ignitibility test results (5). 


Table 1 Approximate order of impact sensitivities of 
some common explosives and propellants 


Propellant 
composite high energy binder 
(11) 
composite polysulfide rubber (15) 
composites polyurethane binder 


Explosive 
PETN (11-167) 


RDX (16-33) 


(16-41) polyvinyl chloride 
binder 
polyester-styrene 
binder 


Tetryl (26-50) double-base Type I (28) 
Composition B (45-82) 
TNT (150-215) 

* Numbers in parentheses indicate height in cm, using 


Bruceton Type 12 tools, 2.5 kg weight. 
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Shock Sensitivity of Nondefective Solid Propellants 


A number of different types of double-base and composite 
propellants have been examined for shock sensitivity. The 
results are listed in Table 2; values for some better known 
military explosives are given for comparison. It is seen that 
conventional composites (polyurethane, polyvinyl chloride or 
polysulfide rubber binders) cannot be detonated in the diam- 
eters tested. The double-base propellants detonate readily, 
and the use of an explosive binder increases the sensitivity al- 
most to that of our common insensitive high explosives. 

The data in Table 2 were taken at or close to 24C. It is 
interesting to examine the effects of temperature. Typical 
data are presented in Table 3. Obviously marked changes in 
temperature do not induce any comparable change in sensi- 
tivity, and such change as does occur is generally in the ex- 
pected direction; i.e., rising temperatures increase detona- 
bility. All propellants which failed to detonate at 24 C also 
failed at —59C and +66C. Furthermore, repeated tempera- 
ture cycling between these extremes had no apparent effect on 
any of the propellants. An indication of the insensitivity of at 


Table 2 Shock sensitivity of selected solid propellants 
and standard military explosives (temperature = 24 C = 
75 F 


Sensitivity 
(cards at 50% 
Material point) Remarks 
Propellant 


polysulfide rubber binder <0? 
polyvinyl chloride binder <0 


polyurethane binder <0 
polyester-styrene binder <0 
double-base Type I (me- 
dium energy ) 93 
double-base Type II 
(high energy) 35 
composite with high en- 
ergy binder 72 
hybrid (double base plus 
inorganic oxidizer ) 88 
Explosive 
ammonium nitrate <0 pressed 1.61 g/cc 
Tritonal 114 cast 
TNT 1393 cast 
Explosive D 1503 pressed 1.59 g/cc 
Composition B 1793 cast 
Tetryl 272 pressed 1.63 g/cc 


* Did not detonate. 


Table 3 Effect of temperature on propellant sensitivity 


(—59 C) (—32C) (—4 C) (24 C) (38 C) 


Propellant 


double-base 

Type I <0 <0 23 94 12} 
double-base 

Type II 263 28 34 35 374 
composite 

with high 

energy 

binder 49 72 973 
hybrid 

(double 

base plus 

inorganic 

oxidizer ) 57 Pit see 88 82 


“In cards 0.0254-cm (0.01 in.) thick. 


least two of the composites is indicated by the failure of two 
polyurethane binder systems to propagate a detonation in a 
modified gap test, i.e., a charge 50.8 cm (20 in.) in diameter 
and 203-cm (80 in.) long confined in a 0.635-cm (0.25 in.) thick 
steel tubing and boostered by a 340.5-kg (750 lb) Composition 
B charge 50.8 em (20 in.) in diameter and 101.6-cm (40 in.) 
long (6). 


Shock Sensitivity of Porous Solid Propellants 


Defective grains (containing pores, cracks, etc.) are not:- 
rious for exploding during or subsequent to ignition (7). In 
order to further our understanding of the process so that con- 
trol might be accomplished, the shock sensitivities of a numb:r 
of porous propellants were measured. Two types of por:s 
were treated: Connected and nonconnected. 

The connected pores were made in either of two ways. With 
the double-base propellants small cylindrical granules (1 x 4 
mm and 2 X 8 mm) were extruded, and these were poured in‘o 
the steel gap test tube until, after moderate shaking, the tu!ve 
was filled. For the composite propellants, charges were pre- 
pared by pressing fine lathe turnings into the steel tube. ‘n 
each case the charge density was determined simply ly 
weighing the charge used and knowing its dimensions. Tlie 
resiit was compared with the density of a solid charge and ex- 
pressed as “% M.D.,” i.e., percentage manufactured density. 

For double-base and composite charges containing noncon- 
nected pores, samples were obtained by manufacturing under 
conditions conducive to bubble formation. 

The results for a typical double-base propellant (Type [, 
see Table 2) are shown in Fig. 2. It is seen that sensitivity in- 
creases with increasing porosity (i.e., decreasing % M.D.). At 
ambient temperature and above, pore type does not seem to be 
important. As the temperature is lowered pore type becomes 
increasingly significant. For sampleswith connected porosity, 
sensitivity decreases with temperature. For those with non- 
connected pores, sensitivity increases remarkably with de- 
creasing temperature. 

A different type of behavior is noted with the composite pro- 
pellants so far tested. Here, as with the double-base systems, 


SENSITIVITY CARDS 


20+ 

TEMPERATURE (°F) 

40 2 60 
TEMPERATURE (°C) 


Fig. 2 Effect of porosity and temperature on the shock sensi- 
tivity of a medium energy double-base propellant (cp = connected 
pores; ncp = nonconnected pores) 


ARS JouRNAL 


il 
t] 
d 
n 
is 
Vi 
cl 
ra 
te 
ra 
; Al 
or 
=o be 
pl 
200+ d 
563% MD. 
: po 
170+ ed 
to 
nic 
140 
ou 
80) | 
pal 
50- tes 
sh 
tra 
wh 
fro 
afi 
ste} 
962 


wae 


the introduction of connected porosity increases sensitivity; 
temperature changes have a secondary effect—increasing tem- 
perature increases sensitivity slightly. Samples containing 
nonconnected pores are, within the range tested, nondetonable 
regardless of % M.D. or temperature. Although they 
have not yet been tested, one would expect that porous com- 
posites containing an inherently explosive binder would behave 
in the same fashion as the double-base systems. 

Is there an explanation for the marked increase in sensitivity 
induced by decreasing the temperature of the foam-like 
double-base propellants containing nonconnected pores? Can 
this behavior be reconciled with the insensitivity of the foam- 
like composites even at low temperatures? We believe both 
these questions can be answered. Let us examine a postulated 
mechanism for the shock induced transition to detonation in a 
confined granular material (8). Presumably a deflagration 
begins at the exposed surfaces of the propellant fragments. 
Since the sample is confined, the product gases Co not escape, 
ani the rising pressure causes an increase in the linear burning 
rate. The process continues until a stable detonation de- 
velops. Obviously the rate at which the pressure rises de- 
pends upon the total area of the propellant exposed to the hot 
gases initially present. In a system containing connected 
pores the gases can diffuse readily and a large area can ignite; 
the result is an easily propagated deflagration; i.e., the sample 
is sensitive. If the pores are nonconnected, the hot gases do 
not readily diffuse forward (away from the reaction zone into 
unreacted material), and, until high pressure gases rupture the 
interpore membrane, the reaction is restricted. If, however, 
the sample is brittle the same shock which serves to initiate the 
deflagration may be strong enough to produce an almost im- 
mediate rupture of these membranes, and the sample will now 
behave in an extremely sensitive fashion. For now, not only 
is gas flow possible, but the total surface area available for 
combustion is increased manyfold. The brittleness of any 
viscoelastic material (such as a propellant) is known to in- 
crease with decreasing temperature and/or with increasing 
rate of stress application (9). The double-base systems 
tested are much less elastic at ambient temperature and slow 
rates of stress application than are the composites tested. 
This brittleness apparently is the cause for the sensitivity. 


Absolute Scale for Expressing Sensitivity 


As convenient as the card gap values are for expressing det- 
onability on a relative scale, it certainly is more desirable to 
be able to express sensitivity in terms of a more fundamental 
scale, e.g., rate at which energy must be transferred to an ex- 
plosive to initiate detonation. Prerequisite to the complete 
determination of these data, and perhaps the most significant 
portion thereof, is knowledge of the pressure at the leading 
edge of the shock wave. Although it has not yet been possible 
to measure such high pressures directly, some indirect tech- 
niques do exist (10). One such method has been applied to 
our gap test where we were interested in determining the pres- 
sure in the leading edge of the shock wave as it left the Lucite 
gap (attenuator) and was about to enter the acceptor. The 
following steps are involved: 

| High speed photographic studies are used to obtain the 
particle velocity (u) vs. shock velocity (U) relation for Lucite. 

2 Using the same donor-attenuator geometry as in the gap 
test, simple probe measurements are used to measure the 
shock velocity as a function of the distance the shock has 
travelled from the tetryl-Lucite interface. 

3 Using the relation for shock pressure 

P = pul 
Where pp is the density of the unshocked material, and the data 
from steps 1 and 2, we can then calculate the shock pressure as 
a function of distance into the gap. 

The data required for step 1 are available (11), and those for 
step 2 have been obtained at this laboratory and will be re- 
ported in detail at a future date. Of particular interest and 
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Fig. 3 Pressure vs. attenuator thickness in the shock sensitivity 
test 


importance to the propellant sensitivity problem are the data 
from the completed calculations. These are shown in Fig. 3, 
in which the pressure is shown as a function of distance into the 
attenuator. The sensitivities of typical materials are indi- 
cated. 


Conclusion 


It is demonstrated that impact test data do not correlate 
with detonability measurements. Despite the inexpensive 
and, therefore, attractive nature of the test, impact test 
results should not be used as criteria for the rejection of in- 
teresting systems. 

Shock sensitivity data demonstrate that propellants tend to 
be detonable with relative ease (approaching that for high ex- 
plosives) if they contain a continuous matrix which is itself 
explosive. For a given class of high energy matrix, sensitivity 
tends to increase with energy. The search for more energetic 
propellants of that type will be complicated by this fact. 

Propellants with high energy matrixes presently seem to 
have a small specific impulse advantage over those with low 
energy binders. Perhaps careful formulation studies with 
medium energy binders might partially bridge this specific 
impulse gap and retain the valuable insensitivity of the low 
energy binder types. 

Defects, which are of such a nature as to permit the flow of 
gaseous combustion products through the matrix, significantly 
increase the shock sensitivity. A shock, or rapid pressure rise 
such as may accompany an improper ignition, can convert an 
otherwise satisfactory brittle propellant into an easily deto- 
nable system. The possibility must be considered in evaluat- 
ing the viscoelastic behavior of interesting formulations. 
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Servo Control of a 


Variable Thrust Rocket 


J. J. RODDEN! 
and R. J. POLLAK? 


Lockheed Aircraft Corp. 
Sunnyvale, Calif. 


Completed tests of a variable thrust engine operated in a closed loop formed by feeding back a 
voltage proportional to chamber pressure to an electrohydraulically actuated injector are described. 


In effect, a 


“thrust”? servomechanism was created which produces thrust proportional to a voltage 


input. These tests were conducted at the Naval Ordnance Test Station, China Lake, Calif. 
They demonstrated closed loop thrust response to a voltage command equivalent to a second- 


order system flat to about 6 cps. 


It was determined that through proper mechanical design 


a continuously controlled engine with full-range-throttlability could be produced with a bandpass 
limited only by the chemical reaction time for fuel ignition. 


N INCREASING need is apparent today for variable 

thrust liquid rocket engines. Requirements exist for 
controllable power to propel aircraft, spacecraft and missiles 
and for reaction controls in all of these (1).4 Engines with 
shutoff and restart capability and with two levels of non- 
zero thrust are in operation at present. The continuously 
variable thrust engine is being developed by many agencies 
(2,3). 

A variation in thrust is achievable practically by varying: 
Propellant flow rate, injector pressure differential, and nozzle 
throat area either singly or in combinations (4). When 
such an engine is operated in an unmanned vehicle or in an 
automatic system in a manned vehicle, some form of feed- 
back signal is mandatory to indicate when the thrust is at its 
commanded level. This signal, compared to the command, 
generates the error signal characteristic of a servomechanism 
and may be electrical, mechanical or pressure in nature. 
Convenient feedback parameters (other than thrust itself) 
are combustion chamber pressure, fuel flow rate and the 
position of a moving part, such as a flow regulator or injector. 
Reference (5) describes a design of a servo controlled bipro- 
pellant rocket utilizing two feedbacks, chamber pressure 
and mixture ratio. This motor successfully demonstrated 
the ability to accurately maintain the pre-set commanded 
levels of these parameters. 

Over the past two years, NOTS at China Lake, Calif., has 
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developed a variable thrust liquid rocket engine utilizing 
UDMH and RFNA (5,6). It employs only one moving 
part, a variable area injector, and has a demonstrated capa- 
bility to operate at any level from zero to rated thrust 
upon demand. Its basic simplicity renders it particularly 
applicable to storable systems. 


Description of Injector Operation 


A simplified schematic of the variable thrust injector is 
shown in Fig. 1. The flow of both the fuel and oxidizer, 
maintained essentially at a constant ratio, is controlled 
within the injector assembly by the movement of a single 
part, a cylindrically chaped pintle. The fuel flows along the 
inner surface of this pintle, and the oxidizer around the outer 
surface. The oxidizer to fuel ratio is established by the 
relative diameters of the two annular passages and the angles 
of the conical surfaces at the orifices. The basic elements 
of this design are the following: 

1 Hydraulic pressure input which acts against the pintle. 

2 The pintle, whose movement permits the mixing of 
fuel and oxidizer and meters the flow rate of the mixture. 

3 The spring which opposes the pintle movement. 

4 The combustion chamber and nozzle. 

The 100 per cent throttlable characteristic of the motor is 
attained by utilization of the hypergolic fuel UDMH with 
RFNA oxidizer which ignites on mixing. The sliding pintle 
regulates the variable orifice through which the fuel and 
oxidizer flow. 

Chamber pressure is a conveniently available indication 
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of engine thrust and was used as the feedback in the servo. 
A block diagram of the servo configuration tested is shown 
in Fig. 2. It consists of the injector and combustion cham- 
ber, hydraulic pressure control servo valve and pressure 
supply, an operational d-c analog amplifier, and a pressure 
transducer. 

The NOTS injector used in these tests has a rated fuel flow 
rate corresponding to 3500 lb generated thrust. For the 
purpose of these tests, the injector pintle travel was restricted 
mechanically to limit the thrust level to 400 lb, and the in- 
jector was used with a 400-lb combustion chamber. 

Operating at this reduced thrust level conserved fuel and 
thus allowed more firings per test stand fueling. The results 
gained from operating at this low thrust level are applicable 
to larger rocket motors, since the operation of the servo feed- 
back depends on the variation of the chamber pressure and 
not the generated thrust. Thrust is proportional to the 
product of chamber pressure, the thrust coefficient and nozzle 
throat area, or 


T = CrA; Pe 


The coefficient and area are dependent on the combustion 
chamber design. With different motor chambers different 
thrust can be attained even with no change in maximum 
chamber pressure. 


Results of Injector Tests 


The open loop tests consisted of determining the transfer 
function of chamber pressure to hydraulic control pressure 
inputs. The response was determined for ramp, step and 
sinusoidal inputs. For ramp inputs there was a high sensi- 
tivity of the chamber pressure and thrust to changes in the 
input current and hydraulic pressure. This open loop 
touchiness severely hampered testing, since a very small 
change in input current would swing the injector full open 
or closed. From these tests an approximate injector transfer 
function was estimated as a low frequency (about 1 cps) 
first-order lag. Further refinement of the injector function 
data was later found from the closed loop runs. 

Prior to closing the loop, tests of the pressure control valve 
were taken to determine the hydraulic system’s influence on 
closed loop operation. The frequency response of the valve 
indicated that the hydraulic pressure supply system would 
limit performance, since the pump had the characteristic of 
a 6-cps second-order system. It was unfortunate that the 
pressure supply limited the closed loop response, but this 
did not detract from the objective of the servo demonstration. 

After determining estimates of the injector and hydraulic 
transfer functions, an amplifier gain for closing the loop was 
estimated for stable closed loop runs. 

Comparison of closed loop response to ramp input with 
the open loop data showed improved sensitivity and linearity 
achieved with the stabilized system. The feedback also 
reduced burning roughness. 

Response to a sinusoidal input is shown in Fig. 3. The 
limitation in frequency response results from the pump 
dynamics. The frequency response of the closed loop in- 
jector to sinusoidal inputs was taken about an operating 
thrust level of about 200 lb. The experimental data from 
those runs for amplitude response of chamber pressure, hy- 
draulic pressure and error current are shown in Fig. 4. 

These responses were combined to find the hydraulic and 
injector transfer functions. By curve-fitting the latter 
functions with templates, idealized linear transfer functions 
were found. The results indicated the injector to be a first- 
order lag with a break frequency at 1.6 cps. With these 
linear transfer functions, a calculated closed loop response 
was determined which is shown with the experimental points. 
It is felt that reasonable agreement was found. 
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Fig. 2 Block diagram of servo configuration 
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Fig. 3 Closed loop response to sinusoidal input 


Analysis of Injector 


A theoretical analysis was performed to obtain a mathemati- 
cal representation of the injector. The validity of this model 
is confirmed by the experimental results. The model is to 
serve for greater understanding and as an aid for design. 
The block diagram of Fig. 5 is shown in order to correlate 
the experimental results with the mechanical design of the 
injector. The net force resulting from the hydraulic pres- 
sure, damping, spring force and chamber pressure acts to move 
the pintle. There are basically two kinds of damping: 
Stiction or coulomb and viscous. The coulomb friction 
adds or subtracts an offset force to the hydraulic input de- 
pendent on the direction of pintle travel. This effects a 
hysteresis in the closed loop operation. The force from the 
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chamber pressure is regenerative, since an increase in cham- 
ber pressure forces the pintle open further to increase flow 
rate and thus produce more thrust. The pintle spring tends 
to balance this chamber force. The chemical reaction time 
between the pintle displacement and the formation of the 
thrust and chamber pressure was found to be small (of the 
order of 50 millisec) compared to the dynamic time constants 
of the rest of the loop, and can be neglected here. This reac- 
tion time can be further reduced with the use of splash plates 
to aid fuel oxidizer mixing. 

With the assumption that the ratio of pintle mass to vis- 
cous damping M/B is small, the injector transfer function can 
be approximated to 


Po_ An Ke/B 
Py S + (K-Ke Ac)/B 


| NG 
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Fig. 4 Closed loop frequency response 
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Fig. 5 Variable thrust injector block diagram 


where 


Pe = chamber pressure 

Py» = input hydraulic pressure 

Ac = effective chamber pressure 

Ay = hydraulic pressure active area 

Ke = slope of chamber pressure vs. pintle displacement 


function 
B= viscous damping coefficient 
K = pintle spring constant 


This first-order system agrees with the form of the experi- 
mentally found transfer function. 

Some of the injector parameters shown were determined 
readily. The measurements of the hydraulic and chaml.er 
pressure areas Ay and Ac, were taken from an assem!)ly 
drawing of the injector. The effective chamber pressiire 
area Ac acting on the pintle had to be estimated, since ihe 
fluid fuel and oxidizer have influence in modifying the effec- 
tive area. The pintle spring constant K was measured 
experimentally on the “cold” injector. 

A numerical check of the gain and break frequency of ihe 
injector frequency response cannot be made directly, since 
there was no measured data for damping coefficient B and the 
chamber pressure vs. pintle displacement slope Ke. 

In an effort to estimate the injector gain Ke, the approxi- 
mate chamber pressure vs. pintle position function was 
constructed from oxidizer and fuel flow rate data. Examina- 
tion of the slope of this function showed that Ke varies over 
a range of 10.1 (from 17,000 to less than 1400 psi per in.). 
From the injector transfer function it can be seen that this 


-nonlinearity changes the injector’s first-order break fre- 


quency and gain. 

An estimate of the viscous damping coefficient was found 
with the frequency response data. From the injector trans- 
fer function the steady-state gain and break frequency can be 
found. These expressions can be solved for the damping 
coefficient as a function of break frequency and gain 
(Ac/AnG + 1) 


Substituting the measured and estimated parameter values of: 


Damping coefficient B = 


Spring constant K = 11,0001b/in. 
Hydraulic pressure area Ay = 3.68 in.? 
Chamber pressure area Ac = 1.82 in.? 
Gain (from freq. resp.) G = 5.6 psi/psi 


1.6 eps — 10 rad/sec 


Break freq. (from freq. We 
resp.) 
into the foregoing expressions gives an approximate damping 
coefficient of 
B = 366 lb-sec/ft 


This viscous friction force is distributed between O-rings 
and fuel, oxidizer and hydraulic fluid film surfaces. At 1.6 
cps, the peak viscous force on the pintle is of the order of 10 
to 20 lb at maximum amplitude. 

The characteristic root of the injector first-order lag 
changes as the injector gain Ke varies. At the critical value 
of Ke = K/Ac the injector root crosses to the unstable region. 
For the injector tested at NOTS the critical value was about 
Ke = 8000 psi per in. which occurred at about half full 
thrust. The shifting of the injector root is not readily de- 
tected in the frequency response, since the average Ke value 
of the thrust excursion determines the basic response. ‘The 
effect of this root’s presence near the origin was noted in the 
open loop response to a ramp input in that the chamber 
pressure and thrust went from 100 per cent to zero with prac- 
tically no change in hydraulic pressure input showing the 
“integrator” effect of the injector transfer function. ‘lhe 
presence of an unstable root was evident during open loop 
tests in the difficulty of trying to operate between saturated 
and zero thrust. 
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In the NOTS test a control loop was formed by connecting 
the output of the chamber pressure transducer to the servo 
amplifier which regulates the hydraulic servo valve. The 
dynamics of this loop is strongly influenced by pintle gain. 
The range of closed loop for the gains and parameters of the 
NOTS test if pump dynamic effects could be neglected would 
be from 6 to 10 cps. Unfortunately, the pressure supply 
pump used in these tests limited the overall response. 


Design Recommendations 


The foregoing analysis is useful in determining considera- 
tions for future design of servo controlled injectors. In a 
high performance servo the break frequency is as large as 
possible. Examination of the injector closed loop break fre- 
quency expression can indicate how to improve performance. 
The closed loop break frequency is 


K KcAu > Ac 
We B + B 


= amplifier gain, v/psi chamber pressure 
Ky = servo gain, psi hydraulic pressure per volt 


| 


It is to be noted that this expression is valid only when the 
bandpass of the other elements of the closed loop, such as 
the hydraulics and electronics, is large compared to the 
overall closed loop break frequency. 

The first term of the frequency of both the open and closed 
loop injector is the ratio of spring constant K to viscous damp- 
ing coefficient B. This term can be increased by stiffening 
the spring and minimizing damping. The maximum spring 
rate is limited by the amount of force available from the con- 
trol pressure source to deflect the pintle open. The damping 
can be reduced by limiting the sliding surface area to the 
pintle and minimizing O-ring seals. 

To reduce the effect of chamber pressure feedback which 
lowers the natural frequency, the ratio of net area of the pintle 
exposed to the chamber pressure Ac to the net active hy- 
draulie pressure area Ay should be made a minimum. 

The feedback gain consisting of the product of the gains 
of the chamber pressure transducer, electronic and hydraulic 
systems KpKy should be made large. The limit of these 
gains must be determined from considerations of the pintle 
gain Ke and any dynamic lags in the system for stability. 
Noise considerations are also influential in determining limits 
of the feedback gain. 

Further need to eliminate O-ring seals exists to reduce 
hysteresis in the closed loop system which is due to stiction 
forces on the sliding pintle. This hysteresis effect can be 
partially compensated, at least at low frequency, by incor- 
porating an integral-lead network in the servo loop. 

The injector gain Ke is probably the most difficult param- 
eter to design, since the chamber pressure vs. pintle opening 
will be fundamentally quite nonlinear. Attempts should be 
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made to make the slope of the flow characteristic near- 
linear and the magnitude of the slope Ke as large as is feasible, 
while maintaining the mixing ratio and flow rates required 
for combustion. This consideration might be achieved by 
minimizing the total travel of the pintle from open to close. 

The injector design discussed here requires control of the 
absolute value of input hydraulic pressure to provide actua- 
tion force. This generally requires the hydraulic supply 
system to circulate large flow rates of fluid—which involves 
considerable power consumption. It was this taxing of 
the hydraulic pump that limited the performance in these 
tests of the NOTS injector. A pintle design that can be 
actuated by a differential in hydraulic pressure would elimi- 
nate any large pump flow requirement. In addition such a 
design would allow elimination of the restraining spring, 
since the feedback loop would regulate the actuation force to 
position the pintle. 

Commonly available high performance servo valves which 
produce differential pressure or flow would be readily adapt- 
able to this alternate pintle design. 


Conclusions 


The feasibility of incorporating the NOTS variable in- 
jector into a servo controlled system was demonstrated. 
It was found that the servo control effectively linearized 
the injector nonlinear characteristic in addition to augment- 
ing stability. 

The analysis of the injector indicated a number of con- 
siderations applicable to design of injectors for fast response 
servo systems. 

Since the servo controlling of an injector of the NOTS type 
is not dependent on the thrust produced, it appears that 
the thrust capability of such a servo controlled injector is 
limited only by the size of fuel and oxidizer flow rates that 
‘an be satisfactorily mixed for combustion. 
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Technical Notes_ 


Effect of Different Types of Video 
Filters on PDM/FM Radio Telemetry’ 


A. T. BUBLITZ? and M. H. NICHOLS? 
University of Michigan, Ann Arbor, Mich. 


GLOW (1)* has compared analytical results of PDM/FM 
telemetry based on two different video filter character- 
istics. The results compared by Uglow were taken from Feld- 
man and Bennett’s analysis (2) with a Gaussian video filter 
(the model used by Feldman and Bennett for PPM/FM 
applies equally well to PDM/FM) and from Nichols and 
Rauch (3) with an ideal sharp cutoff low pass filter hereafter 
abbreviated by ILPF. This paper corrects the method of 
optimization used in chapter 11 of the book by Nichols and 
Rauch, and extends the analysis of Feldman and Bennett. 


Comparison of Filters 


Both the Gaussian and the ILPF characteristics are im- 
possible to realize. In practice, the ILPF is usually approxi- 


proportional to the ratio of the slope of the step response of 
the video filter at the slicing level to the rms noise voltage 
at the output of the filter (which includes the second and 
third of the three criteria for low pass filters). 

The response of an ideal vertical cutoff low pass filter |ias 
an initial overshoot of 9 per cent to a unit step input, and the 
oscillations die off approximately inversely proportional 
to time (8), a slope of 2.0 F. at the slicing level (slicing level 
taken equal to half height), where F. is the cutoff frequency 
and the equivalent noise bandwidth is F.. A Gaussian filter 
has no overshoot in its response to a unit step input; it /ias 
an effective noise bandwidth of 1.34 F, with an FM video noise 
spectrum, where F,. is the 3-db frequency of the Gaussian 
filter, and the slope of the step response is 3.0 F, at the slicing 
level. The ratio of the maximum slope of the step response 
to the three halves power of the equivalent noise bandwidth 
of an ILPF (used in a PDM/FM system) is 2.0 F.~'”2, and for 
the Gaussian filter it is 1.94 F.-’/* > 2.0 F.-’”. As shown 
later in a numerical illustration, the use of the Gaussian filter 
in practice eliminates the problem of cross-talk in PDM/FM. 


Output Signal to Noise Ratio of a PDM/FM System 


The rms output signal to noise ratio S/N of a PDM/FM 
system when well above threshold is, approximately‘ 


N 


mated by a maximally flat gain or Butterworth type of filter 
characteristic; the Gaussian low pass filter is most often 
approximated by a maximally flat phase or Thomson type 
of filter characteristic. See, e.g., reference list in (4). From 
the point of view of using such a filter at the output of the 
demodulator in a PDM/FM system, three characteristics of 
the filter performance are of interest: 

1 The rate of decay of the transient in the response to a 
step input, which determines the amount of cross-talk be- 
tween adjacent pulses. 

2 The rate of rise (slope) of the response to a step input 
at the slicing level. 

3 The equivalent noise bandwidth of the filter (using 
as basis of comparison an ILPF). 

When a low pass filter is used at the output of an FM de- 
modulator, the portion of the noise power spectrum inter- 
cepted by the filter will be proportional to the cube of the 
bandwidth of the filter; i.e., the rms video noise voltage at the 
output of the filter will be proportional to the three halves 
power of the filter bandwidth. For any low pass filter other 
than the ILPF, a definition of bandwidth F. has to be chosen— 
such as the 3-db point of the amplitude response curve which 
is used in this paper. The rms output signal to noise ratio 
of a PDM/FM system using a particular low pass filter is 
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D. C., Nov. 16-20, 1959. 
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Air Development Center. 

* College of Engineering. 

* Numbers in parentheses indicate References at end of paper. 


S _ (rms full modulation time) (pulse slope at slicing level) ( F ) Ve il] 


rms fluctuation noise voltage in video fm 
where 
F = number of samples per channel per sec 
fm = cutoff frequency of the low pass output filter in each 
channel 
pulse slope at slicing levelé _ 4 V3 kSF fp 
rms fluctuation noise voltage in video 
where 
S = rms value of signal in IF, in volts (referred to output) 
ks = constant defined so that the slope of the unit step 
response of the filter is 2 k;F. at the slicing level 
F, = 3-db point of filter, cps 
k,; = rms fluctuation noise voltage per root unit band- 
width in IF channel (referred to output) 
By = equivalent noise bandwidth of video filter for FM 
= keF., where ke is a constant 
fo = peak frequency swing of FM link when fully modu- 
lated 
full modulation time® = 
rms full modulation time® = T 
where 
nm = number of channels 


‘ See, for example (3), - 60. 
. oe level assumed to occur at half steady-state pulse 
eight. 
® The factor 1/2*/: appears, since this derivation assumes sinus- 
oidal modulation; thus the rms value is 1/2*/: of the peak-to- 
peak swing. ; 


Eprror’s Note: The Technical Notes and Technical Comments sections of ARS Journal are open to short manuscripts describing 
new developments or offering comments on papers previously published. Such manuscripts are usually published without editorial 
review within a few months of the date of receipt. Requirements as to style are the same as for regular contributions (see masthead page). 
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7 = guard time = spacing between pulses to reduce cross- 
talk plus minimum pulse length 


The factor (F/2f,)'/* in Equation [1] takes into account 
the frequency components of the noise lying between F/2 and 
f, taken out by the individual channel low pass output filter 
which is assumed to cut off sharply at f,,. 

By substitution into Equation [1], the rms individual 
channel output signal to noise ratio at carrier threshold (S/N), 
is obtained 


where 
By) = IF channel noise bandwidth 
p = constant = Bo/2fp 
k; ks/ke’/? 
The constant 6 is the FM carrier threshold factor given by’ 


S; = 6k, Bo /? [2b] 


when S; is the rms value of the carrier at threshold. 

The constant p can be estimated from the calculations of 
Salinger (6) on the response of an ideal vertical cutoff band 
pass filter to a step change in frequency. In a PDM/FM 
system, the carrier frequency fo is modulated by shifting it 
from fo — fo to fo + fo which, according to the results of 
Salinger, gives a pulse at the output of the discriminator ap- 
proximately of the form Si(Bot) plus a damped ringing os- 
cillation equal in frequency to 2fp. The apparent damping 
as well as the magnitude of this ringing is dependent on the 
value of the ratio Bo/2fn = p, and as the value of p approaches 
1.0, the magnitude of the ringing approaches a maximum 
and the damping approaches zero. For deviation ratios 
fo/l. = D greater than 0.5, this ringing oscillation will be 
attenuated to an amount depending on the actual character- 
istics of the IF bandpass and of the video low pass filter. 
Curves giving the time response of the carrier frequency to a 
sudden shift are given in Salinger’s paper and also in Appendix 
13 of (3). From these curves it appears as though p = 1.4 
would be a conservative choice. If the deviation ratio is 
greater than 0.5 and the action of the video low pass filter 
is considered, then the value of p can be allowed to approach 
closer to 1.0. [See (3), p. 65.] 


Output Signal to Noise Ratio of a PDM/FM System 
With a Gaussian Video Filter 


The guard time 7 depends on the characteristic of the video 
filter and must be chosen to satisfy cross-talk criteria and the 
minimum allowable pulse height. The minimum allowable 
pulse height determines the minimum pulse length, and the 
minimum pulse gap between the end of a fully modulated 


’ This assumes for purposes of analysis that a well-defined 
threshold exists for all deviation ratios. Some idea of the 
validity of this assumption can be obtained by reference to the 
results published by Stumpers (5). 


Fig. 1 Schematic showing reduction of pulse height 
OcroBer 1960 


pulse and the beginning of the next pulse. Since the step 
response of a Gaussian filter has no overshoot, it turns out that 
7 is controlled by the minimum pulse height requirement. 
Fig. 1 illustrates schematically the pulse height effect for a 
particular spacing with a Gaussian filter where 7 = 7 + 
72. If 7, = 72, then hi; = he. The values of h; and he deter- 
mine the leeway for noise discrimination at minimum modula- 
tion and full modulation, respectively. For, if the instan- 
taneous video noise voltage exceeds he, the pulse plus noise 
will exceed the slicing level for the long pulse, and if the in- 
stantaneous video noise voltage is less than —/, the short 
pulse plus noise will not reach the slicing level. Thus h, 
and h, determine the video threshold for minimum and maxi- 
mum modulation.’ Also, Fig. 1 shows that the slopes at the 
slicing levels are reduced, which increases the output noise 
at full and minimum modulation when above threshold.° 
These considerations determine the minimum permissible 
values of 7; and 72. The quantities 7; and 72 so determined 
can be written as 


wie 


Ti T2 
where k; and ky are constants. Equation [2a] then becomes 
for this case 


Following the procedure of Feldman and Bennett (2), 
(S/N), can be maximized by varying F., holding the other 
variables constant, i.e., for constant IF bandwidth By and 
constant nF. This gives 


F./nF = 3(ks + ks) [4] 
Substitution of Equation [4] into Equation [3] gives 


§ 
N max 3F. nF Pp 


The required IF bandwidth Bo to give (S/N); max at carrier 
threshold is obtained from Equation [5] 


Bo _ [1.22 (p\(2im\ (%)" 
nF ky (?)( nF max (6a 


For purposes of illustration, we will assume that 7 = 72, 
which makes h, = he, and that the threshold requirement is 
that h, = h, = h/4. The total guard time is r = 7 + 72. 
Writing the transfer function of the Gaussian filter as exp 
(—c?f?/F.2), where F, is the frequency for 3-db attenuation 
and c? = 0.35, the guard time + for this case is computed by 
use of error function tables to be 0.61/F. sec. From Equa- 
tion [4], F./nF = 1.8, and from the discussion of the com- 
parison of filters, k; = 1.0, so Equation [5] becomes 


It is now necessary to check to see that cross-talk is at a toler- 
able level under the modulation condition illustrated in Fig. 1. 
This is done by dividing the value of the “tailoff’”’ of the long 
pulse at the time the trailing edge of the short pulse crossed 
the slicing level by the slope of the short pulse at that time. 


8 Assuming that the leading edges of the pulses are fixed, the 
modulation sensitivity near full modulation will decrease _be- 
cause of the presence of the leading edge of the next pulse. This 
effect can be calibrated out. 

9 It is conceivable that in some cases the reduction of the slope 
at the slicing level is of more concern than the reduction of the 
pulse height, so that a minimum allowable slope becomes the 
criterion for minimum allowable pulse width. The general pro- 
cedure is the same in either case. 
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To obtain the cross-talk relative to full modulation, we divide 
by the full modulation time, that is, by (1/nF — 7). For 
this illustration, the relative cross-talk is approximately 
0.4 X 10~, which is certainly adequate! It is also necessary 
to check to see that the video threshold is satisfied when the 
carrier is above threshold, i.e., for Equation [6b] to be 
valid, S, > S,, where S; is the rms value of the carrier required 
for carrier threshold and S, is the rms value of the carrier re- 
quired for video threshold. Using the principles for video 
threshold from (3), p. 66, and the expression for carrier thresh- 
old from Equation [2b] and substituting for Bo from Equation 
[6a] and F. from Equation [4], using the relation 2pfp = Bo, 
the ratio of the two threshold values of S becomes 


11.3ks(ks + kas) 
S: (2fm/F)'/(S/N)t max 
Using this numerical illustration, k; = 1.5 and ks + ky = 0.61; 


also (2f,,/F)'/? is the order of unity, and we choose (S/N); max 
= 100. With these values 


Thus for reasonable values of the variables, the video thresh- 
old is satisfied for a lower signal strength than the carrier 
threshold and the results are valid. 


Output Signal to Noise Ratio of a PDM/FM System 
With an ILPF 


In the case of the Gaussian filter, the consideration of mini- 
mum allowable pulse height determines the guard time 7 re- 
quired. Since the envelope of the step function response of 
an ILPF dies off approximately inversely with the time, it 
turns out that cross-talk reduction rather than minimum 
pulse height determines the guard time 7 required. Following 
the procedure in (3), p. 245, the maximum relative cross-talk 
7 is given by 


y = 0.023 nF/F27(1 — nF [7] 
Solving for + and choosing the root which gives 7 > 0 as 
1 2 
2nF 2nF 


Following the same procedure as in the previous section, 
this value of 7 is substituted into Equation [2a] and (S/N), 
maximized against F., holding nF and By constant, giving 


0.041 


This is the same result as obtained in (3), p. 105. 
Substitution of these values and k; = 1 into Equation 


[2a] gives 


Thus if 60-db cross-talk suppression (relative to full modu- 
lation)" is required under the modulation conditions of 
Fig. 1, the factor 1.6y'/* becomes 0.3, which may be compared 
with 0.6 in Equation [6b]. The required IF bandwidth 2, 
to give (S/N); max at carrier threshold is obtained from 


Equation [10] 


By _ 0.73 [? (*) 
nF F 6 N/t max 


It is necessary to check to see that the video threshold is satis- 
fied when the carrier is above threshold. Using the sare 
procedure as with the Gaussian filter, one obtains 


1.9 
S, (*) (#-)" [12] 
N max F 


100 and (2f,/F)'/? ~1, this 


Using y = 10-3, (S/N): max 
relation gives 


= 0.6 


This shows that for these values, the video threshold is satis- 
fied. Equation [10] is valid only for values of (S/N)¢ max > 40 
if y = 10-4, ete. For values such that S, < S,, the value of 
S/k, obtained from the video threshold expression must be 
used instead of Equation [2b] and the optimization carricd 


forward. 


Conclusion 

Based on the special numerical examples, the output signal 
to noise ratio for the Gaussian type of video filter is approxi- 
mately twice that for the ILPF. Also, in practice, it is pos- 
sible to approximate very closely the Gaussian or maximally 
flat phase behavior (4). For example, a tenth-order low 
pass Bessel filter has a step function response very close to 
that of a Gaussian filter with only 0.12 per cent overshoot. 
A tenth-order Butterworth filter has nearly 20 per cent over- 
shoot, or twice the overshoot of the step function response of 
an ILPF. 
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10 Tt might be argued that 60 db is overly stringent, since 
cross-talk of this sort occurs only when one channel is near 
full modulation and the following channel is near minimum modu- 
lation. Questions of this sort can be settled only by detailed 
analysis of a particular measurement program. 
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Gyro Time Constants 


ALLEN LANGE! 
Bendix Corp., Systems Div., Ann Arbor, Mich. 


N (1),? page 89, a discrepancy between the computed and 
measured values of an integrating rate gyro time constant 
is noted. It is suggested in (1) that the elastic compliance 
of the gimbal creates a time delay which may contribute to 
this observed discrepancy. The purpose of this note is to 
suggest an additional reason for this discrepancy, namely, 
the use of a familiar approximation for the angular momen- 
tum of the gyro rotor. Two expressions for the gyro transfer 
function are developed, one using the familiar approximation 
for the angular momentum and one using a more exact rela- 
tion. It is shown that the gyro time constants associated 
with these two transfer functions are different, and an attempt 
is made to estimate the numerical difference in the two 
time constants. 


Angular Momentum of the Gyro Rotor® 


let a gyro rotor be contained in a single gimbal which is 
mounted on a platform rotating with respect to a Newtonian 


frame. The application of Newton’s second law to the spin- 
ning rotor may be expressed in vector form as follows 
M, Di(J,W ir) [la] 
where 
JJ, = moment applied to the gyro roto: 
D{ ) = vector operator indicating the time rate of change 
of the vector with respect to a Newtonian frame 
J0 0 | 
he = inertia tensor of the gyro rotor = |0 J 0 
00 Jsp | 
Vie = angular velocity of the gyro rotor with respect 
to the Newtonian frame 
Equation [la] may be written as 


where H, = J,W;, is the angular momentum vector of the 
gyro rotor. 

J, is invariant only with respect to a set of coordinates 
rotating with the rotor. Consequently, Equations [la] 
or [2a] may be expressed as 


M, = + Wie X (J-Wiz) (1b] 
MW, = DA(H,) + Wi, X A, [2b] 


D,{_) is a vector operator indicating the time rate of change 
of the vector with respect to the rotor coordinates r. Now, 
H, may be expressed as follows, in terms of components in a 
coordinate set g attached to the gimbal on which the gyro 
rotor is mounted 


J ( Wit =~) 


Ai,’ = Hsp (=) [8a] 
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where 
Hsp = Jspw 
w = rotational speed of the gyro rotor. with 
respect to the gimbal 
DG = rotational speed of the gimbal with respect 


to the platform on which it is mounted 

W,, W2, W3 = x, y and z components, respectively, of the 
rotational speed of the platform with 
respect to the Newtonian frame 2, in 
gimbal coordinates 


The superscript g indicates that the vectors are to be ex- 
pressed in terms of their components in gimbal coordinates. 
Similarly, the superscripts r and p will be used to indicate 
rotor and platform coordinates, respectively. 

In a practical gyro, w is a large number, so that 


W3/w< 1 J/Hsp 0 


Then Equation [3a] may be written as 


0 
0 
1 


fi? ~ Hsp [3b] 


Consequently, it is common practice to consider that the 
term D,(H,) in Equation [2b] is zero. The assumptions 
leading to the approximate Equation [3b] are physically 
reasonable; however, it is demonstrated in the following that 
the use of Equation [3b] leads to an expression for the gyro 
reaction torque different from that resulting from the use of 
Equation [3a]. The two expressions resulting from the use 
of Equations [3a and 3b] are developed and compared in the 
following section. 


Gyro Rotor Kinetics 
Combining Equations [2b and 3b] leads to the following 
result 


Wit DG) 
= |M2| = Hsp W, x = 
| Ws 


2 
Hsp |—(Wi + DG@)| [4] 
0 


The first component of Equation [1b] may be expanded as 
follows 
DAT ie”) J,D,(W i?) 
But 
D, Wi? D; Wi? + x Wi? = D; Wi? 
Now, Wi, may be written as the vector sum Wi, = Wip + 
Woo + Wor so that 


= iv?) + Do( Wor’) + Wie? X Wor’ [5] 


where 7',, is a scalar operator transforming I;,? into W;,’. 
For purposes of the present analysis, let the gyro output axis 
be aligned with the z axis of the platform. Then 


1 0 0 
0 cosG sinG 


To = 
0 —sin G cos G 


and the first term on the right-hand side of Equation [5] be- 
comes 


Di(T ip?) D(T Wis? + T opD W in”) = 


Ay 
Az:+ [6] 
A; — DG W:2 
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d 
| 
| 
| 
M = D,(H,) [2a] 
| 
: 
Ws | 
| 
| 


where A;, Az and A; represent the z, y and z components, 
respectively, of the rotational acceleration of the platform 
with respect to the Newtonian frame 7 as expressed in gimbal 
coordinates. 

The second and third terms on the right-hand side of 
Equation [5] may be expressed as 


+ 
DW") + X We? = |WsDG — (Wi + DG)w 
'—W2DG 


Note that the synchronous motor driving the gyro rotor has 
the characteristic that Dw = 0. Then the first component 
of Equation [1b] can be written by combining Equations 
[6 and 7] to give 


[7] 


J(Ai + D*G + Ww) 
+ 2DG W; DG) w) 
J sp(As — 2 DG 


DJ [8] 


The second component of Equation [1b] may be written 


(Jsep — J)(Ws + w) We 


X Je = | sre — DG)(Ws + | [9] 


so that, combining Equations [8 and 9], the exact expression 
for D,(H,”) is given by 


M, | + D*G) + Hsp We + 
(J sp J)W.W3 

M,| = |J(A2 + 2DG Ws) + [10] 

(J — Jsp)W3(Wi + DG) + Hsp(W, + DG) 


| 
| 
We are concerned here with only the x component of Equa- 
tions [4 and 10]. Denoting M, from Equation [4] by M. 
(approximate), and M, from Equation [10] by M, (exact), 
we have 


= + D*G) HspW, + (J sp J)W.W; {11] 
M. HspW, [12] 


Equations [11 and 12] represent the effect of a torque about 
the gimbal axis on a spinning rotor treated as a free body. 
It follows that the reactive torque on the gimbal about its 
xz axis due to the spinning rotor is given by —M, or —M.. 
In the following section, the gimbal is treated as a free body, 
and its motion as the result of the net applied torque is com- 
puted. 


Jsp(A3 — 2DG W2) 


Gyro Gimbal Kinetics 


Newton’s second law applied to the gyro gimbal may be 
expressed as 


M., 
M ve 
M. 


where J, is the inertial tensor of the gimbal, which may have 
the form 


=> DJ ig”) [13] 


We are interested here only in the x component; M,, is given 
by 


M., = M.—- CDG +M, [14] 
where 
M. = applied control torque 
C DG = damping torque due to the product of C the 


damping factor and the relative angular ve- 
locity between the platform and the gimbal DG 


M, = gyroscopic reaction torque, where M, = —M, or 
—M.,, corresponding to the use of the exact 
expression (Eq. [11]) or the approximate 
expression (Eq. [12]) for the gyro reaction 
torque 


The right-hand side of Equation [13] can be examined in 
the same fashion as the right-hand side of Equation [1]. 
Performing this expansion, and equating the xz component to 
Equation [14] gives 


M. — CDG + M, = J.(Ai + D°G) + (WJ. — Jy)WeWs = [15] 


Rearranging, and denoting the terms resulting from the 
use of Equation [11] by the subscript e, and those resulting 
from the use of Equation [12] by the subscript a, we have 


1 M. HspW2 
wr. +1)¢= 5 
(1.4: + wavs) | [16] 
_ — HseWs 
(DT.+1)G = D 4 
(1.4: wars) | [17] 
where 
IAC 
= J)/C 


Equation [16] is identical to Equation [10] of Derivation 
Summary 2 of (1) with the following exceptions: 

1 Equation [16] does not contain a torque due to un- 
certainty. 

2 Equation [16] contains the term [(J, — J2z)/C]W2Ws, 
due to inertial coupling of the gimbal, which disappears if 
J, = J., and in any case may be quite small if J, ~ J, and 
the product is small. 

However Equation [17] differs from Equation [16], and 
hence from the equivalent expression developed in (1), by 
the inclusion of the inertia of the gyro rotor about an axis 
normal to the spin axis with the gimbal inertia. Therefore, 
the ‘‘exact” gyro time constant T, is given by (Jz + J)/C, 
whereas the “approximate” gyro time constant 7. is given 
by J./C. (As before, the notation “exact” refers to the use 
of Eq. [11] for expressing gyro reaction torque, and the nota- 
tion “approximate” refers to the use of Eq. [12].) On page 
89 of (1), the observation is made that the measured gyro 
time constant is greater than the computed gyro time con- 
stant, which, presumably, is computed from an expression 
similar to Equation [16]. Since 7, is greater than 7, it is 
pertinent to examine the possibility that the measured time 
constant exceeds the computed time constant, because the 
computed time constant is based on the reasonable approxi- 
mation D,(H,) = 0. A quantitative evaluation of this 
hypothesis, based on some simplifying assumptions, is given 
in the following section. 


Rotor and Gimbal Inertias 


The ratio of the exact to the approximate gyro time constant 
may be written as 


J 
— 8 
[18] 


From (1), we know that the ratio of 7, the measured time 
constant to TJ, the computed value is greater than 1 or 


T/T, =1+A [19] 
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Using the values given on page 89 of (1) 
Tm = 0.0027 sec T. = 0.0017 sec 


The hypothesis is that 7’, the measured time constant is to 
be associated with 7’, the exact time constant, and 7’, the 
computed time constant is to be associated with 7, the ap- 
proximate time constant. Then A~J/J., and we seek a 
numerical value for J/Jz. The ratio J/J, may be written as 


J z J SP J z 
Since J/J, represents a degradation of the dynamic perform- 
ance of the gyro (because it increases the gyro time constant), 
it is desirable to keep J/J, small. Also, since a large Jsp 
corresponds to a highly sensitive gyro, and a small J, corre- 
sponds to a gyro with good dynamic performance, Jsp/Jz 
should be large. Therefore J/Jsp should be as small as 
possible. 

For a solid homogeneous spherical rotor, J/J sp is equal to 

1. Forasolid homogeneous cylindrical rotor 


A = 0.59 


where n = h/r, the ratio of the height of the cylinder to its 


radius. In the case of a hollow cylindrical rotor 


or 1 n? 


where m is the ratio of the inner to the outer radius, with 
<2 < I. 

For n = 1,m = 0, J/Jsp = 0.67. Therefore, it seems rea- 
sonable to assume that J/Jsp < 1.0. From the table on 
page 76 of (1) 

Jsp 11.97 | 
so that 
J/J, 037 


and it can be concluded, therefore, that the use of the approxi- 
mation given by Equation [3b] accounts for only half, at 
most, of the observed discrepancy A between the measured 
and computed gyro time constant, the remainder, presum- 
ably, due to the elastic effects suggested in (1). 
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Satellite Librations 
C. F. BERNINGER! 


Aveo Research and Advanced Development Div., 
Wilmington, Mass. 


HE PREVIOUSLY published work on satellite libra- 

tions assumes the mass center of the satellite to be 
traveling on a Keplerian orbit. The “librations’’ or oscilla- 
tions of the vehicle about its center of mass are due to differ- 
ences in the gravitational and centrifugal forces acting on the 
various mass elements of the satellite—the centrifugal force 
terms determined from the specified Keplerian orbit. Cou- 
pling between the orbital and libration motions, which appears 
obvious from considerations of conservation of momentum 
or energy, may be illustrated by writing the equations of 
motion in general form. 

For simplicity, the equations of motion are written for a 
dumbbell-shaped satellite consisting of two point masses 
connected by a weightless rigid rod of length 21. The radius 
between the center of mass of the dumbbell system and 
center of Earth is denoted as p (Fig. 1). The radii from the 
center of Earth to each mass are 7; and rz. The angular posi- 
tion to the center of mass of the satellite is given by 6 meas- 
ured from a fixed reference line 0A and the angular attitude 
of the satellite is specified by ¢. 

By employing Lagrange’s equations with p, @ and ¢ as the 
generalized coordinates one obtains 


2mp — 2mpés? = R 


d 
at (2mp*6) = 0 


ll 


d 
(2ml?*¢) 
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where R, 0 and @ are the generalized forces. 
The gravitational forces acting on the masses are designated 


as G, and G; and are given by 


m2 + + 26 cos 2) 
= 
+ — 26 cos Q) 
where 
6 = I/p 
Q= 
u = positive gravitational constant 


The gravitational force system G, and G, may be resolved 
into components acting at the center of mass of the satellite 
consisting of a radial force R, a force 7’ perpendicular to R, 
and a torque ® acting about the mass center. 

The sum of the radial components of G, and G» parallel to 


pis given by 
cos [sin~! (J sin 2/7) ] 
R=-—-—ypm 


2 


cos [sin— (J sin 2/rs) [3] 


The torque @ is readily determined by resolving G, and G2 
into components parallel and perpendicular to the longi- 
tudinal axis of the dumbbell (Fig. 2). The perpendicular 
components are 


F, = G, cos \ 

F, = G2 cos d2 

and, from trigonometric relationships 
cos \y = psin Q/n 


cos kz = p sin 2/re 


n 
J sp 6 
1 


CENTER 
: OF EARTH 


Fig. 1 Force resolution parallel and perpendicular to radius of 
satellite center of mass 


The torque about the dumbbell center of mass is then given 
by 
1 
© = umopl sin Q (4 +) [4] 


23 


The generalized force 6 corresponds to the torque about the 
center of Earth. Since G; and G; pass through 0 


[5] 


Substitution of Equations [3, 4 and 5] into Equations [1] 
yields 


1 
cos sin 6 
(6) 
(d/dt)(p*6) = 0 [7] 
d 
(2ml?¢) = umlp sin Q (4 =) [8] 


Expansion of the trigonometric expressions of Equation [6] 
and expansion of the sums and differences of the reciprocal 
powers of r; and rz permit Equations [6 and 8] to be written as 


— (u/2p?) (2 + + 8 cos? Q — 2sin?Q)] [9] 
65? u ml sin 2 cos Q 
p 


[10] 


d 
2 = 


| CENTER 

OF EARTH 

Fig. 2 Force resolution parailel and perpendicular to dumbbell 
axis 


wherein powers of 6 greater than 2 are discarded. 
Differentiation of the left side of Equation [10] enables 
one to rewrite the equation 


pé = sin cos {11} 


In this form a comparison of Equations [11 and 9] reveals 
that the left sides of both equations are “of the order p.” 
The right side of Equation [11] is proportional to 6, whereas 
the non-Keplerian term in Equation [9] is proportional to 
5°. Since 6 is, for reasonable Earth satellites, very small com- 
pared to unity, a reasonable approximation is given by neglect- 
ing the 6* terms. Equations [9 and 7] then become Kepler’s 
equations of motion. Equation [10] is an auxiliary equation 
giving the librations and is equivalent to the equation used in 
(1,2and 3).2. It should be noted, however, that (4), in effect, 
uses Equation [8] for the investigation of the librations. 
The increased accuracy afforded by the use of this equation 
is questionable, since it includes terms comparable in magni- 
tude to those neglected by the approximation of Equations 
[6 and 7]. 
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Dynamic Stabilization of a Small Body 
Between Two Equal Rotating Masses' 


HOWARD D. GREYBER? 


Aerosciences Laboratory, General Electric Co., 
Philadelphia, Pa. 


N INTERESTING and simple illustration of the systems 

concept of dynamic stabilization is the stability of a 
small body lying on the line of centers and equidistant from 
two equal rotating masses which interact by gravitation but 
are held at some arbitrary distance by a weightless bar. 
Elsewhere, we have applied the dynamic stabilization con- 
cept to problems of hydrodynamic and hydromagnetic insta- 
bilities (1). Other applications of dynamic stabilization 
are the upside down pendulum with an oscillating support 
aud strong focusing in charged particle accelerators. 

Let us consider uniformly rotating frames of reference, 
Sy and + Bie the origins of which coincide. The acceleration 
in unprimed system is related to that in the primed rotating 
system by (2) 


+ + wx (0X4) 


Consider two equal masses as in Fig. 1 with w as shown, 
and the x, y coordinate system rotating with the masses. 
Examining the forces on a small body m, close to the ‘ridge 
point” 0 in Fig. 1, we may write 


Fy = mj + 2mwx — mw*y [2a] 


also 


1 1 
Assuming x and y small, we have 
F, = G(MmA4y/R?) [2b] 
Similarly 
F, = mé — 2mwy — = [3] 
Gathering terms 
+ + y(—w? — 4GM/R*) 
— Qwy + 2(—w* + 2GM/R*) = 
Letting = and y = be* and S = GM/R?', we obtain 
ba? + a(2wa) + b(—w? — 48) = 0 [5a } 
aa? — b(Qwa) + a(—w? + 28) = 0 [5b] 


If these are considered as a set of linear equations in a and 
b, they may be solved by setting the determinant equal to zero 
and obtaining 


Il 


[4a] 
[4b] 


| 


§—w+ 
(1/2)V/4(w? — 8)? — 4(w? — 28)(w2 + 48) [6] 


We note that taking the negative sign means that a? is more 
in the stable region, and thus to find stable regions we con- 
sider only the plus sign case. 

Our criterion for stability is the usual mathematical defi- 
nition that stability is implied when the dependent variables 
(in this case the coordinates of the test particle) remain finite, 
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i.e., bounded, as time goes to infinity. In addition, here 

we claim stability only if both solutions remain finite for the 

parameters assumed. The reason for this is the belief that 

the physical system follows the unstable solution if one exists. 
Reducing Equation [6] we find 


a? = S — + — [7] 


We now wish to discover if the unstable equilibrium of the 
particle on the “ridge point” with no rotation can be made 
stable for some value of w. Three cases are to be 
distinguished. 


Case I 


w? is less than S, a? is greater than 0, and a@ is real. 
Hence, the system is unstable. 


Case II 


S is less than w?, and w? is less than 9/4 S. In this region 
a’ changes so that the system changes from an unstable regime 
to a stable one. To find the boundary between stable and 
unstable regions, we set a? = 0, and find w? = 2S, — 48. 
Since the negative value for w? is meaningless physically, 
w = V258 is the boundary between stable and unstable re- 
gions. Thus, stability is from 2S < w? < 2.25.8. Thus, we 
see that a small body m near point 0 can be stable against 
small perturbations if w? is properly chosen. 


Case III 


w? is greater than 9/4 S. Here one can show that w? is 
complex; one finds a real part for a so that we have instability 
throughout. 


One should be careful to distinguish our problem from the 
usual problem treated in books on celestial mechanics, such as 
Moulton’s (3). In the celestial mechanics problem, the rota- 
tional frequency and distance of the two heavy bodies is 
determined, i.e., #2 = S/4, so no weightless bar is required. 


Fig. 1 Relative positions of the small body m and the rotating 
masses 


W 
m 
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Case I shows that for w? = S§/4, the solution is unstable, 
which agrees with the result in celestial mechanics for the 
straight-line solutions for the centers of libration. 

It is amusing to note that the two equilateral triangle 
point solutions are stable if u is less than 0.0385. Since for 
the Earth-moon system, » = 0.0121, the equilateral triangle 
points for the Earth-moon system are stable (4). If one 
ignores the many perturbing effects, one may anticipate finding 
a relatively high density of interplanetary dust at these equi- 
lateral triangle points for the Earth-moon system. It appears 
worthwhile to design an experiment with a space vehicle to 
travel to and investigate the space in the region of these 
points. 


Conversations with W. A. Newcomb, who suggested this 
illustrative problem, are gratefully acknowledged. 
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Fifth Harmonic of Earth’s Gravitational 
Field 


HERMAN F. MICHIELSEN! 


Lockheed Aircraft Corp., 
Missiles and Space Div., Palo Alto, Calif. 


Available information on the changes in perigee height 
of 1958 Delta Two seem to indicate an effect opposite to 
that which would result from the third harmonic of the 
gravitational field, as it was determined from observations 
on 1958 Beta. This effect is tentatively explained by the 
assumption of a fifth harmonic of slightly lower magni- 
tude than the third. The available evidence is not con- 
clusive but indicates the need for further study. 


N REFERENCE (1),? O’Keefe et al. arrive at the con- 
clusion that certain long period oscillations of the orbit 

elements of Vanguard I, specifically of the eccentricity, clearly 
indicate the presence of an odd harmonic in Earth’s gravita- 
tional field. Assuming, tentatively, this to be the lowest 
possible harmonic, that is, the third one, a specific coefficient 
for this anomaly was derived. Higher odd harmonics 
were surmised to be considerably lower in magnitude and 
insignificant. 

In this note, it is proposed to consider a fifth harmonic on 
top of the third. The words “proposed” and ‘consider’ 
are used here intentionally, since from the information avail- 
able to the author, he is not able to prove conclusively the 
existence of a sizable fifth harmonic. But the indications 
are considered sufficiently strong to warrant dissemination of 
the findings. In addition, it is shown analytically that a 
fifth or higher harmonic of such a magnitude that the orbit 
elements of certain satellites are substantially affected cannot 
be discarded a priori. 

It is well known (1) that the third harmonic causes the 
perigee of a satellite to move up and down, concurrently 
with the latitude of the perigee. In the case of Vanguard, this 
motion amounts to about 5 miles and in such a fashion that 
it always reaches its highest location at extreme southern lati- 
tude and its lowest location at extreme northern latitude. 
Because of the slow decay of the Vanguard orbit, the motion 
could be observed very clearly and repeatedly and substan- 
tiated beyond doubt. 

An interesting feature of the effect of the third harmonic 
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on the perigee radius is that, for all practical purposes, th: 
rate of change of this parameter is proportional to the rate o/ 
change of the argument of perigee. Both are strongly d:- 
pendent on the orbit inclination, but in essentially the sam: 
fashion. As a result, the cumulative effect of the thir! 
harmonic on the perigee radius is independent of the inclina- 
tion and only a function of the instantaneous latitude of th: 
perigee. This is true whether the inclination is greater or 
smaller than the critical inclination of 63°4, at which th 
argument of perigee is stationary. 

Such a proportionality does not exist, however, between the 
rate of change of the perigee radius which results from any highe’ 
odd harmonic in the gravitational field, and the rate of change 
of the argument of perigee. The cumulative effect of 
such higher odd harmonies on the perigee radius can, there- 
fore, be surprisingly high, if the orbit inclination happens to 
be close to 63°4, when the perigee changes extremely slowly 
in latitude and, as a result, accumulation continues in one 
direction for a considerable length of time. 

This fact was recognized when an effort was made to ac- 
count for certain indications that the perigee of 1958 Delta 
Two (Sputnik III, satellite) was lower than expected in 
southern latitudes, an effect opposite to the one caused 
by the third harmonic, as established from 1958 Beta 
(Vanguard I). Assuming the third harmonic to exist 
beyond doubt, it was conceived that an additional effect 
existed, which would be very small in the case of Vanguard, 
but relatively large in the case of Sputnik, so that it could 
override the effect of the third harmonic. The fifth harmonic 
was found to fulfill this requirement. 

Neglecting terms in ¢? and short period terms, the perigee 
height can be expressed as 


= [ig +% | sine 


where 
5 8 — 28 sin? y + 21 sin‘ y 
ae 4 — 5sin? y [2] 
where 
Mp = (r,/Te) | 
Pe = 


re = equatorial radius of Earth 
perigee radius 
semilatus rectum 
6, = latitude of perigee 
y = orbit inclination 
A», A; and A; are gravitational constants in the following ex- 
pression for the gravitational potential 


35 


Ae (2 gins ¢ + sin‘ e)+ 


4s (2 sin — » sin? + sin® e) | [3] 
p 9 
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Fig. 1 Perigee radius vs. number of revolutions 


Numerical evaluation of Equation [2] then yields 


&, = +0.061 for y = 33°25 (in the case of Vanguard) 
&, = +0.832 for y = 65°15 (in the case of Sputnik) 


Consequently, the effect of a fifth harmonic of Earth’s gravi- 
tational field on the perigee height is more than 10 times 
stronger on the Sputnik than it is on the Vanguard. 

Fig. 1 shows the change in perigee radius of 1958 Delta Two 
according to the best information known to the author at 
the time of writing. Unfortunately, hardly any values are 
yet published on the first 4000 revolutions of this satellite. 
The initial perigee radius, however, seems to have been at least 
1.0326 r, and possibly considerably more. The solid line repre- 
sents values for r,, computed from the known period decay 
and a model atmosphere which satisfactorily fits current 
satellite data (see Fig. 2). In this computation neither the 
third nor a fifth harmonic was considered. The points and 
crosses are the best available actual data. Even if consider- 
able allowance is made for uncertainties in these data, the 
impression is gained that the perigee tends to be lower in 
southern latitudes and higher in the northern ones. A satis- 
factory fit could be obtained in the second half of the satel- 
lite’s lifetime if a somewhat lower scale height is used in the 
atmospheric model. The initial value for r, would then, 
however, be approximately 1.0315 r., which is very much 
lower than any published value. Of course, the possibility 
that the scale height was considerably higher during the first 
5000 revolutions and lower during the last, cannot be ex- 
cluded, but there does not seem to be any specific reason for 
such an assumption. 

A critical evaluation of the period decay, as derived from 
Space Track bulletins, indicated a higher relative decay in 
those time periods when the perigee was in moderate southern 
latitudes (2). This indicates that the perigee was relatively 
low in southern latitudes. Again, explanations based on 
atmospheric changes are fully admissible, but not specifically 
indicated. 

Some numerical considerations with respect to the third 
and fifth harmonic may now follow. Assuming Az = 1.6234- 
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10~* and rg = 1.312 r, for the Vanguard and 1.110 r, for the 
Sputnik (as an average), it is found that the displacements 
Au of the perigee height in the southern apex are 


Au = +0.000624 for Vanguard 
Au = +0.001006 for Sputnik 


if 4s = —6.0-10-* and A; = 0, in accordance with present- 
day assumptions. If, however, we would assume tentatively 
A; = —6.5-:10* As = +4.5-10* [5] 


then the values of Au are 
Au = +0.000621 for Vanguard 


Au = —0.000609 for Sputnik 
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Fig. 3 Deviations of equipotential at sea level 


Consequently, the assumption [5] would fit both the observed 
perigee motion of the Vanguard and the suspected inverse 
motion in the perigee of the Sputnik. 

Finally, then, one may ask how this fifth harmonic would 
affect the equipotential surface of Earth, the so-called “pear- 
shape.” Fig. 3 shows the radial displacements, due to A; = 
—6.0-10~* and to A; = —6.5-10-* and A; = 4.5-10~*. The 


difference is surprisingly small. In fact, it is so small that it is 
almost harder to believe that a random anomaly in the poten- 
tial would exactly follow the solid line, than that it would de- 
viate from it by an amount of the order of the difference indi- 
cated by the dotted line. The conclusion must be then that 
there may be a number of harmonics of order higher than fou,, 
slowly decreasing in magnitude with increasing order numbe,, 
which substantially affect the long period orbit changes of 
satellites, inclined at angles reasonably close to 63°4. The 
purpose of this note is primarily to emphasize this finding 
and stimulate further investigation of higher harmonics and 
their effects. 

Observations between June and September 1960 on 1900 
Epsilon One (Sputnik IV) yielded a confirmation of the fin:'- 
ings reported here. In addition, however, observations on 
1960 Gamma One (Transit I, rocket) seem to indicate th:t 
also higher order odd harmonics (seventh or higher) are i: 
volved. 
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High Temperature Equilibria Involving 
Metallic Halides 


L. J. GORDON! 


Aerojet-General Corp., Solid Rocket Plant, 
Sacramento, Calif. 


HEN high temperature equilibria are calculated for a 

system containing a multivalent metal, such as Be, Mg, 
B, Al, C, Si, Ti, ete., and more than one halide, F, Cl, etc., a 
considerable number of mixed halides must be assumed as 
products. It is shown here that the mixed halides need not 
be considered individually, as they may be included in the 
equilibrium gas composition simply as multiplicity factors of 
the primary halides. Depending on the iterative technique 
and on the memory limitations of the automatic computer 
used for the calculation, this method can provide a saving 
of calculating efforts. The method is illustrated as fol- 
lows: For a trivalent metal M and the halides F and Cl, 
regardless of what other elements are present, the metal 
halide vapors which normally would be assumed are 


MF, 


MF; 
MF.CI\ molecular dissociation 


species VF products 
1cl 


One can safely assume, and in fact the data for the mixed 
halides were probably obtained from the assumption, that 
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AH and AF for the following reactions are zero 


2 MF, + + MCI; = MF.Cl [la] 
1 MF, + 3 MC); = MFCI, {1b] 


Therefore, the quantity of /F.Cl and MFC, present can be 
written as functions of the quantities of primary halides 


MF.CI = (MF3)*/2 (MCI3)'/* [2a] 
MFCl, = (MF3)'/* (MCl3)*/* [2b] 


But now as for the contribution of the mixed halides to the 
H and S sums as well as to the number of moles of gas present, 
Equations [1] may be taken as if the equilibrium signs were 
equal signs, yielding 


MF; = MF; + 2 (MF;)* (MCI) + 
3 (MF;)'* [3] 


where 
MF;, MCl; = quantities of these species calculated from 
other equilibria 
MF; = effective total quantity of that species 


Writing a similar expression for MCI; and then rearranging 
terms yield 


MC); = MFI; [1 3 (Fa) + 3 ] [4b] 


Similarly for the divalent dissociation products one obtains 
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the expressions 


MF, 


1 


1 (MF, \'4 
+3 | [5b] 


It is not immediately obvious that calculation of the com- 
plicated expressions within the brackets presents savings over 
straightforward calculation of FCI, and MF-Cl, ete., but 
since in most computer programs H — § data and equilibrium 


MCI, 


constants are calculated as needed from stored C, — T data, 
the saving is obtained by eliminating both storage space and 
mathematical integration of the C, data for the mixed halides. 
When more than one multivalent metal is considered in the 
system, the expressions within the brackets may be written 
as generalized computer sub-routines realizing that the ex- 
pressions used for a bivalent metal are the same as those 
used for the dissociation products of tri- and tetravalent 
metals, ete. 

Development of expressions for cases such as that of a 
tetravalent metal with three halides is left to the reader. 


Heat Storage System for 
Satellite Powerplants 


Cc. L. WALKER,' J. D. MATCHETT? and 
R. J. HANOLD® 


Allison Div., General Motors Corp., Indianapolis, Ind. 


in satellites using a solar energy source, the energy must 
be stored for use during the shaded portion of the orbit. 
A system has been designed for storing energy as heat, 
making use of the high heat of fusion of lithium hydride. 
This heat storage system, for use specifically with heat 
engines or thermoelectric energy conversion systems, will 
potentially store enough energy to supply up to 75 w-hr 
of usable electrical power per pound of storage system 
weight. 


AX EARTH satellite may spend a large portion of orbital 
time shaded from the sun. This is also true of satellites 
orbiting other planets in our solar system. Several studies 
(1-4)4 have shown that for satellite power requirements 
up to a few kilowatts, solar energy provides the best energy 
source for long missions. Thus it is imperative that energy 
be stored for use during the time the vehicle is shaded from 
the sun. 

The conventional energy storage device has been the electri- 
cal storage battery (5,6). Other concepts have included 
regenerative fuel cells (7) and even flywheels. However, 
many of the solar powered satellite designs use a thermal to 
electrical energy conversion device, such as a heat engine 
or thermoelectric converter. Therefore, it seemed feasible 
that heat storage might be used rather than some other 
method of energy storage. An obvious advantage would be 
that a heat engine could run continuously without stops and 
starts. This study has not only verified the feasibility but 
also has indicated that a decided weight advantage also exists 
over other energy storage components. 

The storage of energy as heat is, of course, not a new idea; 
a considerable amount of work has been conducted in con- 
junction with solar home heating research. Nevertheless, 
until recently heat storage as applied to space vehicles has re- 
ceived only slight attention in the literature (2,3). However, 
heat storage is not without its problems. Thermal energy 
may be stored and released by several mechanisms: Sensible 
heat, heat of fusion, heat of vaporization and heat of trans- 
formation. These possibilities must be balanced against 
system requirements which include high storage capacity 
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per unit weight, operating temperature in a range compatible 
with conversion system, thermal stability, materials com- 
patibility, ete. 

The heat of fusion method was chosen as the most desirable 
because it avoided the bulk of vaporizing and condensing 
but still accepted and rejected energy at a constant tem- 
perature. This choice was strengthened when an investiga- 
tion of several elements and compounds disclosed that the 
compound lithium hydride had a heat of fusion higher per 
pound than any other substance. In addition, its melting 
temperature was usable with several thermodynamic cycles. 
Data on the properties of lithium hydride have been obtained 
from the survey reports of Gibb and Messer (8), and from 
communications with the Lithium Corporation of America 
(9,10). A summary of the thermal properties of lithium hy- 
dride are shown in Table 1. 

The heat reservoir — conceived as a group of indi- 
vidually encased capsules of LiH — takes heat from a trans- 
fer fluid (probably one of the alkali metals or their alloys) as it 
comes from the solar collector and rejects the fluid at nearly 
a constant temperature. During the shaded portion of the 
orbit this heat will be returned to the fluid, so as to maintain 
a constant fluid temperature or a constant heat rejection 
rate. A proposed solar absorber, heat reservoir and heat 
transfer loop are shown in Fig. 1. 

The individual capsules were designed with a conical shape 
and were placed so that the thickest portion was exposed to 
the fluid inlet where heat transfer rates would be the greatest. 
With this arrangement the fusion front will approach the 
centerline of the capsule uniformly for maximum utilization 
of the storage material. Subsequent theoretical calculations 
have supported this hypothesis, although the cone angle will 
be smaller than expected. Initially it was feared that, in the 
absence of gravity, the hydrogen “bubble” which is necessary 
to prevent dissociation would locate itself in such a way as to 
impede heat flow. It was also feared that cracking upon 
solidification would also interfere with the heat flow. For- 
tunately, recent casting tests indicate that the void locates 
itself in the center of the capsule in the presence of gravity, 
apparently from surface tension. Very few cracks developed, 
and those were in a radial direction. 

In the design shown in Fig. 1 the lithium hydride accounts 
for approximately 60 per cent of the total heat reservoir 
weight. The energy storage capacity of the LiH is 310 


Table 1 Properties of lithium hydride 


Melting point 1256F 
Heat of fusion 1065 + 135 Btu/Ib 
Density at 20 C 48.5 lb/ft® 
Density of melt 36.1 lb/ft 
Specific heat 1.0 Btu/IbF 
Thermal conductivity 4.2 Btu/hr ft F 
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Fig. 1 A proposed heat reservoir system 


w-hr per lb but the effective capacity of the reservoir is then 
only 190 w-hr per lb. A typical mercury vapor turbine 
engine may have an efficiency of 10 per cent, whereas a re- 
cently developed heat engine operating on the Stirling cycle 
has exhibited efficiencies as high as 40 per cent. Therefore, 
in this range of conversion system efficiencies a heat reservoir 
will store from 19 to 75 w-hr of usable energy per pound of 
reservoir weight. 
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Shock Tube Design for Producing High 
Gas Temperatures 


HUGH N. POWELL' 


Cornell University, Ithaca, N. Y. 


The ideal performance of a new shock tube design is 
analyzed on an ideal gas constant cy basis and compared 
to the standard design. It is shown that within the same 
experimental limitations, the shock temperature may be 
increased by factors up to 3, or the shock density by fac- 
tors of 10 or larger with smaller but sizable gains in both 
temperature and density under most conditions of interest. 
Most standard shock tube modifications are also applica- 
ble to the new design. 


INCE the paper in 1952 by Resler, Lin and Kantrowitz 
(1)? which demonstrated the feasibility of using shock 

tubes to produce high temperatures, the shock tube has be- 
come an essential tool for high temperature gasdynamic and 
physical chemical research. Experimentally, the maximum 
shock strength (and therefore shock temperature) which is 
attainable with any particular driver and driven gas com- 
bination is limited—via the driver to driven gas acoustic 
velocity ratio—by the means of generating and sustaining 
a reasonably uniform high temperature in the driver gas 
chamber; and—via the driver to driven gas pressure ratio— 
by the chamber’s pressure limitations. The latter also 
directly limits the attainable shock density at a given shock 
T. One or more of the following techniques are frequently 
used to augment the performance within a given set of ex- 
perimental limitations: 

1 An area contraction between the driver and driven gas 
sections. 

2 The ignition of a combustible mixture to produce the 
desired high T and p driver gas. 


Received May 23, 1960. 
1 Associate Professor, Sibley School of Mechanical Engineering. 
2 Numbers in parentheses indicate References at end of paper. 


980 


3 The use of multiple diaphragms and a buffer gas be- 
tween driver and driven gases. 

Reference (2) gives a convenient set of ideal performance 
curves for the basic shock tube and the foregoing modifica- 
tions. 

The purpose of this paper is to outline a new shock tube 
configuration which permits quite substantial increases in 
the values of shock T and p to be achieved without increasing 
the driver T or p. For the sake of simplicity, ideal perform- 
ance is assumed and ideal, constant specific heat driver and 
driven gases are used. Although not completely accurate, 
this basis does provide a good measure of the energies involved 
and constitutes a simple means of comparative performance 
evaluation. 

Figs. 1 (a) and (b) show schematically the conventional 
basic shock tube design, termed configuration I, which will be 
used as a basis of comparison. Fig. 1 (c) shows the pro- 
posed configuration II before the diaphragm has ruptured. 
The test (i.e., driven) gas flows isentropically through a con- 
verging-diverging nozzle from state (0) to state (1) where it 
has a Mach number M,*. An asterisk denotes quantities 
which are defined relative to laboratory (shock tube) co- 
ordinates. Before entering the exhaust line (say to a dump 
tank), it passes through an oblique shock attached to the 
crown formed by the edges of two metal flaps which are held 
together, say by a small piece of adhesive. These flaps are 
hinged on their left edges near the diaphragm. When the 
diaphragm ruptures, the expanding driver gas blows apart 
the flaps which then close off the exhaust ducts. The driver 
gas then drives a shock into the oncoming test gas as is shown 
in Fig. 1 (d). An alternative picture of the transition from 
Fig. 1 (c) to (d) is to assume that the flaps part and close 
off the exhaust ducts shortly before the diaphragm ruptures. 
Then the test gas stagnates between the diaphragm and a 
shock moving into the state (1) gas stream. When the 
diaphragm ruptures, a second shock is formed which over- 
takes and coalesces with the first one, again yielding the 
situation pictured in Fig. 1 (d). Insofar as shock formation 
is concerned, the flaps act somewhat like a second weak 
diaphragm in a multiple diaphragm tube; other mechanical! 
arrangements are possible. 
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Fig. 1 Schematic representations of configurations I and II 


For convenience let the subscripts is and sh identify the 
isentropic and normal shock steady-state flow functions of 
the stream and shock Mach numbers, M,* and M,,, respec- 
tively. A combination of these functions with the equation 
for isentropic unsteady flow in a constant area duct yields 


Po 
(p1/Po)is (P2/ Pr) sr 


1 2 ( 1 ) 

p— | (  — — ) — 
2(a4/ ao) yot 1 sh 

[1] 


which expresses the overall pressure ratio p4/o as a function 
of M., and M,*, the driver to driven gas acoustic velocity 
ratio a4/do, and the driven and driver gas specific heat ratios 
7, and yo, respectively. Equation [1] reduces to the stan- 
dard equation for configuration I when M,* = 0. 

Fig. 2 shows 72/7» as a function of ps/po with a@s/ao as a 
parameter for the case yo = 1.4 and y, = $ with M,* = 2 
in the two configurations. The configuration II curves 
originate at the ps/po such that the M@,* = 2 stream is just 
stagnated. In the steeper portions of the as/a) curves, con- 
figuration II is seen to increase T,/T) by factors up to 2 
or so relative to configuration I for the same experimental 
(To, po) and (7's, ps) states. At ps/po = ©, the degree of 
improvement can be determined from the denominator of 
Equation [1] which then becomes zero. The resulting values 
of T2:/T> are shown in Fig. 3 as a function of M,*. It is 
seen that the degree of improvement is substantially less 
than at normal (p4/po). 

An expression for the shock density p2 can be obtained 
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Fig. 2 Comparison of configurations I and II with respect to 
shock temperature. Quantities pertinent to a representative 
shock state (2) are also shown 
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Fig.3 Limiting shock temperature (at infinite driver pressure) as 
a function of driven gas nozzle Mach number 
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from Equation [1] and the gas law 


Pa Ps T? 


—1 2 ( 1 ) 
——— Ma —— ) — 
{! 2(a4/a) (7). 1 Msn 


(2) 


The 72/7 values follow from M,, and M,* via Equation 
[1]. Fig. 4 shows 72/To vs. with as/ao and 
Ps/Po a8 parameters for configurations I and II with ,* = 
2. It is seen that configuration II may also yield substan- 
tially greater shock densities, p2 for the same experimental 
(To, po) and (7, ps) states or alternatively, a markedly in- 
creased p2 at the same T2. This feature is especially useful 
in connection with optical measurements, the sensitivity 
of which are limited by density. 

The general increase in accessible 7. and pz is the result of 
two effects, both of which stem directly from the use of a 


nozzle. 


100 
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Fig. 4 Comparison of configurations I and II with respect to 
shock temperature and density 


The first effect is simply that beginning with any given (po, 

0) state a nozzle expansion followed by a shock to a given 
T. is more irreversible than a shock by itself to 72 (as can be 
clearly scen on a T-S diagram). Consequently, for the same 
Toay = Tot), Peay is less than pea. Alternatively, for the 
same = pear, is greater than Toa). Therefore, 
less driver pressure is needed for the same 72, or the sam: 
driver pressure gives a higher 7». 

The second effect stems from the strong cooling effect o! 
the expansion on the acoustic velocity of the state (1) gas a 
compared with the corresponding state (0) gas in configura- 
tion I. The increase in /,, which is necessary to compen 
sate the cooling to attain a given shock T; is largely offse: 
by the —M/* term in the denominator of Equation [1] sv 
that the VW(7T;/T>)is factor still causes a net increase in th: 
quantity within the braces. Although this increase is not larg: 
the effect on ps/po is considerably amplified by the hig! 
power, 2y4/(7y4 — 1) = 5, to which this quantity is raised 
Because of these effects, both 72 and especially p2 are marked|: 
greater than if the driven gas pressure in configuration I wer: 
merely reduced by a factor equal to the pressure ratio acros: 
the nozzle, p:/p. = 0.128; i.e., regardless of the po valuc 
performance can always be increased by a nozzle expansioi: 
from po (assuming that pp is not so low that boundary laye: 
effects are dominant). 

If /,* is increased, the corresponding Fig. 2 plots continue 
to show sizable increases in 72/7 > relative to configuration 
I (e.g., by factors of 3 for 4,* = 3), but these gains are in- 
creasingly at the expense of ps, so that the matching Fig. 4 
plots show only small improvements in density over those 
already given for M,* = 2. 

Certain additional factors should be noted: 

1 It can be shown that the highly disruptive type of 
shock-boundary layer interaction caused by spillage of high 
pressure shock gas into the boundary layer will not occur 
when M,, > 5.98 (for 4/,* = 2)—a condition which is satis- 
fed under most conditidns when such a tube would be used. 

2 Ifthe tube L/D is large so that it must be made slightly 
divergent to allow for frictional effects on the supersonic 
stream, then the shock experiences a converging cross section, 
and a compensating effect for frictional attenuation of the 
shock is also realized. 

3 Configuration II introduces no modifications in the 
operation of the driver gas; hence, the standard shock tube 
modifications mentioned at the beginning of this note are 
equally adaptable to configuration II with equivalent. in- 
creases in performance. 

In summary, it is concluded that for the same driven and 
driver gas states, the configuration II design permits sizable 
increases in 7 or ps or both to be achieved relative to con- 
figuration I. In effect, these performance increases are pur- 
chased by the use of more apparatus, e.g., configuration I] 
can be regarded as a combination of a shock tube driver with 
a small supersonic blowdown wind tunnel. However, there 
appear to be no outstanding difficulties of principle or prac- 
tice in applying such a design. 
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Optimum Contours for Propulsion 
Nozzles! 


ELDON L. KNUTH? 
The Marquardt Corp., Van Nuys, Calif. 


PTIMUM designs of missiles and space vehicles fre- 

quently require maximum thrust to weight ratios for 
the nozzles. Foelsch (1)* describes an approximate method 
o! constructing the shortest nozzle producing maximum 
thrust for a given area ratio; Guderley and Hantsch (2) 
present an exact method for designing the nozzle giving 
maximum thrust for a given length; whereas Rao (3,4) gives 
an exact method for constructing the nozzle producing maxi- 
nium thrust for a given length and a given throat contour. 
Although Rao (3) states “length is a fair indication of nozzle 
\cight,” the decision of these authors to use a fixed length 
as a constraint was influenced also by the difficulty of exactly 
formulating the problem, using the methods of the variational 
cleulus, if some other constraint (such as fixed surface 
area) were used. In the present paper, relatively simple 
approximate methods for constructing nozzles producing 
maximum thrust for a given length and nozzles producing 
maximum thrust for a given surface area are developed. It 
is believed that, in many cases, surface area supersedes length 
asa measure of weight. 

If the fluid flow distribution at a given cross section of a 
nozzle is independent of the wall contours at all other cross 
sections of this nozzle (a condition satisfied most nearly by 
slender nozzles), then the optimum contour at the given 
cross section also is independent of the contours at other 
cross sections, and the optimum contour of the nozzle may 
be obtained by optimizing independently the contour for 
each differential element of a nozzle. Since the thrust 
contribution dz by an element of the nozzle is given by 


dr = (Pw — Pa) 2rR tan 0 dx 
one may write for thrust per unit length 

dr/dx = (pw — pa) 27R tan 9 
and for thrust per unit surface area 


dr/(2m R sec 0 dx) = (pw — pa) sin 


Pw = pressure at nozzle wall 
Pa = ambient pressure 


R = radius 
6 = angle between nozzle wall and axial direction 
x = distance along axis 


The thrust per unit length (or per unit surface area) is a 
maximum if its variation with respect to a change in 6 van- 
ishes, i.e., if for thrust per unit length 


(dpw/d@) tan 6 + (pw — pa) sec? 6 = 0 
and for thrust per unit area 
(dpw/d0) sin 0 + (po — pa) cos 6 = 0 


The variation dp../d# may be written 
_ ( dpw\(dMw 
do \dM../\ de 
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with M, = Mach number adjacent to nozzle wall. For the 
case in which the flow is isentropic and the isentropic exponent 
+ is constant 


If the slope of a differential wall element is altered, then the 
Mach number corresponding to this differential wall element 
is changed according to the Prandtl-Meyer relation 

dM. _ Mefl + — 


do (M,? — 


Substituting into the equations describing conditions neces- 
sary for thrust maxima, one finds that maximum thrust per 
unit length is obtained if 


Pu — Pe\{ (M.* — 1)'” 
99 = é 
sin 20 = 2 ( [la] 


and that maximum thrust per unit surface area is obtained if 


w — Pa 
Pw yM 


Equation [la] is the same equation obtained by Guderley 
and Hantsch (2) and by Rao (3) for the diterential element 
at the downstream end of the nozzle; Equation [lb] is be- 
lieved to be new. 

The wall slope immediately downstream of the nozzle 
throat may be computed exactly for the case in which the 
Mach number is uniform across the throat section and equal 
to unity. Then, immediately downstream of the throat, the 
angle turned by the fluid since leaving the point at which the 
Mach number is unity is given by the Prandtl-Meyer function 


(y — 1) (++ 1) 


tan-!(M,2 — 1)' [2] 


Solving this equation simultaneously with one of the two 
forementioned equations describing conditions necessary for 
thrust maxima, one obtains the wall slope immediately down- 
stream of the throat. (E.g., for the case in which the thrust 
per unit surface area is maximized, if y = 1.230 and p, = 0, 
then immediately downstream of the throat, 1, will be 
equal to 1.723 and @ will be 21.0 deg.) 

To the approximation that fluid properties are constant 
over any given cross section of the nozzle, an extremely 
simple and yet useful, approximate procedure for the rapid 
design of optimum contours for propulsion nozzles results. 
The following steps are found to be convenient: 

1 Compute values of wall slope and Mach number 
immediately downstream of the throat, using Equations 
[1 and 2]. 

2 For several values of WV greater than the Mach number 
immediately downstream of the throat, compute correspond- 
ing values of R, using the relation for one-dimensional flows 
in nozzles, and corresponding values of p, using the isentropic 
relation. 

3 For the several values of 7, compute corresponding 
values of the wall slope R’, using Equations [1]. 

4 Compute distances between the cross sections corre- 
sponding to the several values of M, using the relation 

= AR/R' avg. 

Several optimum contours calculated by this procedure 
are presented in Fig. 1. Curves A and B indicate the de- 
pendence of the optimum contour on the isentropic ex- 
ponent, whereas curves B and C compare the result of opti- 
mizing thrust per unit length with the result of optimizing 
thrust per unit surface area. Curve B appears to agree 
with the curve given by Guderley and Hantsch (2) in their 
Fig: 8. 
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This design procedure may be used also for nozzles with 
plugs. If one desires to use a throat with a given contour, 
then it is suggested that one begin with a contour calculated 
by the design procedure described and replace the sharp 
throat by the desired throat, the throat section and the re- 
mainder of the nozzle being joined with continuous slope. 
For cases involving real fluids, Equations [1] are unchanged. 
but Equation [2] and the area-Mach number relation used 
in step 2 must be replaced by the appropriate differentia! 
equations; wall contours are calculated using numerical! 


LENGTH 
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Fig. 1 Optimum nozzle contours for two isentropic exponents 
and two constraints (ambient pressure = zero) 


1955, pp. 305-313. 
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Thermodynamics, by Herbert B. Callen, 
John Wiley and Sons, Inc., New York, 
1960, 376 pp. $8.75. 

Reviewed by SercE GratcH 
Northwestern University 


It is a pleasure to find a textbook in 
thermodynamics which combines clarity 
of presentation and logical rigor with a 
modern point of view. Dr. Callen has 
been extremely successful in demonstrat- 
ing that his subject can be presented in a 
mature, rigorous manner without an 
unusual amount of mathematical com- 
plexity. 

This textbook develops thermodynamics 
from a few basic postulates, chosen on the 
basis of the resulting simplicity in the 
subsequent derivations. The first postu- 
late essentially establishes the existence 
of equilibrium states and the existence of 
internal energy. The ideas usually ex- 
pressed in the statement of the first law 
of thermodynamics are then developed 
easily as part of the discussion of the 
measurability of energy. The second 
postulate is essentially the second law of 
thermodynamics, stated in the form of the 
entropy maximum principle. The third 
postulate completes the usual content 
of the second law by asserting the addi- 
tivity of the entropy of composite systems 
and by requiring entropy to be a mono- 
tonically increasing function of energy, 
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thus ruling out negative temperatures. 
The fourth postulate is essentially the 
third law of thermodynamics. 

Dr. Callen’s presentation is essentially 
deductive. The postulates are not de- 
rived inductively from experimental ob- 
servations, but rather their validity is 
established a posteriori by demonstrating 
that they invariably lead to conclusions 
consistent with experimental observations. 
This approach leads to a considerable 
simplification of the mathematics and to a 
sharper emphasis on fundamental ideas. 

Dr. Callen is extremely careful to rest 
the quantitative development of the 
theory completely on macroscopic argu- 
ments, thus preserving ‘‘thermodynamics 
as a self-contained macroscopic discipline, 
autonomous and distinct from statistical 
mechanics.’’ However he uses the ideas 
of statistical mechanics to advantage by 
very lucid elementary qualitative discus- 
sions designed to make the postulates of 
thermodynamics intuitively more attrac- 
tive. 

The first part of the book is devoted to 
the presentation of the postulates and the 
derivation of some of their most general 
consequences. In particular, in this part 
are discussed some of the general relations 
among properties of typical substances, 
some of the basic considerations on proc- 
esses, an excellent coverage of equilibrium 
and stability of thermodynamic systems, 


and a rigorous and clear discussion of 
phase transitions. The second part of the 
book covers representative applications to 
chemical thermodynamics, elasticity, and 
magnetic and electric systems. The third 
part is devoted to fluctuations and irrevers- 
ible thermodynamics. The elegance, clar- 
ity and thoroughness of this section is not 
surprising in view of Dr. Callen’s many 
original contributions in this field. 

There is no doubt that the method of 
presentation of thermodynamics followed 
in this book is a tremendous improvement 
over the traditional Carnot-Clausius- 
Kelvin approach. Dr. Callen’s postula- 
tional approach unquestionably leads to a 
more straightforward, logical and simple 
presentation. The usual conclusions con- 
cerning efficiencies of cycles are obtained 
very easily as corollaries of the basic 
postulates. This reviewer, however, dis- 
agrees with Dr. Callen’s veiled implica- 
tion, on page 43 of his book, that the same 
weaknesses of the conventional treat- 
ments are shared by Caratheodory’s 
method. I believe that the latter method 
leads to just as logical and straightforward 
a presentation of thermodynamics as the 
postulational approach chosen by Dr. 
Callen. However, at present, there are 
no books on _ thermodynamics from 
Caratheodory’s point of view comparable 
in scope to Dr. Callen’s book. Thus the 
latter is unique in removing all the hocus 
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pocus traditionally found in books on the 
subject. This book is admirably suited 
for the advanced undergraduate students 
or first-year graduate students for whom 
it is intended, and, with appropriate omis- 
sions and with the addition of a little in- 
troductory material, it is also an excellent 
textbook for elementary courses. 


AGARD Flight Test Manual, edited by 
Courtland D. Perkins and Enoch J. 
Durbin, Pergamon Press, New York, 
1959, 2865 pp., 4 vols. $100. 

Reviewed by Erno K. Latvata 
and Associates 
ARO, Ine. 


This, the second edition of the “AGARD 
Flight Test Manual,” is another outstand- 
ing example of the North Atlantic Treaty 
Organization (NATO) technical coopera- 
tion. Over 60 authors representing 
Canada, France, the Netherlands, the 
United Kingdom and the United States 
have contributed to make a comprehen- 
sive handbook covering performance, 
stability and control, and instrumentation 
related to flight testing of aircraft. 
Previously published for official restricted 
distribution only, this new edition has 
been revised and enlarged and made avail- 
able for public sale. The manual is 
published in four loose-leaf volumes in 
special binders with annual supplements 
to be made available in 1960. Due to the 
diversity of subjects covered in the vol- 
umes, this reviewer enlisted the aid of 
four of his associates to review the volumes 
individually. 


Volume I ‘“Performance’’ contains 12 
chapters on the test and analysis proce- 
dures for performance flight testing of 
piloted aircraft. The opening chapter is 
devoted to the important topics of the 
definition of standard terms and quantities 
(such as standard atmosphere, altitude 
and air speed) and methods of measure- 
ment of pressure and temperature. Means 
for minimizing the effects of non-ideal 
instrument time response are discussed, 
and correction and calibration techniques 
are shown. Throughout this initial chap- 
ter and, in fact, throughout the entire 
volume, very careful attention is given to 
explicit definitions of quantities. When 
it is realized that the published material 
may be used in several different languages, 
the absolute necessity for clear and non- 
ambiguous definitions is obvious. 

Determination of the inflight thrust 
and power produced by turbojet, turbo- 
prop, ramjet and reciprocating engines is 
discussed in Chapter 2. Several methods 
of determining aircraft powerplant per- 
formance are given. These include direct 
measurement of thrust or torque, calcula- 
tion of thrust based on measurement of 
gas properties, and determination of 
engine performance based on climb per- 
formance of an aircraft. The errors in 
overall powerplant performance resulting 
from typical values of instrumentation 
errors are shown for a wide variation of 
engine operating conditions. 

The important problem of establishing 
methods for determining aircraft per- 
formance at “standard conditions’ from 
data obtained on the aircraft operating at 
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nonstandard conditions is discussed in 
Chapter 3. Both analytical and experi- 
mental methods for reducing test data to 
standard conditions are discussed. 

With Chapters 1 through 3 as an ap- 
propriate background, the data reduction 
techniques and performance test methods 
for turbojet, turboprop and reciprocating 
engine powered aircraft are discussed in 
Chapters 4, 5 and 6, _ respectively. 
Methods for the determination of such 
characteristics as level flight performance, 
range and endurance, rate-of-climb and 
descent are presented. The procedures 
for the determination of the performance 
of turbojet and reciprocating engine 
powered aircraft are much better de- 
veloped than those dealing with the turbo- 
prop powered aircraft. The development 
of test procedures for flight tests of turbo- 
prop powered aircraft has been hampered 
because only relatively few such planes 
have been tested and also because the 
power is produced both by a propeller and 
by direct jet thrust. Throughout these 
chapters, as well as throughout the entire 
volume, the applicable equations are de- 
rived in detail, and final equations are 
clearly presented. For all aircraft, per- 
formance testing recommended instru- 
mentation requirements are shown, and 
test pilot technique is discussed. 

Optimum climb theory and the tech- 
niques of determining climb schedules 
from flight tests are discussed in Chapter 
7. Three separate analyses of the opti- 
mum energy climb schedule are presented ; 
however, the analyses do arrive at the 
same fundamental conclusions. 

The takeoff and landing performance 
of aircraft weighing more than 4400 lb 
operating from hard surface runways is 
considered in Chapter 8. The organiza- 
tion of this chapter is similar to the 
earlier chapters in that, after the terms 
are defined and the methods of measure- 
ment discussed, the basic equations for 
the process are given, and the methods 
of determining performance at standard 
conditions are shown. 

Several unrelated special test tech- 
niques are discussed in Chapter 9. These 
special techniques include: Inflight air 
flow visualization, calibration of angle-of- 
attack and sideslip angle measuring sys- 
tems, and airbrake evaluation. 

Flight test techniques for helicopters are 
discussed in Chapters 10 through 12. The 
test techniques to establish level flight 
performance, climb and descent perform- 
ance, and takeoff and landing performance 
are presented. In addition to these, the 
specialized performance characteristics of 
helicopters relating to ground effect and 
transition performance following a sudden 
loss of power are analyzed for single rotor 
and tandem rotor helicopters. Theo- 
retical analysis of ground effects and of 
performance following sudden loss of 
power are presented and are compared to 
experimental data. 

Several appendixes, which contain fre- 
quently used tables of atmospheric prop- 
erties and graphs of some of the more 
frequently needed relationships of Mach 
number, calibrated air speed and altitude, 
are included in this volume. 

The value of this volume is somewhat 
reduced by the fact that most of the 


material is based on information pub- 
lished no later than 1953. A great many 
of the present-day high performance 
turbojet powered aircraft operate at con- 
ditions for which the assumptions given 
in this manual are not completely valid. 
Care should be exercised to re-evaluate 
the assumptions made if these techniques 
are to be applied to aircraft operating at 
Mach numbers of 2 and 3. Undoubtedly 
the primary reason for reliance on ref- 
erence material which is several years old 
was to avoid the use of classified data. 
The unique provisions made for ease of 
revision and addition of material to the 
volume should minimize the time re- 
quired to incorporate newer techniques 
for high performance aircraft as soon as 
such techniques are available. 

The authors and editors have provided 
complete listings of symbols at the begin- 
ning of each chapter in all four volumes to 
facilitate reading and understanding. 
Unfortunately, however, the symbols for 
several parameters vary from chapter to 
chapter. Since the manual is designed 
to serve primarily as a reference text for 
the study of individual sections, the nom- 
enclature variation may not create as 
much confusion as when the manual is 
read cover-to-cover. 

—Rosert E. Smita 


Volume II “Stability and Control” is 
divided into 14 chapters as follows: 
Introduction, Airplane Motions, Static 
Longitudinal Stability and Control, Flight 
Testing Aircraft for Longitudinal Maneu- 
vering Characteristics, Lateral and Direc- 
tional Control and the Measurement of 
Aerodynamic Coefficients in Steady Asym- 
metric Flight and in Flight on Asymmetric 
Power, Flight Techniques Used to De- 
termine Adequacy of Lateral Control, 
Stalling and Measurement of Maximum 
Lift, Spinning, High Speed Stability and 
Control Problems, Dynamic Response 
Techniques, Measurement and Interpreta- 
tion of Flight Test Data for Dynamic 
Stability and Control, Routine Handling 
Tests for Helicopters, Study of the 
Handling Characteristics of Helicopters 
by Harmonic Methods, A Survey of 
Flight Flutter Testing Techniques. 

Each chapter includes the flight test 
procedure and the basic equations from 
which it was derived, the pilot’s procedure, 
the required instrumentation, the data 
reduction method and the data analysis 
method. Typical flight test data are 
presented in graphical form. Detailed 
mathematical derivations are avoided 
throughout. Criteria representing desir- 
able stability and control characteristics 
are stated, and numerous references are 
included at the end of each chapter. 

The technical content of this volume 
requires that the reader be at least 
generally familiar with the technique of 
stability and control analysis. Thus, the 
volume would be useful to the engineer 
having limited experience in this field and 
of lesser value to the more experienced 
engineer except as a reference. 

—Rosert C. BAvER 


Volume III ‘Instrumentation Catalog”’ 
is intended to serve only as a guide in the 
preliminary design of flight test instrumen- 
tation systems. Although it contains well 
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over 600 pages, this ‘‘catalog’’ will defi- 
nitely not serve as a detail design hand- 
book. It was desired that current general 
usage be reflected in this compilation, and 
to this end material was solicited from 
users rather than manufacturers. 

The material is organized under the 
following broad categories: Sensing de- 
vices, electromechanical indicators, re- 
corders, radio link systems, timing de- 
vices; electrical signal modifiers, power 
suppliers, ground tracking systems, cali- 
bration equipment and data reduction 
aids. The information presented for 
each item is arranged in a concise format 
on a single loose-leaf page, which will 
facilitate the revisions which are expected 
at regular intervals. On one side of this 
page is a photograph of the particular item 
and on the other side are tabulated the 
technical characteristics, which include 
mode of operation, ranges, types of output, 
sensitivity, accuracy, linearity, resolution, 
dynamic response, cost, availability, power 
requirements, size, special features and 
typical applications. 


The most comprehensive coverage is 


given sensing devices, recorders and cali- 
bration equipment. The rather sparse 
coverage of telemetry and radio control 
equipment indicates the book was in- 
tended mainly for the aircraft field rather 
than the missile field. Obviously, a 
tremendous amount of information has 
been sorted, condensed and organized by 
the compilers of this volume, yet it should 
be stressed that it will find its greatest 
usefulness as a source of leads to more 
detailed information. 
N. MacDermorr 

Volume IV “Instrumentation System’’ 
rounds out the series by bridging the gap 
between methods of flight testing (Vols. I 
and II) and the instrumentation catalog 
(Vol. III). Volume IV is devoted to the 
concept of the total data systems, i.e., 
data acquisition, data storage, data proc- 
essing and data analysis systems. The 
material is divided into four parts. Part 
I discusses the importance of systems 
planning considering the type of data re- 
quired, number of channels, dynamic re- 
sponse, storage, data reduction capa- 
bilities, and accuracy and_ reliability 
problems with only a brief mention of 
cost problems of manpower, time and 
hardware. 

Part II is a very excellent and compre- 
hensive survey of the various types of 
data collection systems used in aircraft 
and missile flight testing. Every major 
type of data system is described and 
analyzed. These include displacement, 
temperature, pressure, acceleration, force, 
altitude, flow and icing systems, including 
signal modification, telemeter techniques 


and data storage. The manual describes 
in detail various French telemeter systems. 
Since the basic techniques described are 
similar to those used in the United States, 
the reviewer feels that the description 
given carries too much detail. Similarly 
detailed descriptions are given for re- 
corder systems though a majority of these 
products described are manufactured in 
the USA. 

Part III describes data reduction sys- 
tems, i.e., data handling and processing. 
There follows a brief description of the 
various systems and hardware available 
with a passing glimpse into analog com- 
puter techniques. This is followed by a 
reasonably thorough description of digital 
computers and their application to data 
handling. 

Part IV starts with a good description 
of the various types of errors that are 
incurred in the overall systems, and meth- 
ods of defining and eliminating these 
various errors. The second portion of 
Part IV brings forth the very important 
concept that all data is merely a sampling 
of occurrences and that the data must and 
should be treated statistically. Although 
some examples of data fairing and dis- 
tribution are given, the basic concepts of 
data distribution, systematic error, ran- 
dom error, sigma and confidence levels 
could be expounded further. 

Taken in its entirety, Volume IV is an 
excellent reference covering the broad 
ranges of instrument systems. 

—ITIra R. HEIMLIcH 


The excellent manner in which the 
contributions of many authors have been 
integrated into a single, well-organized 
reference work demands that credit be 
given to the editors. If the manual 
serves to provide a uniform basis for 
standardizing flight test techniques and 
for quoting of aircraft performance char- 
acteristics among the member nations of 
NATO, the tremendous amount of time 
and effort which has been expended to 
prepare the manual! will have been well 
spent. The broad scope of flight testing 
technology which is covered should make 
the manual a valuable reference to the 
engineer working in this and related 
fields. 

—Er1no K. Latvana 


Symposium of Plasma Dynamics, General 
editor, Francis H. Clauser, Addison- 
Wesley Publishing Co., Inc., Reading, 
Mass., 1960, 369 pp. $12.50. 

Reviewed by Rosert A. Gross 
Columbia University 


This book is a forerunner of what 
promises to be a flood of books concerning 


the newly popular subject of electrically 
conducting gases interacting with electro- 
magnetic fields. This large subject is 
literally exploding now with new results 
in many fields of science and engineering. 
Recognizing the interdisciplinary nature 
and rapid growth of this subject, the Air 
Force Office of Scientific Research and 
the National Academy of Sciences organ- 
ized a symposium in June of 1958 which 
brought together many leading scientists 
and engineers working in this field. The 
four-day meeting at Woods Hole, Mass., 
had no formal papers presented, but, 
instead, introductory speakers on each 
topic and then informal and often pro- 
vocative discussions. A transcript of the 
discussion was taken, and 10 technica! 
section editors put the proceedings in 
book form. 

The technical material covered is ver) 
broad and reflects the many disciplines 
represented at the meeting. High tem- 
perature plasma experimental research 
results are described and concern mostly 
thermonuclear fusion efforts. Sections on 
gaseous electronics and electron beams 
present the experience and viewpoint of 
electrical engineers. Sections on the 
statistical and continuum approach to 
analytically describing the behavior of 
plasmas show the power and limitations 
of each approach. Applications to high 
speed flight are briefly covered. Solar and 
finally cosmical magnetohydrodynamics 
are presented to complete the readers’ 
picture of the subject. 

This book offers a broad and advanced 
view of the field of plasma dynamics. It 
is not for the beginner in this field, since 
it shows specialists describing and defend- 
ing their work before other experts in the 
field. Nomenclature is varied and fre- 
quently unexplained, and the text often 
lacks detail which would be expected in 
informal blackboard-type discussions. 
This reviewer found it difficult reading 
at times, but when augmented by work 
and supplementary reading, it was often 
rewarding. The book does bring the ex- 
perienced reader nearly up to the present 
(i.e., to mid 1958) and gives him a glimpse 
of other areas of this big subject. One 
very unusual and worthwhile feature is an 
extensive bibliography which takes 68 
pages and includes many references from 
the “subliterature”’ of industrial and gov- 
ernment reports. This book is not recom- 
mended for the technical reader who wants 
to learn for the first time what plasma 
dynamics is about. It does present a 
remarkable chance for the specialist to 
learn what other experts in his own and 
allied fields are presently doing and 
thinking. Hence it should provide a 
stimulus for persons engaged in plasma 
research. 
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Nitrogen Tetroxide provides immediate, hypergolic igni- 
tion with amine fuels. Even after 9 years of storage, 
nitrogen tetroxide is ready for instant use. Combustion 
efficiencies closely approach the theoretical 99%. 
N.2Oy—Requires no refrigeration—can be stored indefi- 
nitely in missiles at launching sites——Can be used with 
most fuels—including those containing carbon.—Elimi- 
nates rough starts—fast reaction prevents accumulation 
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of propellants in thrust chambers.—Allows throttleable 
control of motors. 

We'll gladly supply technical literature, including: a 
59-page Product Bulletin, and a brochure entitled “Large 
Scale Handling of Nitrogen Tetroxide.” 

For specifications and local offices, see our insert in 
Chemical Materials Catalog, pages 475-482 and in Chemi- 
cal Week Buyers Guide, pages 37-44. 
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Propellants and Combustion 


Polymerization of Solid Propellants, b 
W. H. Groetzinger III, ARS Journat, vol. 
30, no. 5, May 1960, pp. 498-499. 

Propagation of Flames Supported by a 
Zeroeth-order Chemical Reaction, by 
Joseph O. Hirshfelder and Sara S. Van 
Domelen, Phys. Fluids, vol. 3, March- 
April 1960, pp. 210-215. 

Shock-initiated Detonations, by Roger 
A. Strehlow and Arthur Cohen, Phys. 
Fluids, vol. 3, March-April 1960, p. 319. 

Mass-spectrometric Studies of Ioniza- 
tion in Flames, I: The Spectrometer and 
Its Application to Ionization in Hydrogen 
Flames, by P. F. Knewstubb and T. M. 
Sugden, Proc., Royal Soc., London, vol. 
225, no. 1283, May 10, 1960, pp. 520-537. 

The Formation of a Detonation Wave 
During the Combustion of Gas in Pipes, 
by G. D. Salamandra, T. V. Bazhenova, 
Soviet Physics-Tech. Phys., vol. 4, no. 11, 
May 1960, pp. 1244-1249. 

The Interaction of Weak Shock Waves 
with the Flame Front, by G. D. Salaman- 
dra and I. K. Sevast’yanova, Soviet 
Physics-Tech. Phys., vol. 4, no. 11, May 
1960, pp. 1250-1257. 

Propellant Vaporization as a Design 
Critierion for Rocket-engine Combustion 
Chambers, by Richard J. Priem and Mar- 
cus F. Heidmann, NASA Tech. Rep. R-67, 
1960, 127 pp. 

The Interactions Between Nitrogen and 
Oxygen Molecules, by Willard E. Meador 
Jr., NASA Tech. Rep. R-68, 1960, 52 
pp. 

American Chemical Society, Literature 
of the Combustion of Petroleum (Ad- 
vances in Chemistry Series, no. 20), Wash- 
ington, D. C., The Society, 1958, 295 pp. 

Cool Flames and the Organic Reaction 

Mechanisms Involved in Their Forma- 

tion, by Cecil E. Boord, pp. 5-14. 

Fundamental Principles of Flammabili 

and Ignition, by rnard Lewis an 

Guenther von Elbe, pp. 15-21. 

Application of Chemical Reactor Theory 

to Combustion Processes, by E. C. 

Woodward, pp. 22-38. 

Detonative and Deflagrative Combus- 

tion, by Richard B. Morrison, Thomas C. 

Adamson Jr. and Alexander Weir, pp. 

69-84. 

Explosion Limit Phenomena and Their 

Use in Elucidation of Reaction Mech- 

anisms, by Walter Roth and Milton 

Scheer, pp. 85-103. 

Role of Vaporization Rate in Combustion 

of Liquid Fuels, by Gilbert S. Bahn, pp. 

104-115. 

Burning of a Liquid Droplet, by Henry 

Wise and George A. Agoston, pp. 116- 

135. 

Production and Measurement of Single 

Drops, Sprays, and Solid Suspensions, 


Epitor’s Note: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged. 


988 


Technical Literature Digest 


M. H. Smith, Associate Editor 


The James Forrestal Research Center, Princeton University 


by James A. Browning, pp. 136-154. 

Drop-size Distributions of Fuel Sprays, 

by J. Mason Pilcher and Ralph E. 
homas, pp. 155-165. 

Flame Propagation in Premixed Gases, 

by Frank E. Belles, pp. 166-186. 

Combustion of Atomized Liquid Pro- 

pellants, by C. C. Miesse, pp. 243-262. 

Nitrogen Oxide Formation in Autoigni- 
tion of * Fuel Sprays, by G. J. Mull- 
aney, Ind. Engng. Chem., vol. 52, June 
1960, pp. 529-532. 

American Society for Testing Materials, 
Symposium on Materials Research Fron- 
tiers, Boston, June 25, 1958, Philadelphia, 
The Society, 1959, 48 pp. (ASTM Special 
Technical Publication 243.) 

Modern Liquid Fuels, by Arthur L. 
Lyman, pp. 17-31. 


Materials and Structures 


Puncture Problem at Cosmic Velocities, 
by M. A. Lavrent’ev, ARS JournaLt, vol. 
30, no. 4, April 1960, Russian Suppl., pp. 
386-388. (Also NASA Tech. Transl. 
F-40, May 1960, 6 pp.) 

Satellite Pressure Losses Caused by 
Meteoroid Impacts, by M. Kornhauser, 
ARS Journat, vol. 30, no. 5, May 1960, 
pp. 475-479. 

Procedures for Electroplating Coatings 
on Refractory Metals, by John G. Beach 
and John A. Gurklis, Battelle Mem. Inst., 
Defense Metals Info. Center, DMIC Memo 
35, Oct. 1959, 26 pp. 

Temperature Flexure of Elastic Ele- 
ments, by G. A. Slomyanskii, NASA Tech. 
Transl. P34, May 1960, 20 pp. (Translated 
from Some Questions Relating to Modern 
Instrumentation Technology, edited by 
Poliakov, pp. 5-21.) 

An Investigation of Metal-ceramic Com- 
— for High-temperature Applications, 

y Sherwood W. McGee, Armour Res. 
Found., ARF 2175-1, May 15, 1959, 6 pp. 
(ASTIA AD 217, 445.) 

Oxidation of Tungsten and Tungsten 
Based Alloys, by E. A. Gulbransen, K. F. 
Andrew, P. E. Blackburn et al., Westing- 
house Electric Corp., Res. Labs., Sept. 1959, 
78 pp. 

Application of Rate-temperature Param- 
eters to Tensile Data for Magnesium 
Alloys and a Relation Between the Larson- 
Miller Constant and the Activation Energy, 
by Charles R. Manning Jr., NASA T 
D-172, May 1960, 23 pp. 

American Society for Testing Materials, 
Symposium on Materials Research Fron- 
tiers, Boston, June 25, 1958, Philadelphia, 
the Society, 1959, 48 pp. (ASTM Special 
Technical Publication 243.) 

Materials in the Nuclear Age, by 
John J. Antal, pp. 7-16. 

Some New Types of Neutron-shielding 
Materials, by Anthony N. Fasano, Wright 
Air Dev. Center, W A DC-T N-59-227, June 
1959, 12 pp. 

Vibration and Near-field Sound of Thin- 
walled Cylinders Caused by Internal Tur- 
bulent Flow, by Paul F. R. Weyers, NASA 
TN D-430, June 1960, 58 pp. 


Creep Bending and Buckling of Thin 
Circular Cylindrical Shells, by Burto: 
Erickson, Francis W. French, Sharad A. 
Patel, N. J. Hoff and Joseph Kempner, 
NASA TN D-429, June 1960, 49 pp. 

Strength Properties of Reinforced Plas- 
tic Laminates at Elevated Temperatures, 
by K. H. Boller and K. E. Kimball, 

right Air Dev. Center TR 59-229, Aug. 
1959, 54 pp. 

Temperature Distribution and Thermal 
Stresses in Structures with Contact 
Resistances, by Frederick V. Pohle, T. J. 
Lardner and Francis W. French, Brooklyn 
Polytechnic Inst., Dept. Aeron. Engng. & 
oe Mech., PIBAL Rep. 557, May 1960, 

pp. 

An Improved First - approximation 
Theory for Thin Shells, by J. Lyell Sanders 
Jr., NASA Tech. Rep. R-24, 1959, 11 pp 

A Design Philosophy for Repeated 
Thermal Loading, by E. W. Parkes, 
NATO, AGARD Rep. 213, 1958, 24 pp. 

Buckling of Conical Shells Under Ex- 
ternal Pressure and Thermal Stresses, by 
Josef Singer and Nicholas J. Hoff, Technion 
Res. and Dev. Foundation, Haifa, Israel, 
Tech. Final Rep., Dec. 1959, 78 pp. 

Some Experiments Concerning Subse- 
quent Yield Surfaces in Plasticity, by 
Harvey G. McComb Jr., NASA TN D- 
396, June 1960, 33 pp. 

Missile Structures and Materials, by 
E. Z. Gray, NATO AGARD Rep. 218, 
1958, 22 pp., 8 refs. 

Some Effects of Internal Heat Sources 
on the Design of Flight Structures, by G. E. 
A. Thomann and R. B. Erb, NATO, 
AGARD Rep. 208, 1959, 62 pp., 7 refs. 

Yield Conditions for Rotationally Sym- 
metric Shells Under Axisymmetric Load- 
ing, by P. G. Hodge Jr., J. Appl. Mech. 
(ASME Trans.), vol. 82E (vol. 27, no. 2), 
June 1960, p. 323. 

Integral-differential and Integral Equa- 
tions for the Equilibrium of Thin Elastic 
Shells, by N. A. Kiltschewskij, Zeztschrift 
ftir Angewandte Math. und Mech., vol. 40, 
April 1960, pp. 153-160. (In German.) 

Symposium on Surface Effects on Space- 
craft Materials, Palo Alto, Calif., May 
1959, ed. by Francis J. Clauss, John Wiley 
& Sons, N. Y., 1960, 404 pp. 

Effect of Surface Thermal-radiation 

Characteristics on the Temperature- 

control Problem in Satellites, by Walter 

G. Camack and Donald K. Edwards, 

pp. 3-54. 

Temperature Control of the Explorers 

and Pioneers, by T. O. Thostesen, A. R 

Nibbs and E. P. Buwalda, pp. 55-91. 

Coatings for Space Vehicles, by Joseph 

C. Richmond, pp. 92-116. 

Spectral Reflectance Measurements, by 

Robert V. Dunkle, pp. 117-140. 

Emissivity of Materials Near Room 

Temperature, by Ralph McDonough, 

pp. 141-151. 

Vanguard Emittance Studies at NRL, by 

Louis F. Drummeter Jr. and E. Gold- 

stein, pp. 152-163. 

Emissivity, Absorptivity, and High- 

temperature Measurement at Armour 


ARS JourNAL 


fe 
Bi 
1 


Research Foundation, by O. H. Olson 
and Sidney Katz, pp. 164-181. 

Some Methods Used at the National 
Bureau of Standards for Measuring 
Thermal Emittance at High Tempera- 
tures, by Joseph C. Richmond, pp. 182- 
194. 

A Method for Measuring Total Hemi- 
spherical Emissivity of Metals, by C. P. 
Butler and Edward C. Y. Inn, pp. 195- 
211. 

High-tem Emissivity Measure- 
ments at The Martin Company, by Mel- 
vin Crompton, pp. 212-219. 

Apparatus for the Measurement of 
Spectral and Total Emittance of Opaque 
Solids, by Clyde C. Shaw, pp. 220-246. 
Materials Sublimation and Surface 
Effects in High Vacuum, by G. F. 
Vanderschmidt and J. C. Simons Jr., pp. 
247-262. 

Vacuum Volatility of Organic Coatings, 
by Lt. George F. Matacek, pp. 263- 
285. 

Effects of High Vacuum on Mechanical 
Properties, by M. R. Achter, pp. 286- 
306. 


Surface Phenomena and Friction, by T. 
Bruce Daniel, pp. 307-318 
Panel Discussion: High-vacuum En- 
vironmental-test Chambers and Equip- 
—_ by M. R. Achter, T. B. Daniel, 
L. Jepsen, Lt. G. F. Matacek, N. 
Milleron: E. Perkins, J. P. Reynolds, 
J. C. Simons and G. F. Vanderschmidt, 
pp. 319-346. 
Measurement of Extreme Ultraviolet 
Solar Radiation and Associated Photo- 
electric Emissions from Solids, by 
Hans E. Hinteregger, pp. 347-363. 
The Effects of Ultraviolet Radiation on 
Organic Film-forming Polymers, by J. 
E. Cowling, A. Alexander, F. Noonan, 
R. Kagarise and 8. Stokes, pp. 364-377. 
Interplanetary Dust Distribution and 
Erosion Effects, by David B. Beard, 
pp. 378-390. 
Atomic and Molecular Sputtering, by 
Robert P. Stein, pp. 391-404. 


Fluid Dynamics, Heat 
Transfer and MHD 


Plasma Turbulence, by Leslie S. G. 
Kovasznay, Space Tech. Labs., Inc., Phys. 
Res. Lab., STL/T R-60-0000-AE038, Feb. 
1960, 22 pp. 

Transverse Pinch Oscillations, by Bur- 
ton D. Fried, Space Tech. Labs., Inc., Phys. 
Res. Lab., STL/T N-60-0000-09059, March 
1960, 8 pp. 

Hardness of the Axially-symmetric 
Magnetic Piston, by Allan B. Schaffer, 
Space Tech. Labs. Inc., Phys. Res. Lab., 
STL/T R-60-0000-09011, Jan. 1960, 14 pp. 

Surface Disturbances in Magnetohy- 
drodynamics, by S. I. Syrovatskiy, NASA 
Tech. Transl. F-23, May 1960, 45 pp. 


Compressible Flow Tables, k = 1.30, 
Isentropic Process, by C. Hoebich, 4rmy 
Rocket & Guided Missile Agency, Ordn. 
Missile Labs. Div., Design Dev. Lab., 
Redstone Arsenal, ARGMA TN 1H1N- 10, 
April 1960, 87 pp. 

Calculation of the —— of Shells at 
High Speeds, by I Sokolova, ARS 
JOURNAL, vol. 30, no. 4, Apa 1960, Russian 
Suppl., pp. 375-379. 

Concerning the Problem of the Influence 
of Internal Degrees of Freedom of Parti- 
cles on the Transport Coefficients of a 
Multi-component Mixture of Gases, by 
E. V. Samuilov, ARS JourNAL, vol. 30, no. 
Fe 1960, Russian Suppl., pp. 388- 
Q¢ 


OcToBER 1960 


Method of Calculation of the Transport 
Coefficients of Air at High Temperatures, 
by E. V. Stupochenko, B. V. Dotsenko, 
I. P. Stakhanov and E. V. Samuilov, ARS 
JOURNAL, vol. 30, no. 4, April 1960, Rus- 
sian Suppl., pp. 394-402. 

Investigation of a Stationary Discon- 
tinuity Surface in an Electro-magnetic 
Field with a Gas Conductivity Jump, by 
G. A. Lyubimov, ARS JourNAL, vo 
no. 4, April 1960, Russian Suppl., pp. 
416-420. 


Heat and Mass Transfer in Aeronautical 
Engineering, by D. B. Spalding, Aeron. 
Quarterly, vol. 11, May 1960, pp. 105-136. 

Recent Progress in Rarefied Gasdy- 
namics, by S. A. Schaaf, ARS JourNaL, 
vol. 30, no. 5, May 1960, pp. 443-447. 

Hypersonic Shock Wave - Boundary 
Layer Interaction and Leading Edge Slip, 
by H. T. Nagamatsu and R. E. Sheer Jr., 
ARS Journal, vol. 30, no. 5, May 1960, 
pp. 454462. 

International Symposium on Rarefied 
Gas Dynamics, Nice, 1959, Rarefied Gas 
Dynamics: Proceedings, edited by F. M. 
Devienne, Pergamon Press, N. Y., 1960, 
442 pp. 

Techniques of Measurement in Rarefied 

Gases, by I. Estermann, pp. 38-54 

Techniques of Measurements, Flow 

Visualization, by F. C. Hurlbut, pp. 55- 

68. 

The Application of Supersonic Beam 

Sources to Low Density, High Velocity 

Experimentation, by Hermon M. Parker, 

A. R. Kuhlthau, Richardo Zapata and 

John E. Scott, pp. 69-79. 

Theory of Rarefied Gases, by Harold 

Grad, pp. 100-138. 

Recent Investigations of the Boltzmann 

Equation, by Eugene P. Gross, pp. 139- 

150. 

The Boltzmann-Maxwell Equation as 

the Fundamental Equation of Gas Dy- 

namics, by Toyoki Koga, pp. 151-160. 

Experiment Versus Kinetic Theory for 

Rarefied Gases, by F. S. Sherman and 

L. Talbot, pp. 161-191. 


Contributions to the Theory of Almost- 
free-molecule Flows, by V. C. Liu, pp. 
239-245. 

Theoretical Solutions to Some Nearly 
Free Molecular Problems, by D. R. 
Willis, pp. 246-257. 

Aerodynamics of Rarefied Gases, by 
Ronald F. Probstein, pp. 258-275. 

The Study of Boundary Conditions in 
Slip-flow Aerodynamics, by Robert E. 
Street, pp. 276-292. 

Slip Flow over a Short Flat Plate, by J. 
A. Laurmann, pp. 293-316. 

The Drag of a Sphere Which Moves at 
High Speed Through a Rarefied Gas, by 
Pierre Welander, pp. 317-328. 
Aerodynamic Coefficients of a Pointed 
Cone at Angles of Attack in Rarefied 
Gas Flow, by V. A. Langelo and A. 
Lengyel, pp. 329-351. 

On the Fundamental Aspects of Free 
Molecule Flow, by M. Z. Krzywblocki, 
pp. 352-366 

General Considerations on the Aero- 
thermodynamics of Rarefied Gases, 
by Edmond A. Brun, pp. 367-382. 
(In French.) 

Some Remarks on the Free Molecular 
Regime, by F. M. Devienne, pp. 383- 
386. (In French.) 

Determination of Fundamental Atomic 
and Molecular Properties from Gas 
Dynamic Measurements, by Wayland 
C. Griffith, pp. 387-389. 

Study of the Different Flow Regimes in 


Rarefied Gases by the Revolving Arm 
Method, by B. Crave, pp. 390-393. 
(In French.) 

Plasma Dynamics—from Gas Dynamic 
— of View, by Shih-I Pai, pp. 394- 


A Discussion of Magnetic Boundary 

Layers with Boundary Conditions As- 

similating Combustion, Blowing or 

Sublimation at the Wall, by Paul S. 

Lykoudis, pp. 407-415. 

Continuum Theory and Rarefied Hyper- 

sonic Aerodynamics, by Ronald F. 

Probstein, pp. 416-431. 

Numerical Comparison Between Exact 
and Approximate Theories of Hypersonic 
Inviscid Flow Past Slender Blunt Nosed 
Bodies, by Saul Feldman, ARS Journat, 
vol. 30, no. 5, May 1960, pp. 463-468. 

Solar Heating of a Rotating Cylindrical 
Space Vehicle, by A. Charnes and 8. 
Raynor, ARS Journat, vol. 30, no. 5, May 
1960, pp. 479-484. i 


Flight Mechanics 


The Long-time Satellite Rendezvous- 
Trajectory, by Louis W. Spradlin, Aero/- 
pace Engng., vol. 19, June 1960, pp. 32- 


Orbits for Planetary Satellites from 
Doppler Data Alone, by Armin J. Deutsch, 

RS JOURNAL, vol. 30, no. 6, June 1960, 
pp. 536-542. 

Angular Displacement of a Magnetized 
Sphere by the Terrestrial Magnetic Field 
During Vertical Descent, by Marvin H. 
Hewitt, ARS JourNAL, vol. 30, no. 6, 
June 1960, p. 558 

An Approximate Analytic Solution of 
Re-entry Trajectory with Aerodynamic 
Forces, by Kenneth Wang and Lu Ting, 
ARS Journat, vol. 30, no. 6, June 1960, 
pp. 565-566. 

Elliptical Orbit Characteristics, by G. S. 
Gedeon and B. J. Greer, ARS JouRNAL, 
vol. 30, no. 6, June 1960, pp. 570-572. 

The Orbits of Cosmic Rockets Towards 
the Moon, by L. I. Sedov, Astronautica 
Acta, vol. 6, no. 1, 1960, pp. 16-31. 


Vehicle Design, Testing and 
Performance 


Manned Space Flight—Present and 
Future Steps, by Robert R. Gilruth and H. 
Kurt Strass, Aero/Space Engng., vol. 19, 
May 1960, pp. 16-17, 88. 

The Military Test Space Station, by 
Col. Lowell B. Smith, Aero/Space Engng., 
vol. 19, May 1960, pp. 18-19. 

Atmospheric Entry of Manned Vehicles, 
by Carl Gazley, Aero/Space Engng., vol. 
19, May 1960, pp. 22-23, 90. 

Design Criteria and Their Application 
to Economical Manned Satellites, by H. 
Hermann Koelle, Erich E. Engler and 
John W. Massey, Aero/Space Engng., vol. 
19, May 1960, pp. 24-25, 90 

Manned Scientific Orbital 
by M. Stoiko, G. G. Kayten and J. 
Dorsey, Aero/Space Exgng., vol. 19, May 
1960, pp. 26-27, 92-93. 

Some Factors Influencing the Selection 
of a Manned Space Station Concept, by 
Joseph H. Doss and Leo J. Montague, 
Aero/Space Engng., vol. 19, May 1960, pp. 
28-29, 94. 


A Modular Concept for a Multimanned 
Space Station, by Saunders B. Kramer and 
Richard A. Byers, Engnq., vol. 
19, May 1960, pp. 30-31, 

Man’s Contribution to an Operational 
Space Station Concept, by M. A. Grodsky 


989 


ee 
| 
tor 
ity 
‘hi 
hin 
tor: 
ner, 
as- 
res, 
ali 
ug. 
nal 
lyn 
60, 
ion 
ers 
op 
fed 
on 
D- 
es 
E. 
J 
2), 
a- 
tic 
oft 
0, 
q 
e 
ay 
e- 
a 
rs 
3 
mn 
j- 
|| 


and R. D. Sorkin, Aero/Space Engng., vol. 
19, May 1960, pp. 80-81 

The Australian Sounding Rockets Long 
Tom and Aeolus, by D. A. H. Bird, J. Roy. 
Aeron. Soc., vol. 64, May 1960, pp. 277- 
289. 


Studies on Dynamics and Instrumenta- 
tion of the Holloman Track, Air Force 
Missile Dev. Center, Holloman AFB, New 
Mexico, AFMDC-TR-59-8, April 1959, 
624 pp., 84 refs. 

Calculation of the Flexural Rigidity of 

the Holloman Track, by Alfred W. 

Simon, pp. 1-8. 

A Method of the Experimental Deter- 

mination of the Parameters of the 

Holloman Track, by Alfred W. Simon, 

pp. 9-22 

The Experimental Determination of the 

Fundamental Parameters of the High- 

speed Test Track (Supplementary Re- 

port), by Alfred W. Simon, pp. 23-34. 

Calculations Pertaining to a Moving 

Load on a Beam on an Elastic Founda- 

tion, by Alfred W. Simon, pp. 35-54. 

Dynamic Analysis of a High-speed 

Track, by Stephen H. Crandall, pp. 55- 

78. 


The Timoshenko Beam on an Elastic 
Foundation, by Stephen H. Crandall, pp. 
79-106. 

Deflection of the Holloman Track 
Under Static Load, by F. J. Calkum and 
Ray L. Larkins, pp. 107-122. 

The Relation Between the Bending 
Strain and the Load for a Beam on an 
Elastic Foundation, by Alfred W. Simon, 
pp. 123-138. 

Vibrations Set up by Moving Load in 
Track on Elastic Foundation, by Frank 
W. Bubb, pp. 139-200 

Track Vibration Study in Controlled 
Excitation, by Thomas P. Rona, pp. 
201-232. 

Vibration Level Measurements in the 
Supersonic Test Track, by Thomas P. 
Rona, pp. 233-246. 

Transverse Beam Vibrations in Ac- 
celerated Rocket Sleds, by Thomas P. 
Rona, pp. 247-272. 

Recommendations for the Strain Meas- 
urement Program at the Holloman 
pe by Thomas P. Rona, pp. 273- 


Statistical Analysis of Markings: Sled 
Run of 12 July 1956, by Paul R. Rider, 
pp. 281-292. 
High Accuracy Velocity Measurements, 
by Thomas P. Rona, pp. 293-324. 
Notes on Design of Sled Slippers for 
the EE Track, by S. Kyropou- 
los and M. E. Shank, pp. 325-334. 
Some Considerations Concerning Lubri- 
cation of -e for the High-speed 
Track, by M Shank, pp. 335-356. 
Performance Characteristics of Lubrica- 
tion for High-speed Tracks, by S. 
Kyropoulos, pp. 357-370. 
A Study of the Aerodynamics of High- 
speed Sleds, by Joe R. Foote and S. 
oyer, pp. 371-402. 


Some Remarks on the Optimum S 
Attainable with a Two-stage Sled, 

Joe R. Foote, pp. 403-422. 

Ram Wind Tunnel, by Harry L. Gep- 
hart, pp. 423-448. 

Preliminary Investigation of the Aero- 
dynamic Effects of Linear Acceleration 
as Related to the High-Speed Track, by 
Ray Holland Jr. and J. R. Foote, pp. 
449-572. 

Parameters for Rocket Sled Propulsion 
Systems, by H. J. Weigand and Ferdi- 
nand F. Kuhn, pp. 573-612. 


990 


Some Requirements for the Efficient 
Attainment of Range by Air-borne Vehi- 
cles, by M. Arens, J'echnion, Israel Inst. 
Tech., Dept. Aeron. Engng., Haifa, Rep. 3, 
Sept. 1959, 26 pp. 

Investigation of the Low-subsonic Flight 
Characteristics of a Model of an All-wing 
Hypersonic Boost-glide Configuration Hav- 
ing Very High Sweep, by Robert E. Shanks, 
NASA TN D-369, June 1960, 21 pp. 

The Tiros System on the Ground, by 
John E. Kiegler and Charles B. Oakley, 
Astronautics, vol. 5, June 1960, pp. 
44-45, 100-102. 


Guidance and Control 


On Automatic Controls, by S. Lefschetz, 
RIAS, Baltimore, Md., Tech. Rep. 60-9, 
March 1960, 6 pp. 

Comparison of Some Actuation Methods 
for Attitude Control of Space Vehicles, by 
Walter Haeussermann, Aero/Space Engng., 
vol. 19, May 1960, pp. 64-65. 

Attitude Control of a Space Vehicle by a 
Gyroscopic Reference Unit, by Marvin L. 
wad Lewis G. Lynch, Aero/Space 
Engng., vol. 19, May 1960, pp. 66-67. 

Attitude and Flight Path Control System 
for a Space Station Supply Vehicle, by J. 
Stalony-Dobrzanski and O. Imai, Aero/- 
Space Engng., vol. 19, May 1960, pp. 68- 
69, 100, 102. 

Orbital Rendezvous and Guidance, by 
Ernst A. Steinhoff, Aero/Space Engnq., 
vol. 19, May 1960, pp. 70-71, 106. 

Rendezvous in Space—Effects of Launch 
Conditions, by Norman V. Petersen and 
Robert S. Swanson, Aero/Space Engngq., 
vol. 19, May 1960, pp. 72-73, 106-107. 

Ascent Guidance for a Satellite Ren- 
dezvous, by T. B. Garber, Aero/Space 
Engng., vol. 19, May 1960, pp. 74-75, 108. 

A Guidance Technique for Achieving 
Rendezvous, by Philip G. Felleman and 
Norman E. Sears Jr., Aero/Space Engng., 
vol. 19, May 1960, pp. 76-77. 

Manned Control of Orbital Rendezvous, 
by E. Levin and J. Ward, Aero/Space 
Engng., vol. 19, May 1960, pp. 78-79. 

Longitudinal Range Control During the 
Atmospheric Phase of a Manned Satellite 
Reentry, by Arthur Assadourian and 
Donald C. Cheatham, NASA TN D-253, 
May 1960, 46 pp. 

On the Analysis of the Stability of the 
Periodic Modes of Operation in Nonlinear 
Control Systems with Many Degrees of 
Freedom, by V. A. Taft, Automation & 
Remote Control (trans. of Avtomatika i 
Telemekhanika), vol. 20, no. 9, Sept. 1959, 
pp. 1132-1140. 

On Designing Correcting Circuits for 
Automatic Control Systems in Accordance 
with the Mean ~~ Error Criterion, by 
V. I. Kukhtenko, Automation & Remote 
Control (trans. of Avtomatika i Telemekhan- 
tka), vol. 20, no. 9, Sept. 1959, pp. 1151- 
1160. 

Experimental Determination of Auto- 
matic Control System Links’ Transfer 
Functions by Means of Standard Electronic 
Models, by L. N. Darovskikh, Automa- 
tion & Remote Control (trans. of Avto- 
matika i Telemekhanika), vol. 20, no. 9, 
Sept. 1959, pp. 1180-1187. 

Frequency Elements of Remote Con- 
trol, by I. N. Listiskaya, V. N. Mikhailov- 
skii, K. D. Nadtochii and V. N. Okhot- 
skva, Automation & Remote Control 
(trans. of Avtomatika i Telemekhanika), 
vol. 20, no. 9, Sept. 1959, pp. 1230-1234. 

The Principles Involved in the Construc- 
tion of an Automatic Radio Direction 
Finder, by Iu. A. Ierumhimovich, Radio 
Engng. (trans. of Radiotekhnila), vol. 


14, no. 6, 1959, pp. 86-94. 


The Transient Behavior of Nonlinear 
Systems, by Francis H. Clauser, 7he 
Johns Hopkins University, Dept. Aeron., 
April 1960, 56 pp. 

Self-adaptive Control Systems, by 
George J. Thaler, Naval Postgrad. School, 
Tech. Rep. 18, April 1960, 32 pp. 

Dynamic Programming and Numerical 
Experimentation as Applied to Adaptive 
Control Systems, by M. Aoki, Univ. Cali/., 
Los Angeles, Dept. Engng., Rep. 60-—i6, 
Feb. 1960, 215 pp. 

Investigation of Third-order Contactor 
Control Systems with Two Complex Poles 
Without Zeros, by Irmgard Fliigge-Loiz 
and Tomo Ishikawa, NASA TN D-42s, 
June 1960, 69 pp. 

The Theory of Correctional Manoeuvres 
in Interplanetary Space, by D. F. Lawd«n 
and RB. Lang, Astronautica Acta, vol. i, 
no. 1, 1960, pp. 48-60 

On Evaluating the Quality of Nonlinear 
Automatic Systems with Random Inter- 
ference, by E. P. Popov, Automation 
Remote Control (trans. of Avtomatika i 
Telemekhanika), vol. 20, no. 10, Oct. 1959 
(June 1960), pp. 1281-1287. 

L. S. Pontryagin Maximum Principle in 
the Theory of Optimum Systems, by L. |. 
Rozonoer, Automation & Remote Control 
(trans. of Avtomatikai Telemekhanika), vol. 
20, no. 10, Oct. 1959 (June 1960), pp. 
1288-1302; and vol. 20, no. 11, Nov. 1959 
(June 1960), pp. 1405-1421. 

The Problem of Optimizing Systems 
which Contain Essentially Non-linear 
Elements, by E. P. Merkulova, Automa- 
tion & Remote Control (trans. of Avtomatika 
it Telemekhanila), vol. 20, no. 10, Oct. 
1959 (June 1960), pp. 1303-1313. 

Optimal Control of an Object with Two 
Controlling Stimuli, by E. A. Rozenman, 
Automation & Remote Control (trans. of 
Avtomatika i Telemekhanika), vol. 20, 
no. 10, Oct. 1959 (June 1960), pp. 1314- 
1318. 

A Method cf Synthesizing Linear Sam- 
pled-data Automatic Contrcl Systems sg 
Accord with Dynamic Criteria, by L. N. 
Vongin, Automation & Remote Control 
(trans. of Avtomatikai Telemekhanika), vol. 
20, no. 10, Oct. 1959 (June 1960), pp. 
1319-1328. 

AC Servo Systems in Which the Signal 
Depends on the Error and Its Derivative, 
by N. P. Vlasov, Automation & Remote 
Control (trans. of Avtomatika i Telemek- 
hanika), vol. 20, no. 10, Oct. 1959 (June 
1960), pp. 1329-1334. 

Graphicoanalytic Method of Computing 
Steady-state Speeds of Pneumatic Drives, 
by V. N. Dmitriev and V. I. Chernyshev, 
Automation & Remote Control (trans. of 
Avtomatika i Telemelhanika), vol. 20, no. 
10, Oct. 1959 (June 1960), pp. 1335-1341. 

Experimental Determinaticn of Power 
Positioning Servo System Characteristics, 
by M. B. Tumarkin, Automation & Re- 
mote Control (trans. of Avtomatika i Tele- 
mekhanika), vol. 20, no. 10, Oct. 1959 
(June 1960), p. 1396. 


Instrumentation and 
Communications 


Methods for the Control of Interfering 
Currents Originating at the Input of an 
Electrostatic Fluxmeter During Its Opera- 
tion in a Conducting Medium, by I. M. 
Imyanitov and Ya. M. Shvarts, NASA 
Tech. Transl. F-16, May 1960, 10 pp. 

Methods of Combatting Interfering 
Currents That Arise at the Input of an Elec- 
trostatic Fluxmeter Operating in a Con- 
ducting Medium, by I. M. Imyanitov and 
Ya. M. Shvartz, ARS JourNAL, vol. 30, no 


ARS JourNAL 


| 
ag 


4, April 1960, Russian Suppl., pp. 403-406. 

Telebit—An Integrated Space Navigation 
and Communication System, by George E. 
Mueller, ASTRONAUTICS, vol. 5, May 1960, 
pp. 26-27, 88, 90, 92, 94, 96. 

How Good Is the Lunik III Photog- 
raphy?, by Merton E. Davies, AsTRo- 
nautTics, vol. 5, May 1960, pp. 28-29, 
112-113. 

A Solid-state UHF Radar Transponder, 
by Lincoln Cartledge and G. Bradford 
Tiffany, ASTRONAUTICS, vol. 5, May 1960, 
pp. 30-31, 114-115. 

A Recoverable Interplanetary Space 
Problem, by Milton B. Trageser, AsTRo- 
nautics, vol. 5, May 1960, pp. 32-35, 116, 
118-120. 

Television Cameras for Space Explora- 
tion, by M. H. Mesner and J. R. Stanis- 
zewski, ASTRONAUTICS, vol. 5, May 1960, 
pp. 36-37, 126-130. 

Maser, Iraser, and Laser, by Harold 
Lyons, ASTRONAUTICS, vol. 5, May 1960, 
pp. 38-39, 100-104. 

Space Age Electronic Materials, by 
John W. Buttrey, Astronautics, vol. 5, 
May 1960, pp. 40-41, 121-130. 

Project Mercury Tape Recorder, by G. 
W. Boyer, Astronautics, vol. 5, May 
1960, pp. 44-45, 72, 74-75. 

X-15 Operations: Electronics and the 
Pilot, by Neil A. Armstrong, ASTRONAU- 
tics, vol. 5, May 1960, pp. 42-43, 76, 78. 

Signal Transmission Through Ionized 
Media; How Man-made Ionization Af- 
fects Communications, by W. A. Green- 
how, Electronics, vol. 33, May 20, 1960, 
pp. 81-85. 

Spectrometric Analysis of Missile 
Flights; Rapid-scan Spectrometer Pro- 
vides Function-of-time Chart, by J. N. Day 
Jr., Electronics, vol. 33, May 20, 1960, pp. 
86-88. 

Measuring the Temperature of a Gas 
Downstream from a Shock Wave, by S. A. 
Losev and N. A. Generalov, Instruments & 
Experimental Techniques (trans. of Pribory 
i Tek. Eksper.), no. 3, May-June 1959 
(published April 1960), pp. 454-455. 

A Theory of Radar Scattering by the 
Moon, by T. B. A. Senior and K. M. 
Siegel, J. Res., Nat. Bur. Standards, D. 
Radio Propagation, vol. 64D, May-June 
1960, pp. 217-230. 

Theory of the Stagnation-point Lang- 
muir Probe, by Lawrence Talbot, Phys. 
— vol. 3, March-April 1960, pp. 289- 
298. 

The Design of Automatic Gain Control 
in Transistor Amplifiers, by Ye. P. De- 
ment’ev, Radio Engng. (trans. of Radiotekh- 
nika), vol. 14, no. 6, 1959, pp. 51-59. 

Some Questions on the Dispersion of 
Radiowaves in the Ionosphere, by V. D. 
Gusev, Radio Engng. & Electronics (trans. 
of Radiotekhnika i Elektronika), vol. 4, 
no. 1, 1959 pp. 15-22. 

The Investigation of Winds and Non- 
uniformity in the Ionosphere by a Radio 
Astronomical Method, by V. V. Vitkevich 
and Ia. L. Kokurin, Radio Engng. & 
Electronics (trans. of Radiotekhnika i 
Elektronika), vol. 4, no. 1, 1959, pp. 23-29. 

Signal Detection Against a Background 
of Normal Noise and Random Reflections, 
by V. D. Zubakov, Radio Engng. & Elec- 
tronics (trans. of Radiotekhnika i Elektron- 
tka), vol. 4, no. 1, 1959, pp. 42-59. 

Maximum Noise-proofing of the Radio- 
reception of Signals with Amplitude Fre- 
quency and Phase Manipulation, by N. P. 
Teplov, Telecommunications (trans. of 
Elektrosvyaz’), no. 1, 1959, pp. 34-47. 

Bulletin of Stations for Optical Obser- 
vations of Artificial Earth Satellites, NASA 
Tech. Transl. F-19, May 1960, 50 pp. 


OcToBER 1960 


(Translated from Iskusstvennye Sputniki 
Zemli, no. 7, Moscow, Acad. Sci. USSR, 
1959. ) 

Artificial Earth Satellites, vols. I & II, ed. 
by I. V. Kurnosova, transl. from Russian, 
Plenum Press, N. Y., 1960, 120, 107 pp. 

Optical Methods for the Observation of 

Artificial Earth Satellites, by I. S. 

Shklovskii, vol. I, pp. 56-63. 

Use of the Doppler Effect for Determin- 

ing Orbit Parameters of an Artificial 

Earth Satellite, by V. A. Kotel’nikov, 

V. M. Durbrovin et al., vol. I, pp. 64-78. 

An Acoustical Method for the Measure- 

ment of the Mechanical Parameters of 

Meteorites, by M. A. Isakovich and N. 

A. Roi, vol. II, pp. 105-107. 

Pulsating Flow Measurement—Today, 
Part II, by P. H. Stirling and Henry Ho, 
Ind. Engng. Chem., vol. 52, June 1960, p. 
75-76A. 

Variable Capacitance Liquid Level Sen- 
sors, by L. Bruce Wilner, Rev. Sci. Instr., 
vol. 31, May 1960, pp. 501-507. 

Satellite Astronomical Observatories, 
by Nancy G. Roman, Aero/Space Engngq., 
vol. 19, May 1960, pp. 86-87. 


Atmospheric and Space Physics 


Some Results of the Determination of 
the Structural Parameters of the Atmos- 
~_ Using the Third Soviet Artificial 

arth Satellite, by V. V. Mikhnevich, B. 
8. Danilin, A. I. Repnev and V. A. Soko- 
lov, NASA Tech. Transl. F-13, May 1960, 
20 pp. Also ARS Journat, vol. 30, no. 
4, April 1960, Russian Suppl., pp. 407-413. 

Analysis of the Moon’s Surface by 
Nuclear Reactions, by E. F. Martina and 
Carlton D. Schrader, Space Tech. Labs., 
STL/T R-60-0000-00078, April 1960, 31 
pp. 

International Geophysical Year Bibliog- 
raphy of Literature in the Russian Lan- 
guage for 1958, compiled by R. F. Zatrutina 
and Z. M. Mikhaylova, edited by B. I. 
Silkin, NASA Tech. Transl. F-12, May 
1960, 75 pp. 

Some Aspects of Lunar Geology, by J. 
Green, Astronaut. Sci. Rev., vol. 2, no. 1, 
Jan.—March 1960, pp. 9-12, 16. 

Ice in Space, by Donald H. Robey, J. 
Brit. Interplanet. Soc., vol. 17, Jan.-Feb. 
1960, pp. 205-215. 

Energy Dependence of Transient 
Changes in the Primary Cosmic-ray 
Spectrum, by K. G. McCracken, Phys. 
con vol. 117, March 15, 1960, pp. 1570- 
1579. 


Stability of an Idealized Atmosphere, I: 
Discussion of Results, by K. M. Case, 
Phys. Fluids, vol. 3, March-April 1960, pp. 
149-154. 

Stability of an Idealized Atmosphere, II: 
Zeros of the Confluent Hypergeometric 
Function, by Freeman J. Dyson, Phys. 
Fluids, vol, 3, March-April 1960, pp. 155- 
157. 

Artificial Earth Satellites, vols. I & II, 
ed. by I. V. Kurnosova, transl. from Rus- 
sian, Plenum Press, N. Y., 1960, 120, 107 
pp. 

Cosmic-ray Measurements with Arti- 

ficial Earth Satellites, by S. N. Vernov, 

N. L. Grigorov et al., vol. I, pp. 5-9. 

Determination of Atmospheric Density 

from the Observed Deceleration of the 

First Artificial Earth Satellites, by M. L. 

Lidov, vol. I, pp. 10-24. 

Determination of Density of the Upper 

Atmosphere from Secular Variations of 

Orbital Elements of the First Two 

Artificial Earth Satellites, by P. E. 

El’yasberg, vol. I, pp. 25-29. 


Rocket Measurement of Electron Con- 
centration in the Ionosphere by Means 
of an Ultra-short-wave Dispersion In- 
terferometer, by K. I. Gringauz, vol. 
I, pp. 79-84. 

Distribution of Electron Density as a 
Function of Height from Data of Experi- 
ments with Rockets and Artificial Earth 
Satellites: Effect on Propagation of 
Radio Waves, by L. A. Zhekulin, vol. I, 
pp. 85-101. 

Some Results of Measurements of the 
Thermo-dynamic Parameters of the 
Stratosphere Obtained with the Aid of 
Meteorological Rockets, by E. G. 
Shvidkovskii, vol. II, pp. 12-19. 
Disturbances Occurring in the Gaseous 
Medium During the Flight of a Satellite, 
by B. A. Mirtov, vol. II, pp. 20-31. 
Preliminary Results of a Determination 
of the Density of the Atmosphere Above 
100 km, by V. V. Mikhevich, vol. II, pp. 
32-39. 

An Investigation of the Ionic Composi- 
tion of the Earth’s Atmosphere Using 
Rockets and Satellites, by V. G. Isto- 
min, vol. II, pp. 40-44. 

Soviet Ionospheric Studies Using Rock- 
ets and Artificial Earth Satellites, by 
V. I. Krasovskii, vol. II, pp. 45-62. 
Preliminary Report on Geomagnetic 
Measurements Carried out from the 
Third Soviet Artificial Earth Satellite, 
by S. Sh. Dolginov, L. N. Zhuzgov and 
N. V. Pushkov, vol. II, pp. 63-67. 
Micrometeorite Studies Using Rockets 
and Satellites, by O. D. Komissarov, 
T. N. Nazarova et al., vol. II, pp. 68-74. 
The Search for Corpuscles with the 
Help of the Third Artificial Earth Satel- 
lite, by V. I. Krasovskii, Yu. M. Kush- 
nir et al., vol. II, pp. 75-77. 

The Study of the Soft Component of 
Cosmic Rays Outside of the Earth’s 
Atmosphere, by S. N. Vernov, P. V. 
Vakulov et al., vol. II, pp. 78-89. 

The Heavy Nuclei of the Primary Cos- 
mic Radiation, by L. V. Kurnosova and 
L. A. Razorenov, vol. IT, pp. 90-95. 


Human Factors and 
Bioastronautics 


Selection of Personnel for Manned 
Space Stations, by Romney H. Lowry, 
Aero/Space Engng., vol. 19, May 1960, 
pp. 32-33. 

Work and Living Space Requirements 
for Manned Space Stations, by Fred A. 
Payne, Aero/Space Engng., vol. 19, May 
1960, pp. 34-35. 

Gravity Problems in Manned Space 
Stations, by Carl C. Clark and James D. 
Hardy, Aero/Space Engng., vol. 19, May 
1960, pp. 36-39. 

Atmosphere Control, by Carey B. 
Jackson, George R. Roush Jr. and Robert 
M. Bovard, Aero/Space Engng., vol. 19, 
May 1960, pp. 40-41, 94, 96. 

Design of an Operational Ecological 
System, by Eugene B. Konecci and Neal 
E. Wood, Aero/Space Engng., vol. 19, 
May 1960, pp. 42-43, 96. 

Space and Man, by W. M. Meyer, 
Raketentechnik und Raumfahrtforschung, 
vol. 4, no. 2, June 1960, pp. 37-40. (In 
German. ) 

Artificial Earth Satellites, vols. I & II, 
ed. by I. V. Kurnosova, transl. from Rus- 
sian, Plenum Press, N. Y., 1960, 120, 107 pp. 

Research on Animal Flight in an 
Artificial Earth Satellite, by V. N. 
Chernov and V. I. Yakovlev, vol. I, 
pp. 102-120. 


991 


near 
by 
rical 
ative 
ali, 
x 
ictor 
oles 
428, 
vres Le 
8, 
hear 
959 
e in 
L. 
vol. 
959 
ems 
lear 
ma- 
Det. 
of 
20) 
to 
N. 
trol 
vol. 
al 
4 
nole 
vek- 
une 
yes, 
lev, 
. of 
no. 
wer 
ics, 
Re- 
959 
ing 
lec 
and 
no 


INFERNO 


Boevvenuto Cellini (1500-1571), one of the world’s most revered artists, 

was also a pioneer metallurgist. His scientific methods of casting bronze 

and alloys for his statues are nearly as well known as his great artistry. 

But Cellini had his troubles with his furnace, or manica. In his auto- 

biography he says, “‘...all of a sudden an explosion took place, as though a 
thunderbolt had formed and been discharged amongst us...” 

Here at CEA, our scientists and engineers analyze advanced systems. For 
example, we supply missile and ground support equipment manufacturers 
with detail data on the response of their structures to thermal and aero- 
elastic environments. In our studies, we use an invention that would have 
deeply impressed Cellini: our Direct Analog Computer. 


COMPUTER ENGINEERING ASSOCIATES, INC. 
350 North Halstead 
ELgin 5-7121 


Pasadena, California 
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